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CXXVIT. Nuclear Reactions Produced by the Bombardment 
of Boron by Low Energy Alpha-Particles 


By E. S. Suirz, J. R. Wormatp, G. Linpsay-Jones, A. LunpEN 
and A. G. STANLEY 
Cavendish Laboratory, Cambridge* 


[Received July 9, 1953] 


ABSTRACT 


Boron targets have been bombarded by helium nuclei of energy 
between 1 and 2 Mey. Resonances for production of gamma-rays and 
protons from the reactions 1°B(ap)8C*, of protons from 1°B(«, p)8C, 
and of deuterons from 1°B(ad)!2C were observed at 1-13, 1-24, 1-39, 
1-51, 1-64, 1-68 and 1-83 Mev. Resonances for neutrons from the 
*°B(an)18N reaction were detected at 1-51, 1-64 and 2:05 Mev. Angular 
distributions of the particles emitted were measured at 1-13, 1-51, 1-64, 
1-68 and 1-83 Mev and spins and parities of 4-, 3-, 4+, 4- and 4+ 
tentatively assigned to the corresponding !N levels. Spins and parities 
of 3/2— and 5/2+ were confirmed for the second and third excited levels 
of 1#3C. Reaction cross sections and reduced level widths were obtained 
for alpha-particle energies of 1-51 and 1-64™Mev. A resonance in the 
11B(«p)!4C reaction was found at 1-58 Mev alpha-particle energy. 


$1. INTRODUCTION 


WE have studied the reactions which occur when boron is bombarded 
with singly charged helium ions of energy between 1 and 2Mev. 1°B is 
of particular interest as a target for alpha-particle bombardment in 
this energy region because the compound nucleus formed can break up 
in a number of different ways; by emission of neutrons to form #N, 
deuterons to form !2C, or protons to form 18C, either in the ground state 
or in one of the first three excited states. By measuring the angular 
distributions of the emitted particles we have been able to obtain spins 
and parities of certain states of the 44N nucleus and of the four lowest 
states of #C. 

Talbott and Heydenburg (1953) have studied the gamma-rays from 
18Q* produced in the reaction 1°B(«, py)!%C* and have found several 
resonances at alpha-particle energies below 3 Mev. The production of 
protons and deuterons from the reactions !°B(«p)'8C and 1°B(ad)"C has 
been studied for alpha-particles of 7-45 and 6-64 Mev energy by Creagan 
(1949) and Perkin (1950) who reported a number of proton groups 
Ok css BL A Peet Oe v= eer ee 
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corresponding to various levels in the residual ‘°C nucleus. Neutrons 
from 1B(an)!2N were found by Rogers et al. (1952) who used polonium 
alpha-particles, and steps in the thick target yield were observed by 
Heydenburg (see Ajzenberg and Lauritsen 1952a) at H, =1:5, 2°3 and 2-7 
Mev. 


§2. EXPERIMENTAL METHOD AND APPARATUS 


We made all measurements with a beam of singly charged helium ions 
produced by our electrostatic generator. We did not attempt to obtain 
very accurate absolute measurements of the energy of these ions but 
we believe their energy was homogeneous and constant to about | kev 
during any one measurement. The radio-frequency ion-source gave a 
very steady beam of a few microamps of target current. 

We used thin targets of separated boron isotopes about 5 yg cm? 
thick on thick copper or silver backing. They were normally run at 
about 150°c to prevent carbon deposits. 


2.1. Gamma-ray Measurements 


To obtain gamma-ray intensities we used both a sodium-iodide crystal 
with a photo-multiplier, and a Geiger counter. The former was usually 
employed to obtain gamma-ray yield curves, but the latter was more 
stable over long periods of time and was usually used as a monitor 
during experiments on particle angular distributions. We analysed the 
gamma-radiation by using the crystal and multiplier in conjunction with 
a pulse-analyser of the type described by Hutchinson and Scarrott (1951), 
calibrating the system by radiations from radioactive sources and from 
known (p, y) and («, ny) reactions. 


2.2. Charged Particle Measurements 


Charged particles were observed with three kinds of apparatus. To 
detect long-range protons we used a proportional counter, short-range 
particles being stopped by an aluminium foil equivalent to 18 cm of air. 
For protons of shorter range and for deuterons we used either an analysing 
magnet fitted with a zinc-sulphide screen and photo-multiplier, or a 
double ionization chamber. The magnetic analyser had 54° inflexion 
focusing with an energy resolution of about 8°, and a solid angle of 
acceptance of 4 10-4 steradian. 

The ionization chamber, which we filled with pure argon to about one 
atmosphere pressure, is shown in outline in fig. 1. Electrons, produced 
by protons, P, entering the chamber through the window, W, were 
accelerated through the grid, G. The two sets of small pulses arising 
from electrons intercepted by electrodes C, and C, were amplified about 
twenty times and then mixed in opposition. By suitable adjustment 
of the relative amplification of the two sets of pulses we were able to 
arrange that only protons stopping opposite C, gave an output pulse ; 
longer range particles passed both C, and C, producing no output. 
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We used mica foils to alter the range in the chamber of particles to 
be recorded and moved C, and C, parallel to themselves to provide the 
fine adjustment necessary to bring the end of the range opposite the 
gap between C, and C,. We were sometimes able to separate the second 
fastest group of protons from faster and slower groups by connecting the 
electrodes C, and C, effectively in parallel instead of in opposition, or 
by using one electrode only. The slower groups were then stopped 
before reaching C, while the fastest group of protons gave pulses small 
enough for adjustment of the scaler discriminating voltage to exclude them. 
The advantage of this ionization chamber over the magnetic-analyser 
lay mainly in the larger solid angle which could be used, though its energy 
resolution was not as good. 


Fig. 1 
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Differential ionization chamber. For details see text. 


2.3. Neutrons 


To detect neutrons we used a zinc-sulphide-plastic phosphor and a 
photo-multiplier as described by Hornyak (1952), setting the discri- 
minator bias high enough to prevent detection of gamma-rays of energy 
less than 6Mev. We placed this detector in line with the incident 
alpha-particle beam so that it subtended a solid angle of about three 
steradians. To enable us to compare the neutron yield with the yield 
of other particles, R. Edge found the absolute cross section of the 
10B(an)3N reaction for an alpha-particle energy of 1-51 Mev by measuring 
the production of radioactive manganese when the neutrons were absorbed 
in sodium permanganate solution. 


2.4. Angular Distributions 


We measured the angular distributions of charged particles by mounting 
the proportional counter, the ionization chamber or the magnetic-analyser 
so that each could be moved over an angular range of 20° to 160° relative 
to the incident helium-ion beam without shutting down the electrostatio 


generator. We could not as a rule work at angles forward of 90°. Most 
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distributions were measured with only the most suitable detector for that 
distribution, but a few measurements were made both with the magnetic- 
analyser and with the ionization chamber. 


§3. EXPERIMENTAL RESULTS 
3.1. Resonances and Yields 
(i) Gamma-ray Intensities 


Figure 2 (a) shows the yield curve which was obtained with a 1°B 
target for gamma-rays emitted at 90° to the incoming alpha-particles. 


Yield 


10 HS 


20 MeV. 


Alpha-particle energy. 
Yield curves for °B. The vertical scales are different for each particle 


: For 
absolute yields see table 1. 


Resonances occur at 1-13, 1-39, 1:51, 1-64, 1-68 and 1-83 Mev (all 
-++0-01 Mev). The resonance at 1:13 Mev was previously found by 
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G. A. Jones (private communication) and those at 1-51, 1-63, 1-68 and 
1-83 Mev by Talbott and Heydenburg (1953). We found no gamma-rays 
when we repeated the experiments with a separated 1B target. 


(ii) Protons 


With the magnetic-analyser we found four groups of protons from the 
10B(ap)8C reaction, corresponding to Q values of 4:13, 1-00, 0-41 and 
0-27 Mev (all +0-02 Mev). These groups will be referred to as py p, 
Pe Ps respectively. They arise from reactions in which the residual 
SC nucleus is left in its ground state or in one of the first three excited 
states at 3-13, 3-72, 3-86 Mev. Our values are in agreement with those 
quoted by Ajzenberg and Lauritsen (1952 a). We made a careful search at 
the 1:51 and 1-64 Mev resonances for a proton group which has been 
occasionally reported (Ajzenberg and Lauritsen 1952 b), with Q—=3-3 Mev 
corresponding to a level in °C at about 0-7 Mev, but we could not find it. 
If such a group exists its intensity is less than 7% of the intensity of 
the corresponding py group, which is itself weak. 

Figure 2 shows yield curves for the various proton groups obtained by 
counting at 90° to the incident alpha-particles. We used the propor- 
tional counter for the py group, the ionization chamber for the p, group 
and the magnet for the p, and p, groups. 

Resonances for the py group appear at all the gamma-ray resonances 
except 1-39 Mev and there is an additional resonance at 1-24 Mev, where 
the gamma-rays are weak or absent. We could not detect resonances 
for the p, group at 1-13 and 1-24 Mev but observed marked resonances for 
the energies 1:39, 1-51, 1-64, 1-68 and 1-83 Mev. The solid angle of 
acceptance of the magnetic-analyser was too small for us to measure 
the p, yield curve below 1-4 Mev, but we observed resonances at the 1-51, 
1-64, 1:68 and 1-83 Mev. Resonances for the p, group occur at the 
energies which give gamma-ray resonances. This group is much the 
strongest for energies from 1-5 to 2 Mev. 

Similar experiments with a thin 11B target showed only one pronounced 
resonance for the emission of protons at an alpha-particle energy of 
1-58-++0-03 mev. Figure 3 shows the yield curve obtained by means of 
the proportional counter. The energy of the emitted protons, measured 
by means of the magnetic-analyser, corresponds to a Q value of 
0:85-0:02 Mev which is to be compared with the value of 0:75-0:01 
Mev measured by Frye and Wiedenbeck (1951) from experiments using 
polonium alpha-particles. Since the analysis of the particles emitted 
from our 1°B targets showed no sign of this group, we are satisfied that 
our results for 1°B have not been disturbed by traces of ‘1B in the °B 


targets. 


(iii) Deuterons from the 1°B(ad)#C Reaction 


Figure 2(f) shows the yield curve we obtained with the magnet for 
deuterons at 90°. Resonances occur at all the energies given earlier for 
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protons or gamma-rays except at 1-24 Mev. That this sree, i eas 
of deuterons was confirmed by the change in energy when the a y : 
was moved from 90° to 160° or when a thin foil was placed ities ered 
and analyser. The measured Q value of this group 1s eee i. 
in reasonable agreement with the value 1-35 Mev calculated from masse 
(Ajzenberg and Lauritsen 1952 a). 

Fig. 3 
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{5 2:0 MeV. 
Alpha-particle energy. 
Yield curve for |B(ap)MC. 


(iv) Neutrons from the ©B(an)!8N Reaction 
Figure 2 (g) shows the yield curve for neutrons measured at 0° to the 
incident beam. Definite resonances were observed only at 1-51 Mev 


and at 1-64 Mev in the region below 1-9 Mev and there appears to be a 
third resonance at 2-05 Mev. 


(v) Relative Intensities of Different Particle Groups at each Resonance 


At 1-51 Mev alpha-particle energy the intensities of proton and 
deuteron groups were sufficient for us to measure all of them with the 
magnetic-analyser. We then obtained the intensity of each group at 
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other resonances relative to its intensity at 1-51 Mev from the appropriate 
yield curve. Our results are shown in table 1. The relative values for 
the 1:51 Mev resonance are accurate to about 10%, but the errors for 
the weaker groups may be considerably larger. 


Table 1. Number of Particles, per milli-coulomb of He+ Ions Incident 
on a Nominal 5 pg cm~? Target, Emitted at 90° in Laboratory Space 
into a Solid Angle of 4x 10-4 steradian 


Figures enclosed in brackets indicate that a clear resonance was 
not observed. Figures marked with an asterisk were obtained by 
measurement only at the energies given in column 1 and not from 
yield curves. 


(vi) Gamma-ray Energies 

Gamma-rays may be expected from transitions in the “N nucleus 
formed by capture of an alpha-particle or from transitions between the 
four lowest states of 1C (see fig. 4). We could distinguish only one 
gamma-ray energy above 0-5Mev at 3-73+0-06 Mev. Within the 
accuracy of our measurements this may be ascribed to transitions from 
the second or third excited level of °C, or from both, to the ground state. 
We were not able to detect the small amount of 3-1 Mev radiation which 
the existence of the p, group of protons shows must have been present. 
Radiation due to alpha-particle capture must have been weak in compari- 
son with that from the excited *®C nucleus. 

Below 0-5 Mev pulses due to gamma-radiation of energy 0-21-+0-03 
MeV were superimposed on a broad background of small pulses arising 
from the presence of the 3-7 Mev radiation. By making measurements 
with and without a lead sheet 0-11 cm thick in front of the scintillation 
counter it was possible to obtain the relative intensities of the 0-2 Mev 
and 3-7 Mev radiations. It is reasonable to assume that the 0-2 Mev 
radiation arises from transitions between the second and third excited 
levels in 8C. With this assumption, and after correcting our experimental 
results for the different efficiencies of the counter for 0-2 Mev and for 
3-7 Mev radiation, for the measured anisotropy of the gamma-radiation, 
and for the 3-7 Mev radiation following the emission of p, protons from the 
4N nucleus, we estimated that about 30% of the “C nuclei in the 
3-9 Mev level reached the ground state by way of the 3-7 Mev level. 
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3.2. Angular Distributions 


(i) Charged Particles 

We measured the angular distributions of the charged particle groups 
only at resonances that were sufficiently strong relative to the tails of 
adjacent resonances for us to consider the whole yield to be due, at least 
to a first approximation, to a single level of the compound nucleus, 4N. 
We have analysed our results, after correcting to centre-of-mass 
coordinates, in the form 


N, =1-a cos? 6+5 cost 6, 


except for the case of the py group at 1:83 Mev where an additional 
term k cos @ was necessary. 
Table 2 gives the values of a, b and k we have obtained. 


(ii) Neutrons 

The intensity of these at 1-51 Mev and 90° was 3-5 times that at 0°. 
Errors of measurement were large and this figure is consistent with a 
distribution the same as that for p, protons from the mirror reaction. 


§4. COMPARISON OF EXPERIMENTAL RESULTS WITH 
THEORETICAL PREDICTIONS 


4.1. Theoretical Angular Distributions 
If we assume that only one level of the compound nucleus, 44N, is 
involved for a given type of reaction the theoretical angular distribution 
of emitted particles relative to the incident alpha-particles is given by 
N(A)=22| 221+ 1) (9710 | JI.) |? x 
ie 52/1 L 


Geda4 
(FU) CTS, [79 lim? SY" (6,0) | 2 ee 


j is the resultant spin of alpha-particle and target nucleus, 1°B, and is 
thus always 3. J is the total angular momentum of the compound 
nucleus and j’ the combined spin of the residual particles. f,(l) and 
f,(U’) are complex factors describing the probability amplitudes associated 
with orbital momenta, / and 1’, of the incident and emitted particles. 
(jj10| JJ.) and (JJ,|j'j’J’m’) are transformation coefficients (Condon 
and Shortley 1951) and Y,’”"’(@, 0) is a spherical harmonic (cf. Seed and 
French 1952). Since both the proton and the neutron have spins 1/2 
and the corresponding residual nuclei, %C and 18N, spins not zero, j’ may 
in general, though not always, take two values. When this is so, the 
observed angular distribution will be the sum of two functions N(6), 
one multiplied by a factor, the spin-channel factor, which can only be 
predicted by assuming a specific type of coupling between spin and 
orbital momentum. Table 3 gives theoretical values of N(4) calculated 
on the assumption that only the lowest values of J and I’ permitted by 
conservation of momentum and parity take part. 
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4,2. Assignment of Spins and Parities 
Rotblat (1951) has shown from a study of the 12C(dp)!8C reaction that 
the ground state and the excited levels occurring in C at 3-09, 3-68 and 
3°89 Mev above the ground state have spins and _ parities 1/2- ; Ves 
3/2~ or 1/2-; 5/2+ or 3/2+ respectively, while the ground state of 12C is 
known to be 0+. If we compare tables 2 and 3 for the emission of Po 
Fig. 4 


Energy levels. 


Energies are plotted vertically in Mev. Bombarding energies are plotted to 
scale in centre-of-mass coordinates and indicated in laboratory co- 
ordinates. Figures in brackets are taken from Ajzenberg and Lauritsen 


(1952 a). 
protons and deuterons we can tentatively assign spins and parities of 
4-, 3-, 4+, 4-, 4+ to the 14N resonances formed at 1-13, 1-51, 1-64, 1-68 
and 1-83 Mev bombarding energy. For the 1-68 mev level a better fit 
for the deuterons could be obtained by assuming (J, P)=1~ (P=parity), 
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formed by f-wave alpha-particles, but this is unlikely as such a resonance 
level could emit neutrons and p, protons with /’=0, resulting in a broad 
level and prolific emission of protons and neutrons, contrary to our 
observations. Comparison of the theoretical and observed p, distributions 
confirms the above assignments and definitely rules out the assignment 
3/2+ to the excited level at 3-89 Mev in the *C nucleus. The experimental 
p, distributions for the 1-51 and 1-64 resonances fit the theoretical distri- 
butions quite well if a suitable choice of spin-channel factor is made. 
For the 1-68 and 1-83 Mev resonances the fit is better for (J, P)=4— 
and 4+ respectively than for the other values of J listed in table 3. With 
a suitable choice of spin-channel factor the experimental p, distribution 
for the 1-51 Mev resonance agrees with that predicted theoretically for 
(J, P)=3- for the 14N level and 3/2- for the 3-68 Mev level of #C. The p, 
distribution does not rule out 1/2~ for this level of °C though the agree- 
ment between the predicted distribution and our observations is not as 
good as when 3/27 is chosen. 

There is, however, evidence for preferring 3/2~ to 1/2~ in the existence 
of the 0-2 Mev gamma-radiation which we detected (§3.1 (vi) above). 
The transition from the third excited level (5/2+) of ®C at 3-9 Mev to 
the ground state (1/2~) involves 4J=2 and change of parity, so that 
magnetic quadrupole radiation may be expected. If the second excited 
level at 3-7 Mev were 1/2~ the 0-2 Mev radiation arising from transitions 
between the third and second excited states would also be magnetic 
quadrupole, about 10’ times weaker than the radiation of 3-9 Mev energy 
(Blatt and Weisskopf 1952). The 0-2 Mev radiation would have been 
much too weak for us to detect. On the other hand, if the second excited 
level is 3/2~, the 0-2 Mev radiation will be electric dipole, comparable in 
intensity with the 3-9 Mev radiation. The branching ratio we observed 
was 0-3, in good agreement with this assumption. The evidence from the 
0:2 Mev radiation thus strongly supports our assignment 3/2- to the 
second excited level of C. 

Figure 4 gives the energy level scheme relevant to our experiments and 
summarizes our results. 


4.3. Discussion of Assignments Made 


Comparison of tables 2 and 3 shows that for our assignments of spin 
and parity to the N levels the experimental and theoretical values of 
the coefficients a and b of cos?@ and cost @, though always of the right 
sign and right order of magnitude frequently differ by more than the 
experimental error. We think these discrepancies can be accounted for 
by our neglect, in the theory, of values of / and 1’ higher than the minimum 
values permitted by angular momentum and parity conservation and by 
neglect of interference between levels. Interference is clearly shown in 
the case of the py proton distribution produced by 1-83 Mev alpha- 
particles. We have not included in table 3 distributions for a compound 
nucleus state with (J, P)=3+. If such levels had been formed by 


- 
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alpha-particles of J=0 only, all the emitted particles would have had an 
isotropic distribution. We found no level where this was the case. 
If a mixture of =0 and /=2 alpha-particles formed a 3+ level, anisotropic 
distributions would be possible, but we have not been able to fit such 
distributions to our observations. 

Experiments, which will be reported separately, on the angular 
correlation between gamma-rays and emitted particles agree with 
predictions based on our assignments of spins and parities, though the 
processes concerned are too complex for us to deduce spins or parities 
from such correlations alone. 


4.4, Reaction Cross Sections and Level Widths 
From the yield curves, the angular distributions of the emitted particles, 
and the target thicknesses we were able to estimate the reaction cross sec- 
tions in many cases. Owing to the difficulty of estimating the thicknesses 
of the targets used, and their gradual deterioration under bombardment 


Table 4. Cross Sections and Partial Widths 


Cx ips ya" 2M ayn? Vien 2M ayn? 
Ei, (Mev) | Particle|(10-27cm2)| (kev) |(10-#2 kev| 3A? (|(10-1% kev; 3h? 
(x) om) | (%) | em) | (%) 
1-51 ox “= bey 470 10 240 6-0 
Po 4-7 0-62 4-3 0-022 1-7 0-012 
Pp: 1:3 0-17 340 1-7 4] 0-29 
Pe 5:3 0-70 190 0-94 44 0-31 
D: 42 56 140 0-72 66 0-47 
d 7-0 0-93 20 0-28 18 0:26 
n 32 4:3 68 0-34 26 0-19 
1-64 a a= 1-0 620 14 330 8-2 
Po 0-98 0-18 8:2 0-042 1-7 0-012 
Dy 0-46 0-085 | 6100 31 380 2°7 
Pe 2-4 0-44 3400 18 420 3-0 
Ps 52 9-6 2700 14 680 4-9 
d 11 2-0 330 4:5 270 3°9 
n 3:2 0-59 150 0-75 23 0-16 


Columns 5 and 6 give reduced widths assuming reaction radii 3-41, 4-86 
and 5:06x10-1%cm for emission of nucleon, alpha-particle and deuteron. 
Columns 7 and 8 give reduced widths assuming radii 4-86, 5-43 and 
5:15 x 10-8 cm. 


Table 5. Widths of Resonances 


E, (Mev) | 1-13 1-24 1-39 1-51 1-64 1-68 1-83 
I" (kev) 4344 | 3645 | 50+5 | 1444 | 1444) 542 | 2144 


(ie ae Ce ee ee 
Widths have been corrected for target thickness and centre-of-mass motion. 
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our results, which are given in table 4, are subject to considerable 
error, but are, we believe, correct to about a factor of two. We obtained 
total widths I from our yield curves by correcting the observed widths 
for target thickness and centre of gravity movement. Our results are 
shown in table 5. 

In the case of the 1:51 and 1-64 Mev resonances, for which we are 
reasonably confident of our spin and parity assignments, we have estimated 
partial widths I", for the various reactions on the assumption that the 
cross section, o.,, for a reaction can be obtained from the Wigner one-level 
formula 2 oJ+1 Lr, ' 

"= GGT) Gye) +aeF 
where J is the angular momentum of the compound nucleus, s(=0) and 
i(=3) the spins of alpha-particle and !°B nucleus respectively and I’, 
the partial width for alpha-particles. Since [=I,+2T., we must solve 
an equation of the form 


aI ?(28 +1)(21-+1) 
fj BAG RF i) errs ee 210 gs mera * (3) 
This leads to two values of I", and hence to two values of each of the other 
partial widths. In table 4 we have chosen the more plausible smaller 
value of I’, and hence the larger of the two possible values for each of the 
other partial widths. Since we could observe no capture gamma-radiation 
we have put J’,=0 in the above calculations. 

In table 4 we also give the reduced width y,? of a level from the 
relation I,=2kP,y,? and the corresponding value of 6?=2May,?/3h? 
where M and a are the appropriate reduced mass and reaction radius 
(Teichman and Wigner 1952). We obtained P,, for charged particles from 
the tables of Bloch e¢ al. (1951) and P., for neutrons from the formulae 
given by Blatt and Weisskopf (1952). 

In table 4 we have given y,? and 0,? for two sets of values of a; one 
(columns 5 and 6) for a=1-45x 1313 10-8 em for nucleons, 1-45 
(121/34-1-2) x 10-13 em for deuterons and 1-45 x (101/3+ 1-2) x 10-18 em for 
alpha-particles ; the other (columns 7 and 8) for a=1-45 x (A ,1/3+ 4,8) x 
10- cm where A, and A, are the atomic weights (460—16) of the 
particles concerned. 


4.5. Discussion of Level Widths 


As can be seen from table 4 the value of 0? for a given reaction is 
sensitive to the reaction radius, a. For both sets of radii used in table 4 
however, ©? is noticeably smaller for the py protons than for the other 
groups of protons. It is also considerably smaller than 62 for neutrons 
whereas, since py protons and neutrons are produced by mirror reactions, 
much the same value of ©? might have been expected for the two cases. 
We think that the difference between 62? for py protons and for neutrons 
cannot be accounted for by experimental error, though since the relevant 
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cross sections were measured by very different kinds of experiment some 
possibility of such an explanation remains. If we compare our values of 
© for protons with that of 12% calculated by Teichman and Wigner 
(1952) for protons in the ®C(py)!4N reaction at 550 kev bombarding 
energy and nearly the same reaction radius, we see the decrease in 02 
which is to be expected theoretically (Teichman and Wigner, ibid.) with 
increasing excitation energy, 12-7 Mev in our case compared with 8-06 Mev 
for the C(py)!4N reaction. 


§5. CoNcLUSION 

Our results enable us to confirm the assignments J, P=1/2- and 1/2+ 
made by Rotblat (1951) for the ground state and lowest excited state of 
#8C. They enable us to select 5/2+ from his assignment 5/2+ or 3/2+ for 
the level at 3-86 Mev, and 3/2~ rather than 1/2- for the level at 3-72 Mev. 
These values of spins and parities of the four lowest levels in 13C are in 
accord with a simple application of the shell theory. The ground state 
consists of a 1C core (i.e. a complete 1s shell and a complete ps/. sub-shell) 
plus a single p,,. neutron. We can excite the system with the least 
possible energy by promoting the odd neutron to the next shell, a process 
which would require about 4 Mev and give two levels 2s,;. and ds); 
or else we can break the closed ps,. shell and obtain a nucleon, which we 
then pair with the p,/. neutron, a rearrangement which needs the energy 
for lifting a nucleon from ps;. to py/z, plus the difference in pairing energy 
between (p3/.)? and (p4/2)?.. Thus we would expect the following low- 
lying levels: d5/. and sy/y (5/2 and 1/2*), (Ps/2)"*(Pae)” (3/27); Pasa 
(ground state) (1/2~) and we would expect all other levels to be higher. 
These predictions agree with what we have inferred from our experiments. 
With our assignments the three lowest excited levels correspond exactly 
with those of the mirror nucleus N as reported by Jackson and Galonsky 
(1953) from measurement of the elastic scattering of protons by carbon. 
For a more detailed discussion see Inglis (1953). 

It is not possible to draw any very definite conclusions at present from 
our measurements of reduced widths for the various reactions which occur 
when 1°B is bombarded with alpha-particles, but it would appear to be 
well worth while to measure the widths for competing reactions for more 
cases. We intend to do this and also to repeat under better conditions 
our comparison of 6? for py protons and neutrons. 
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ABSTRACT 


On the assumption that the energy per unit length of a dislocation line is 
proportional to the square of the Burgers vector, the stability of dis- 
locations in simple cubic, face-centred cubic, diamond, hexagonal close- 
packed, and body-centred cubic lattices is discussed. It is shown that in 
each of these lattices, there is only a finite number of Burgers vectors 
which correspond to stable perfect dislocations; these vectors are 
enumerated. Imperfect dislocations are discussed for those lattices in 
which they can occur and the finite sets of Burgers vectors corresponding 
to isolated stable imperfect dislocations are found. Stable extended 
dislocations, consisting of arrays of imperfect dislocations separated by 
areas of translation-twin surface, are discussed but there is no limit to the 
number of these that is geometrically possible. 


§ 1. INTRODUCTION 


Tue Burgers vector of a perfect dislocation may be any lattice vector : 
that of an imperfect dislocation may be any lattice vector, added to a 
vector characteristic of the translation-twin surface or surfaces attached to 
the dislocation. Therefore, in principle, there is an infinity of possible 
Burgers vectors for the dislocations of a given lattice. However, if one 
neglects the crystal anisotropy and the small variations of energy with the 
orientation of the Burgers vector relative to the direction of the dislocation 
line, the energy per unit length of a dislocation may be taken as propor- 
tional to the square of the magnitude of its Burgers vector (Nabarro 1952). 
Therefore, if the Burgers vector, b (of magnitude b), of a dislocation is 
equal to the sum of two other possible Burgers vectors, say b, and by, and 
b2>b,2+5,2, ie. b,.b,>0, then the dissociation of the dislocation, 
according to b+b,+b,, lowers the energy. The original dislocation is 
then unstable against such a dissociation and, when this principle is used 


* Communicated by the Authors. 
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consistently, the number of possible Burgers vectors corresponding to 
stable dislocations becomes finite. -Because the dislocation energy is not 
exactly proportional to 6?, dislocations which are stable in one orientation 
may be unstable in another, particularly when b, . b, is small, but such a 
complication will be neglected in this paper. We shall also ignore the fact 
that even a moderate stress may dissociate a weakly stable dislocation 
since a stress will, in general, exert different forces on the two component 
dislocations. 

In the sections which follow, we treat the stability principle as exact and 
apply it to different lattices to find all the possible Burgers vectors which 
correspond to stable dislocations. For each lattice, we first choose three 
primitive vectors, a,, aj, a,, and form all possible combinations, 


b= 1a,+%a.+N3a3 


where the , are integral or zero. These are the possible Burgers vectors 
for perfect dislocations. For each different vector, we then consider all 
possible decompositions into two other vectors, say b, and b,, with 
b=b,+b, and calculate the maximum value of b?—(b,?+0,"). If this is 
negative the perfect dislocation corresponding to b is listed as stable, if 
zero as doubtfully stable and if positive as unstable. The treatment for 
imperfect dislocations is similar and the details will be given as the various 
Cases arise. 


§ 2. SIMPLE CuBICc LATTICE 


In the simple cubic lattice (e.g. caesium chloride), primitive vectors are 
(1, 0, 0), (0, 1, 0), and (0, 0, 1), the unit of length being the lattice constant. 
From these we may construct : 

(a) 6 vectors: (1, 0, 0) etc. (where etc. implies permutations and 
changes of sign), for which b?=1. These are the displacements from a 
lattice point to its nearest neighbours. 

(6) 12 vectors: (1, 1, 0) etc., for which 6?=2. These are the displace- 
ments from a lattice point to its second nearest neighbours. 


(c) 8 vectors: (1, 1, 1) ete., for which b?=3. These are the displace- 
ments from a lattice point to its third nearest neighbours. 


Dislocations of type (a) are stable but the stability of type (b) against 
dissociation into pairs of type (a), or of type (c) against dissociation into 
one each of types (a) and (6) or into three of type (a), is in doubt. Any 
dislocation whose Burgers vector has a component greater than 1 is 
unstable. 

There are no imperfect dislocations in this lattice. 


§ 3. Face-Centrep Cusic Larricr 


In the face-centred cubic lattice, taking the customary unit cell to 
define a Cartesian reference system and using the lattice constant as unit 
of length, we find a set of primitive vectors to be (0, 3, 4), (4, 0, 4), and 


ag: 
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(3, 3,0). From these the Burgers vectors corresponding to stable perfect 
dislocations will be derived in § 3.1. However, as is well known, the face- 
centred cubic lattice can also be regarded as built of close-packed layers of 
atoms in (111) planes. The projections of these layers along the [111] 
axis occupy three different sets of positions and in this way the layers 
may be labelled as a, b, orc. The perfect lattice is then characterized by 
one of the continuing sequences: ...abcabc... or ...cbacba... 
Though many metals crystallize in this system a comparable number do so 
in the hexagonal close-packed system which is characterized by a layer 
structure such as ...ababab... while others, e.g. Ca, La, Tl, Co, and 
some Cu-rich Cu/Si alloys, occur in either of these structures or indeed in 
structures with random arrangements somewhere between these two 
regular sequences. Furthermore, growth and recrystallization twins 
corresponding to the sequence ...abcbacba... are prevalent in most 
face-centred cubic metals. The common feature of all these structures, 
that no two successive layers are alike, implies that every atom has twelve 
nearest neighbours, six in its own layer and three in each of the layers 
immediately above and below. It is customarily assumed that, to a good 
first approximation, the lowest energy states are determined by this 
simple rule alone. Therefore, any surface at which there is an anomaly in 
the stacking without a breach of this rule satisfies the conditions for a 
translation-twin surface (Frank 1951 a) and in face-centred cubic lattices 
there will be imperfect dislocations associated with the edges of such a 
surface. The actual energy of the surface probably arises from the 
reflection of electron waves at the change of lattice pattern but so far no 
reliable estimate of its magnitude is available. 


3.1. Perfect Dislocations 
From the primitive vectors, (0, 4, 4), (4, 0, 4), and (4, $, 0), we can 
construct : 
(a) 12 vectors: (0, 4, $) etc., for which b?=4. These are the displace- 
ments from a lattice point to its nearest neighbours.. 
(b) 6 vectors: (1, 0, 0) etc., for which b?=1. These are the displace- 
ments from a lattice point to its second nearest neighbours. 
Dislocations of type (a) are stable but the stability of a dislocation of 
type (6) against dissociation into a pair of type (a), e.g. 
(I, 0, 0)>(2, op 0)+(3, mat 0) 
isin doubt. All further combinations are unstable. 


3.2. Imperfect and Extended Dislocations 
3.2.1. Isolated imperfect dislocations on one plane 
We may generate a translation-twin on (111) in two ways, viz. (i) by 


removing a layer of atoms, say a, and closing up the gap by the displace- 
ment normal to (111) which brings the c and b layers into contact or (ii) 


4K2 
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by inserting a plane of atoms, say b, between ¢ and a layers which have 
been separated by the necessary normal displacement. These processes 
give rise to sequences (i) ...abcbcabe... and (ii) ...abcbabe... and 
the translation-twins so formed are intrinsic and extrinsic respectively in 
the sense defined by Frank (1951 a). The first involves one breach of the 
stacking rules and the second two, but both involve the same kind of 
anomaly in second nearest neighbour relationships between two pairs of 
layers. We cannot tell at present which has the higher surface energy, nor 
do we know whether one occurs more commonly than the other. These 
forms of stacking fault are illustrated in figs. 1 (a) and (6). 

If a bounded translation-twin is formed, i.e. if only part of a layer is 
removed or added (fig. 1 (c)), then we must have imperfect dislocations 
at the boundaries. These can be classified into four sets according to 
which type of twin is present and according to which boundary is being 
considered, i.e. the left or right hand one when looking at the layers in 


Fig. 1 
é c a———____ a ————. 4 
‘ e 8 NS ee 
é ro) <— a ° c e a 
—~ 
¢ c ce a ———c 
3 r) 5——“"_ a ee 
a a ee 


—s —4, —>)L (3, $s 3)R 
(a) (6) (c) 

Stacking faults in a face-centred cubic lattice forming (a) an intrinsic and 
(6) an extrinsic translation-twin surface. The arrows show where the 
layer has been removed or inserted to generate these structures. A 
bounded faulted region with its associated imperfect dislocations is 
shown in (c). (In figs. 1-6 the [111] direction is vertically upwards in 
the plane of the paper.) 


a standard way. We shall define this standard way by taking the [111] 
direction as ‘top’. The sets will then be denoted by LZ and R for the 
dislocations associated with twins of type (i) and A and p for those associated 
with twins of type (ii). Other translation twins of the extrinsic type 
may be generated by inserting three or more extra (111) planes, but the 
boundaries of these are necessarily groups of dislocations, not single ones, 
so that a consideration of all possible stable dislocations of the four sets 
above is sufficient. 

We shall use the convention that the circuits involved in defining 
Burgers vectors are taken in a clockwise sense. Since the displacement 
of the b-layer in a translation of type (i) above is (—4, —4, —4), typical 
Burgers vectors of dislocations of the L and R sets are (—3, —4, —3)L 
and (3,4, 3)2. Similarly, for the other sets we have (4, 4, 4)A antl 
(—3, —3, —3)p. All other Burgers vectors of these sets can then be 
found by combining these vectors with true lattice translations. Clearly 
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we need consider only one set, say the L set, and the others are then 
easily deducible. For stability, we must first test each L dislocation 
against dissociation into a perfect dislocation and another L dislocation. 
Then further, since the removal of one layer is equivalent to the insertion 
of two, we must test each L dislocation against dissociation into two 
A dislocations, and conversely, since the insertion of one layer is equivalent 
to the removal of two, we must test each A dislocation against dissociation 
into two L dislocations. Obviously, the same argument holds for the 
Rf and p dislocations. 

For the L set, we can construct the following vectors which correspond 
to all dislocations which are not unstable against dissociations involving 
perfect dislocations : 

(a) 3 vectors: (4,3, —4)L, (4, —4, DL, and (—4, 4, 4)L, with 6?=1. 
These represent displacements from a lattice point to the nearest twin 
lattice positions in the (111) plane. 


(6) 1 vector: (—#, —4, —4)L, with b?=4. This is the displacement 
normal to (111) described above. 

(c) 3 vectors : (8, 4 DL, (4, 3 AL and (8, 4, #)L with b=]. 

Of these, types (a) and (b) are stable while type (c) are stable against 
the dissociation 


(3.6 oh > (3,3, +, —2 eh 
provided (4, $, 0) cannot dissociate further. However, if the dislocation 
line lies in the (111) or (111) planes, as well as in the (111) plane, i.e. in 
[O11] or [101], then the (4, 4, 0) dislocation can further dissociate (as 
described below) and the stability of (%, 4, 4)L becomes doubtful. Further, 
dislocations of this type are doubtfully stable against dissociations into 
two A dislocations, one of type (a) and one of type (5), e.g. 
(2, 4 DL > (b, —t —DAH( EDA 

and this is independent of the orientation of the dislocation line. 

Some of these imperfect dislocations have already been described in the 
literature, type (a) being ‘ half-dislocations ’ of the Heidenreich—Shockley 
type and type (b) being the original ‘ sessile ’ dislocations of Frank (1949). 
Type (c) are a new type of sessile dislocation. 

Figure 1 (c) shows L and R dislocations of type (b) while fig. 2 shows the 
arrangement of layers about dislocations of all the L and A types. The 
dissociations of the type (c) dislocations are also shown. 

In order to enumerate all the imperfect dislocations of the above 
types, we must distinguish between the planes which can form the 
translation-twin surface. There are four essentially different ones, viz. 
(111), (111), (111), and (111), and we will denote these by a, 6, c, and d. 
Then any of these imperfect dislocation sets may be described by a symbol 
of the form D, where D denotes one of L, R, i, or p, and i denotes one of 
a, b, c, or d. The dislocations we have specified above are those of the 
set Ly. If we use (h, k, 1) to denote the Burgers vector of a typical 
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n, then it is clear that, corresponding to the set 


(h, k, IL; we have the sets (—h, —k, —))R,, (—A, —k, —l)A,;, and 
(h, k, l)p;. Further, we can easily see that, corresponding to the set 


(—h, k, —l)Dp, 


and 


(h,k,1)Dz, we have the sets (h, —k, —l)D,, 
(—h, —k,1)D,. Table 1 gives the changes necessary to derive the Burgers 
vectors of any set from those of the Lg set. 
Fig. 2 
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Layer structures about dislocations of L and A types in a face-centred cubic 


structure. 


@ indicates the projection of the dislocation line while the 


sloping lines indicate that a change must be made in the nomenclature 
of the plane on passing from the unslipped to the slipped region (see 


Frank 1950). The dissociations of the doubtfully stable dislocations 
with b?=4 are also shown. 


Table 1 
a b c d 
Or) (h, —k, —l) (—h, k, —l) (—h, —k, l) (h, k, l) 
R,A (—A, k, l) (h, —k, l) (h, k, —l) (—h, —k, —l) 


When all twinning planes are considered we have 64 stable and 48 
doubtfully stable imperfect dislocations of this kind. 

If, on the one plane, we have adjoining regions of intrinsic and extrinsic 
translation-twin surface, the boundary line between the regions will be 
the line of an imperfect dislocation. Such dislocations can be considered 
as combinations of two of the imperfect dislocations described above so 
that they can be denoted by a symbol of the form D,D,. The conditions 
for their existence imply that the only occurring combinations are Lp; 
and #,A;. Considering first the Lp, group we find that a typical Burgers 
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os is (3, 3, 3), i.e. a vector of the R, type. Therefore, the new dis- 
location group can be correlated with the known R, group and it is easily 
seen that the stable dislocations of the two groups will have the same 
vectors. However, such a dissociation as 
(pp SU eal + —3, —3)Lat+(—4, & dpa 

(analogous to the doubtfully stable dissociation of (3,4, 4)L above), 
although it involves no change in dislocation energy, does reduce the 
area of translation-twin surface and hence the total energy. Therefore, 
dislocations of this type will be classed as unstable. In a similar way, 
the R,A; group may be correlated with the L, group. Ultimately 
therefore, on any plane we have four stable dislocations of Lp; type and 
four of R,\; type, making in all 32 stable imperfect dislocations of this 
combined type. Strictly, the dissociations L > A-+-A etc., described above 
(and in fig. 2) should have been described as L > A+Lp etc. Such a 
change does not alter the argument. 


3.2.2. Hatended dislocations 


When imperfect dislocations can exist, stable perfect dislocations in 
some orientations become unstable since the dissociation of a perfect 
- dislocation into two imperfect ones can lead to a reduction in energy. 
In such a dissociation a new area of translation-twin surface is created 
between the two imperfect dislocations and the surface energy associated 
with this limits the separation of the component dislocations. Such a 
dissociation is described as an extension and the combination of imperfect 
dislocations and translation-twin surface as an extended dislocation. 

Each dislocation with a Burgers vector of (4, $, 0) type has eight modes 
of extension, two in each octahedral plane. Those in the two planes 
containing its Burgers vector may occur spontaneously, and permit 
glide. Those in the other two octahedral planes will not occur without 
an additional stress as they do not involve a reduction in energy ; since 
one of the component dislocations is sessile, these extensions will anchor 
the dislocation. For the (4, —4, 0) dislocation the extensions are : 


(3, —§ 3)Lat($, —$, —3)Rq and (§, —3, —3)Aat(é, —& 3) Pa 3 
(3, —$, 3)Lo+(s, —8 —3)R, and (§, —% —3)An+(3, —3 3)Po 5 
(4, —4, HLet+(s, —3, —s)R- and (3, —3, —s)Act($, —& B)Pe 3 
(3. —$, d)Lat(d, —3, —3)Rq and (3, —3, —s)at($, —4, b)pa- 


Figure 3 shows the two modes of extension of the (3, —3, 0) dislocation 
in (111). The extensions in (111) would be similar. The first mode is 
that of Shockley and Heidenreich (1948). The energies involved in the 
two modes are presumably slightly different, but a change from one to 
the other must involve an activation energy so that both should occur. 
Neither type is able to climb but both are free to glide in the (111) plane, 
and cause the same slip by their motions. There is no obvious macro- 
‘scopic observation by which they could be distinguished. Figure 4 
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shows the extensions of a (4, 4, 0) dislocation in (111). These are equiva- 
lent to the extensions of (4, —4, 0) in either (111) or (111). 

The perfect dislocations of (1, 0,0) type can be dissociated, with a 
lowering of energy, into two imperfect dislocations, one of (3, $, §) type 
and one of (4, —1, —4) type, in whatever octahedral plane the dislocation 
line may lie. Therefore, all dislocations of (1, 0, 0) type lying in a {111} 
plane become unstable when imperfect dislocations can be formed. 
Since the dislocations of (2, 4, 4) type are sessile, such an extension renders 


the dislocation incapable of glide. 


The extension in (111) is illustrated 


in fig. 5 
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Fig. 4. Possible modes of extension in (il) of a dislocation with a Burgers 
vector of (4, 4, 0). 
Fig. 5. Possible modes of extension in (111) of a dislocation with a Burgers 
vector of (1, 0, 0). 


There are certain groups of more than two imperfect dislocations, in 
which the elastic interactions between members of the group are repulsive, 
but the stacking faults connecting them prevent an unlimited expansion 
of the group, which as a whole constitutes a perfect dislocation of the 
lattice. Such groups can also be classified as extended dislocations : their 
Burgers vectors are necessarily relatively large. Figure 6 shows examples 
of such extended dislocations with resultant Burgers vectors of 
(1, —2, —3), (1, —1, 0), (1, 3, $), and (1, 1, 1). In the first two examples 
the components are glissile, and so the complete groups can glide. The 
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second two examples are sessile. The most likely processes for production 
of extended dislocations are processes of combination of dislocations in 
complex crystal slip. Those in which the components are sessile can only 
attain their equilibrium configuration at temperatures which permit 
vacancy diffusion to occur. Clearly there is no limit to the number 
and complexity of geometrically possible extended dislocations. 


Fig. 6 
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More complex extended dislocations in a face-centred cubic lattice. 


3.2.3. Combinations of two imperfect dislocations on different planes 


If we have two imperfect dislocations on different {111} planes these 
may combine together to form a new type of dislocation, the dislocation 
line of which must lie in the (110) direction common to the two {111} 
planes. Following a suggestion by Nabarro, Thompson (1953) has 
described these as ‘stair-rod’ dislocations. Obviously we have a 
_large number of possible combinations as we have to define one of L, R, 
X, p, and one of a, 6, c, d, for each combining imperfect dislocation. 
We, therefore, have to denote each such dislocation by a symbol of the 
form D;,D;' where D and D’ may or may not be different but 7 and j 
must be different. Since, however, D,D, and D,’D, are indistinguishable 
there are only 96 different sets. Among these there are only two really 
distinct types, any set being simply derivable from one of these. The 
basic two that we shall consider are LL, and L,f,,. 

Treating L,L, first, a vector typical of the combination is (0, 0, 3). 
We now combine this with true lattice vectors and test each resulting 
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vector for stability against decomposition either into a different L,L, 
vector and a true lattice vector or into one J, vector and one L, vector. 
An additional condition for instability in the latter case is that at least 
one of the dislocations corresponding to the component vectors should 
be glissile. Under these conditions, we find the following combinations : 


(a) 1 vector: (0,0, —4)L,L,, having 6?=1/9. 
(b) 4 vectors: (4,0,4)D,L,, (—4,9,a)L,Ly, (0,4, #) Lay, and 
(0, —4, dL,L,, having b?=5/18. 
(c) 2vectors: (4, —4, —4)L,L, and (—4, 4, —4)L,L,, having b?=11/18. 
Types (a) and (b) are stable while type (c) is doubtfully stable against 
the decomposition 
$5 Sg —})L bys (3, 


since the dislocation line must lie in the intersection of the a and 6b planes, 
i.e. along Lp and hence cannot lie in either the (111) or (111) plane, 
ie. the (4, —4,0) dislocation cannot lower its energy by extension. 
All further combinations are unstable. 

For the type L,R,, a vector typical of the combination is (4, —4, 0) 


—i, 0)+ (0, 0, = eee 


Fig. 7 


Possible modes of dissociation of a combined dislocation of LL, type. A single 
line indicates a translation-twin surface of intrinsic type, a double line 
onerd extrinsic type. The positions of the dislocations are indicated 

y dots. 


and, by combination with true lattice vectors and stability tests similar 
to those above, we find the following combinations : 

(a) 1 vector: (—4, 4, 0)L,R,, having 6?=1/18. 

(6) 1 vector: (4, —4, 0)L,R,, having b?=2/9, 
Types (a) and (6) are both stable but all further combinations are unstable. 

However, we can now find further dissociations for the dislocations 
of L,L,, type since, as shown in fig. 7, such dislocations can dissociate 
either es one LA, and one R, dislocation or into one A,L, and one R, 
dislocation. Then diesoain name such as the following can occur : 


it 0, $) Labs <r (3; 4, a) Ry +( ? —+4; 0)A,Ly 
and 


(by 4, —H)LgLiy > (fs —4, Lado + (4, —, — DR 
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so that the only stable dislocation of the L,L, type is one with a Burgers 
vector of (0,0, —4). We can, of course, have the corresponding dis- 
sociations, LR, > R,+A,R,, etc., but these do not lead to instability of 
dislocations having Burgers vectors of the types (a) and (b) above. 


To find the corresponding vectors for the other sets, we use the following 
rules : 


(1) If (2, &, )L,L, and (H, K, L)L,R, are any Burgers vectors belong- 
ing to either stable or unstable dislocations, then the corresponding 
vectors for other intersections are : 


(h, l, k)L,L,, (2, —k, —h)L,L,, (ath k, —h)L La; (—h, =; k) Ly La; and 
(—h, k, —) LL, ; 
(H, L, K)L,R,, (L, —K, —H)L,R,, (—L, K, —H)L,R;, 
Se, Sehy K)L, Ra, and Goats 1G —L)L Ra. 

(2) If (A, k, )D,D,;' is a Burgers vector, then (h, k, 1) is still a possible 
vector when any L in the symbol is replaced by p and/or any R by 4, or 
conversely. Further, (—h, —k, —l) is a possible vector when every L 
or p is replaced by F or A and simultaneously every R or A by L or p. 


Table 2 lists all the vectors corresponding to stable stair-rod disloca- 
tions. It is interesting to note that since the dislocation lines of these 
must lie in the intersection of the two planes concerned, all the stable 
dislocations are of edge type and all are, of course, sessile. 


Table 2 


ab 


LL, Lp, pL, pp \(0,0, —3)|(0, —¥,9)|(—3,9,0)| (3,9, 0) (0,5, 0) (0,0, 3) 


RR, Rx, AR, AA (0, 0, 4) (0,4, 0) (4, 0,0) (—4$,0,0) (0, a5) (0,0, =) 


LR, LA, pk, pa ( & & 9)I( $9, §)|(0, a) (0, $; @) (4,0, &) (% $, 0) 
(3, =; 0) (3; 0, 3) (0, 3; 3) (0, iS 3) ( S 0, 4) ( # et 0) 

RL, Rp, AL, Ap (% —t, 0) (%, 0, =+) (0, b; é) (0, Lp 5) ( ie 0, ') ( 6> é3 0) 
(—4, 4, 0)|(—4, 0, $)|(0, —3,3)| (0.4.4) (3,0, 3) (3,3, 9) 


3.2.4. Combinations of three or more imperfect dislocations 


When we consider combinations of three imperfect dislocations, various 
cases can be distinguished. First, if all three lie in different planes, then 
clearly the dislocation line can have only point size and any analysis of 
stability is meaningless. Secondly, if all three lie in the one plane, they 
can be considered in pairs and no new results are found. However, 
further investigation is needed in the cases where two of the dislocations lie 
in one plane, one belonging to an intrinsic and one to an extrinsic trans- 
lation twin, and the third lies in an intersecting plane. 
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An example of such a combination is L,p,l, for which a typical vector 
is (—4, 4, 0), ie. a vector of the &,L, type. Thus we have the same 
combinations with lattice vectors in both these cases and it can easily 
be shown that any such combination of three imperfect dislocations has 
vectors identical with those of a two-dislocation combination. This 
correlation implies that we need treat only the combinations LiuPaus 
and Lp, which have vectors identical with the combinations #,L, 
and R,R,, respectively. The stability of any three-fold combination 
must, however, be tested against dissociations not only into a perfect 
dislocation and another dislocation of similar type but also, as shown in 
fig. 8, into a simple imperfect dislocation and a different sort of combi- 
nation. For instability in the latter case, we need the further condition 
that the simple imperfect dislocation is glissile. When account is taken 


Fig. 8 
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Pussible modes of dissociation of combined dislocations of Lapgly and Lapalops 
types (orientation of diagrams as in fig. 7). A dotted line indicates a 
plane across which no fault exists after the dissociation. 


of all such dissociations, the stable combinations in the L,p,Z, and 
Lypak, sets are found to be exactly those of the R,L, and RR; sets. 
The vectors corresponding to stable three-dislocation combinations can 
therefore be found by considering the vectors of the corresponding two- 
dislocation combination in table 2. 

Four imperfect dislocations can meet in a line only when we have one 
of intrinsic and one of extrinsic type occurring on each of two planes, 
e.g. Lapalypy- This combination has Burgers vectors equal to those of 
the &,R, combination and in a similar manner to that described above all 
such four-fold combinations can be correlated with two-fold ones so that 
again we need consider only two cases, viz. L,p,l,p, and L,p,F,Ap- 
These can dissociate in such ways as are shown in fig. 8 and when ‘the 
stability tests are applied, we find the corresponding Burgers vectors to 
be those of the stable R,R, and R,L, sets. A particularly interesting 
dissociation is 


(0, 0, —§)Lapaliyp, Te (4, 3 —$)L,+(-34, ny —s)L,+(0, 0, —})PaPo> 
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in which the total dislocation energy is unaltered. However, as shown 
in fig. 8, the area of translation-twin surface is reduced and hence the 
total energy. We may therefore class this combined dislocation as 
unstable. 

For any given pair of planes, there are four possible types of four- 
dislocation combination, so that we have 24 sets in all. The Burgers 
vectors of the stable dislocations can be found from table 2. 


3.3. Use of Thompson’s Notation 

Thompson (1953) has developed a convenient notation for describing 
the Burgers vectors of dislocations in a face-centred cubic lattice. This 
is based on the geometry of a regular tetrahedron having faces parallel 
to the four {111} planes used above. If we construct a tetrahedron 
with the points (0, 4, $), (4, 0,4), (4, 4, 0), and (0, 0, 0), as the vertices 
A, B, C, and D, then the outward drawn normals to the faces opposite 
these points are in the [111], [111], [111], and [111] directions so that 
these faces represent the a, b, c, and d planes respectively. 

The Burgers vectors of the perfect dislocations can then be described 
as follows : 


(a) those of the (0, $, 4) type are represented by the joins DA, etc., 
the choice of letters and order sufficing to cover the twelve cases ; 

(6) those of the (1, 0, 0) type are represented by the symbol AB/CD 
and its permutations. This symbol means either the sum of the vectors 
AC and BD or, geometrically, a vector equal to twice the join of the 
midpoints of AB and CD. Since AB/CD=AB/DC=BA/CD=—CD/AB, 
etc., there are only six distinct permutations for this symbol. 

To describe the imperfect dislocations it is necessary to introduce the 
points «, 8, y, and 6 at the centres of the faces a,b,c, andd. The presence 
of one of these letters in the symbol for a dislocation implies that the 
dislocation includes an imperfect component on the corresponding plane. 

The single imperfect dislocations of L, type are represented by : 

(a) 5A, dB, dC with b?=}4; 

(6) 6D with b?=4 ; 

‘(c) D8/AB, D8/BC, D3/CA with b?=3. 

Dislocations of the R, type are represented by a reversal of the order 
of the letters but, of course, no distinction can be made between L and 
p or R and A so that this notation is not as descriptive as the one used 
in § 3.2 though it is much more compact. For the dislocations on other 


planes, an obvious change of lettering can be made. 
The dissociations of the type (c) dislocations correspond to the forms : 


Ds/AB > DA+6B, 
and D3/AB +> Cd+D6, 
where the geometrical proof of the second needs the use of the formula 


Aé+ Bb+C5=0. 
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The dissociations of perfect dislocations of (4, 3, 0) type follow one of 


the patterns 
AB > Au+aB 


or AB —> A6d+5B, 


the choice depending on which plane contains the dislocation line. For 
the combinations of imperfect dislocations, we have : 


LL, type: (a) «B/CD with b°=1/9 ; 
(b) «B/AC, «B/AD, «B/BC, «B/BD with b?=5/18 ; 
(c) BA/ax, «B/BB with b?=11/18. 

LR, type: (a) «8 with 6?=1/18 ; 
(6) Aa/BB with 6?=2/9. 


The instabilities of the last two groups of the L,L, type result from 
dissociations such as : 


op/AC > DB+aB 
and BA]aa« > Aa«/BB+ AB. 


3.4. Associated Lattices 


The analysis in § 3.1 was based solely on the assumption of a lattice 
having primitive vectors of (0, 4, 4), (4. 0, $), and ($, 4, 0). It can there- | 
fore be applied not only to metals such as aluminium or copper which 
have the simple face-centred cubic structure but also to substances 
which have other atoms in the unit cell and yet still preserve the symmetry 
of the face-centred cubic structure, e.g. NaCl, «-ZnS (zincblende), CaF, 
(fluorite). Further, it can be applied without change to the diamond 
lattice since this corresponds to a face-centred cubic lattice with extra 
atoms displaced by (+, 4, +) from the original set and the two sets of 
atoms are distinguishable by the orientations of their neighbours. There- 
fore such a displacement as (4, }, }) is not a permissible Burgers vector 
for a dislocation. It would in fact produce a configuration in which one 
carbon atom was linked equidistantly either to six others at the corners 
of a flattened octahedron or to four others in a plane. Therefore the 
stable perfect dislocations of all these structures are of the types (a) and 
(6) described in § 3.1. 

However, these complex structures do differ in their twinning properties 
so that, before applying the analysis of § 3.2 to them, we must first 
consider their layer structures. In diamond the layer structure is 
... aabbecaabbee ...and in fluorite it is ...axabBbeycana... If the 
layer doublets and triplets in these structures are manipulated like the 
single layers of the elementary face-centred cubic structure, each carbon 
atom in diamond is still surrounded by four equidistant carbon atoms 
and each ion in fluorite preserves its number of nearest neighbours of 
like and unlike sign. The prevalence of twins in both these substances 
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confirms the possibility of this manipulation. Zincblende with a layer 
structure of ...aabBcyax ... differs from these in that the layer pairs 
are polar in character but their sequence can still be manipulated in the 
same way. A great range of ‘ stacking-fault ’ structures is known for 
carborundum which has the same structure. The regular periodicity 
of these must be attributed to the dislocation systems on which their 
growth is based (Frank 1951 b). 

On the other hand, in NaCl with a layer sequence of .. . aybacBayb... 
no such manipulation is possible. In this sequence any three successive 
letters are unlike in respect of their positions in their respective alphabets 
and Latin and Greek letters alternate. Hence each layer is determined 
by the two preceding ones. Accordingly, twins are not observed and 
there should be no stacking faults. Therefore, all dislocations in a 
NaCl-type lattice must be perfect. 


§ 4. HEXAGONAL CLOSE-PACKED LATTICE 


In the hexagonal close-packed lattice (e.g. magnesium, B—ZnS (wurtzite)) 
the primitive unit cell contains two like atoms whose environments 
differ in orientation. The displacement from one to the other is therefore 
not a permissible Burgers vector for a dislocation line surrounded by good 
crystal (cf. diamond in § 3.4). We shall use a Cartesian reference system 
with the interatomic distance in the basal plane as unit of length. A set 
of primitive vectors is then (1, 0, 0), (—4, 44/3, 0), and (0, 0, K), where 
K is the axial ratio, c/a, which is equal to »/(8/3) for the packing of hard 
spheres. These primitive vectors are parallel to the first, second, and 
fourth axes when the hexagonal structure is referred to the more usual 
four-axed system. 


4.1. Perfect Dislocations 


From the primitive vectors given in § 4, we can construct : 

(a) 6 vectors: (+1, 0, 0) and (+4, +41/3, 0) with 6?=1. These 
are the displacements from one lattice point to its nearest neighbours 
in the basal plane. 

(b) 2 vectors: (0, 0, +K) with b?=K?~8/3. These are the displace- 
ments from one lattice point to its second nearest neighbours. 

(c) 12 vectors: (+1, 0, 4K) and (+4, +4+/3, +X) with 6?=1+ Kk? 
~11/3. These are the displacements from one lattice point to its third 
nearest neighbours. 

Dislocations of types (a) and (b) are stable while those of type (c) are 
of doubtful stability against dissociation into one of type (a) and one of 
type (b). All other combinations of the primitive vectors correspond 
to unstable dislocations. 

It should be noticed that the displacements to the nearest neighbours 
not in the basal plane are excluded as Burgers vectors for the reason 
given in § 4. 
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4.2. Imperfect Dislocations 


Considered as a structure of layers parallel to the basal plane, the 
perfect lattice has the sequence ...abab... where the a- and b-layers 
are close-packed planes which have projections on the basal plane of 
(0, 0)-+-base lattice vectors and (4, §/ 3) + base lattice vectors respectively, 
the perpendicular distance between layers being }K. To obtain imperfect 
dislocations we must consider c-layers which project onto ($, —$+/3)+ base 
lattice vectors, and analysis on the same lines as for the face-centred 
cubic structure (§ 3) indicates that these can be introduced into the 
structure at intervals without greatly affecting the energy. Obviously 
the sequences ...bcbcbe... and acacac... also represent perfect 
close-packed hexagonal lattices. In order to classify the various types 
of translation-twin that can occur, we will use the convention due to 
Frank (1951 b) of describing any of the successions a to b, b to ¢, or ¢ toa 


Fig. 9 
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Possible stacking faults in a hexagonal close-packed lattice. (The [001] 
direction is vertically upwards in the plane of the paper in figs. 9-11.) 


as 4, and their opposites as Y. On this basis, hexagonal close-packing 
produces a sequence ... AVAVA... and any translation-twin will 
introduce a finite number of extra A’s or V’s. Thus a twin that produces 
two extra 4\ successions is denoted by (24), and so on. 

Four different types of intrinsic translation-twins are possible and 
these are denoted by (14), (17), (2A), and (2V) (see fig. 9). Of these, 
the first two involve only one breach of the rule that next nearest neighbour 
layers should be alike while the second pair involve two such breaches. 
This suggests that the former ones may involve a smaller energy increase 
but, on the other hand, the latter ones can be formed by a single act 
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of basal slip and hence they may occur more frequently. Further, 
translation twins of the extrinsic type may be formed by separating 
two layers by 4K and then inserting a c-layer. These are of the types 
(3) and (3\7) shown in fig. 9. As these involve two breaches of the next 
nearest neighbour rule and cannot be formed by one basal slip they may 
be rare. More complex extrinsic twins may be formed but they can all 
be considered as combinations of these simpler types so that our list of 
imperfect dislocations will include only those which bound the six types 
of twin already introduced. 

Since the Burgers vector of the dislocation at the left boundary of a 
translation-twin surface may be equal to any translation of the upper 
portion which will generate the twin and that of the right boundary may 
be the negative of any of these, we can write down possible vectors for 
the bounding dislocations. A typical set is : 


(3, 3V3, SK)L(1A) 


6V3, —$K)R(1A) 


-— 
is) 
v 


—2 —6V3, —3K)L(1V) (4, 6-3, $K)R(LV) 

(=F, =tv/3; 0)L(2A) (3, v3, 0)R(2A) 

b, 4V3, OLY) (—4, —$v3, 0)R(2V) 

(0, 0, $K)L(3A), L(3V) (0, 0, —}K)R(3A), R(3V). 


The most general boundary vector then consists of a combination of the 
appropriate one of these with a true lattice vector and the dislocations 
corresponding to all such combinations must be tested for stability 
against dissociation. Moreover, as in the face-centred cubic case, in 
addition to the dissociation into a perfect dislocation plus another 
imperfect dislocation of the same type, we can also have dissociations 
into two imperfect dislocations of different types. Such possible dis- 
sociations are : 

D(1A)>D(2A)+D(1V) or DBA)+D2V); 

D(2A)>D(1A)+D(1A) or D(8A)+D(1V); 

D(3A)>D(2A)+D(14). 
For these dissociations to occur without delay at least one of the compo- 
nent dislocations must be glissile, i.e. the third index of its Burgers 
vector must be zero. This implies that at least one of the components 
must be of (2/\) or (2V/) type so that the dissociations given for D(2/) 
must be considered further and this is done below. 

On the above conditions, we find the following stable sets (illustrated in 
me. 10) : 

(a) 3 vectors: (+4, —4+/3, 0), (0, $+/3, 0) as L(2A) or R(2V), together 
with their negatives as L(2\/) or R(2A), making 12 in all having b=. 
These are the displacements to the nearest twin lattice positions in the 
basal plane. 

(b) 2 vectors: (0, 0, +4$K) as L(3A), R(3A), L(3V), or R(3V/), making 
8 in all having b2=}K2~3. These are the normal displacements from one 
layer to the next. 
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(c) 6 vectors : 3,6V 38, £$K), (0, —44/3, +4K) as L(1A) or R(1V), 
together with es negatives as L(1\7) or R(1A), making 24 in all having 
b°=5+3K?~1. These are the displacements to nearest neighbour 
Aeoins not in the same basal plane. 


(d) 3 vectors: (0, —$4/3, 0), (+1, $4/3, 0) as L(2A) or R(2V), together 
with their negatives as L(2\/) or R(2A), making 12 kinds in all having 
b2—4/3. : 

(e) 6 vectors: (+4, —41/3, 4K), 44/3, +K) as L(2A) or R(2V), 
together with their negatives as naeN or "RE A), making 24 kinds in all 
having 6?=44+ K?~3, 

(f) 6 vectors: (0, —34/3, +K), (+1, 4/3, +K) as L(2A) or R(2V), 
together with their negatives as L(2\7) or R(2A), making 24 kinds in all 
having b?=4/3+ K2~4. 


Types (a) and (b) are stable while type (c) are doubtfully stable against 
dissociations of the type: 
(3, V3, ¢K)L(LA) > (0, 0, $K)L(BA)+(5, §V/3, 0)L(2V). 
Types (d), (e), and (f) are all doubtfully stable against dissociations into 
one perfect dislocation and one imperfect dislocation of type (a). Further- 
more, the latter two types are unstable against dissociations of the types : 


(3, —¢v/3, K)L(2A) > (3, V3, 2K) L(LA)+(0, —3-/3, $K)L(1A) 


and 


—$/3, K)L(2A) > 2(0, —3-/3, 3K) L(1A), 
but these Pancreriee can te occur by a process involving dislocation 
climb since the component dislocations are sessile (see fig. 10). Since, 
however, this climb involves only a small amount of local diffusion of 
vacancies or interstitial atoms, it should occur quite easily and we will 
class dislocations of types (e) and (f) as unstable. This analysis gives in all 
20 stable and 36 doubtfully stable imperfect dislocations. However, it 
must be noted that, by a rotation of the original cartesian axes through 
60° about [001], we get an equivalent reference system in which / has 
become \ and vice versa. Therefore, the distinction between the A and 
7 sets, though it has been convenient in the above, really depends only on 
the original choice of reference system. In fact, then, there are only 10 
different Burgers vectors corresponding to stable and 18 corresponding to 
doubtfully stable imperfect dislocations in the hexagonal lattice. Unlike 
the face-centred cubic case, there is here only one plane parallel to which 
we can get translation-twin surfaces so that no complexity arises from the 
interaction of two imperfect dislocations on different planes. 

The ordinarily stable perfect dislocations become unstable against 
extension into two or more imperfect dislocations when the dislocation 
line lies in the basal plane. Possible extensions are given below and 
illustrated in il ee 


(1, 0, 0) > ies pee ene A) or (3, ¢V3, O)L(2V) 
ene 


AL2 
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(0, 0, K) > (4, 4/3, EK)L(LA) + (3, —3V3, ZA) R(LA). 

(1, 0, K) > (4, $3, $K)L(IA) +3, —6-V3, A) R(1LA). 
The first of these occurs spontaneously and presumably such extended 
dislocations are the dislocations involved in basal slip in hexagonal lattices. 
The other extensions involve dislocation climb and hence cannot occur 
spontaneously. Further, these dislocations become incapable of glide 
when they are extended. 


Fig. 11 
Oe ee ee (SS 
eg qQ——»—_+_. 2 
a oe 
Ppa oad cso SER eee 
aa a 5 mais sas 3 2 3 
a (24 a a a a 
», 0) 124 (4 5359) Zev) 
(3-5-9) (24) or ACHE 
1 1 x 
+ (4555-0) BRA) + ( — 575-9) Rev) 
= (1, 0, 0) 


+ (= 355-5) R04) + (-& - 555-5) R04) 
= (0, 0, K) =(0, —V3, K) 


We Ne 
or (3, 2/3 > 9) L(1A) 
Fe 
+ (4 -355-5) R04) 
= (1,0, K) 


Possible modes of extension of various perfect dislocations. 


Other more complex stable or metastable arrays may be produced, 
provided only that the imperfect dislocations form the left and right 
boundaries for some kind of translation-twin and that the scalar products 


of the Burgers vectors of adjacent dislocations are positive, e.g. the 
combination 


(0, —3-V3, 0)L(2A)+(3, —3/3, 4K)R(1A)+(—4, —3+/3, $K)R(1A) 
= (0, —/3, K) 


which is shown in fig, 11. 
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§5. Bopy CENTRED Cusio Lattice 
5.1. Perfect Dislocations 


In the body-centred cubic lattice (e.g. sodium or a-iron), taking the 
customary unit cell to define a Cartesian reference system and a unit of 
length, a set of primitive vectors is (—4, 4, 4), (4, —4, 4), and (4, 4, —}). 
From these we can construct : 43° ea 


(a) 8 vectors: (4, 4, 4) ete., with b2=3. These are the displacements 
from a lattice point to its nearest neighbours. 


(5) 6 vectors: (1, 0, 0) etc., with b2=1. These are the displacements 
to the second nearest neighbours. 


Dislocations with Burgers vectors of types (a) and (b) are stable but all 
further combinations are unstable. This is the only case in the cubic 
system where we have found two different types of perfect dislocations 
with definite stability in the lattice. 


5.2. Imperfect Dislocations 


The body-centred cubic lattice can be described as a layer structure by 
taking layers perpendicular to a <{112) direction. There are then six 
distinct layers (distinguished by their projections onto a {112} plane) 
and the perfect lattice is described by the continuing sequence 
... 12345612... Since the twin of this structure is characterized by the 
sequence .. . 65432165... , Cottrell and Bilby (1951) have proposed that 
translation-twins may be formed by sequences in which the stacking 
order is reversed at only one layer, e.g... . 1234345612... However, there 
is no convincing evidence that such a stacking fault has low energy and 
hence that such a translation-twin is a stable configuration. Experi- 
mentally, growth and recrystallization twins of this form are not observed. 
Theoretically, although the stacking fault does preserve the distances 
between atoms in neighbouring layers, it preserves neither the nearest 
neighbour distances, which are of importance if the binding of the lattice 
is due to central forces, nor the ‘ bond-angles ’, which are important if 
directed valency bonds are involved. However, if we suppose such faults 
to exist, then the displacement necessary to form the translation-twin on 
(112) is (—4, —4, 4) and from this we can derive the following Burgers 
vectors for stable imperfect dislocations : 


(a) 2 vectors: (—4, —4, })L and (4, 4, —$)R, having b?=1/12. These 
are the twinning dislocations proposed by Cottrell and Bilby. 

(b) 2 vectors: (4, $, —$)L and (—4, —4, 4)R, having b°=}3. The 
motion of dislocations with these vectors would cause twinning with twice 
as large a shear as motion of dislocations of type (@). 

(c) 6 vectors: (4, 4, #L, (—3, 4, —§)L, (4, —3, —4)L, and their 
negatives R, having b?=2/3. 
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(d) 6 vectors: (—4, —4, —3)L, (—4, % dL, (8, —% aL, and their 
negatives R, having b?=?. 

The symbols L and R here refer to the left and right boundaries of a 
bounded translation-twin when [112] is taken as‘ top’. Figure 12 shows 
the layer arrangement about dislocations with these Burgers vectors. 
Clearly, dislocations of types (a) and (6) are glissile but those of types 
(c) and (d) are sessile. 

We should now consider all possible {112} planes and the possibilities 
of combinations of imperfect dislocations from different planes but this 
will not be done here since the evidence for the existence of a stacking fault 
at all is so inconclusive. 


Fig. 12 
2 ——_»,—___— 6 pe et 2———, 
| ——_-_— 5 (eae See | 
peas eee) Eins sort ocy Z 4 
je igenenney ie a Ne 
j Wim Vy , enter ——— 
' \ } 2 ————a_anaer 


| 


(—$ —% aL or (4, 4, —H)L (3 3, #)L (—4, 4, —)L or (3. —§, —)L 


(—4, —4, —§L (8 —% dL or (—4,% HL 
Layer structures about the stable L dislocations of the body-centred cubic 
lattice with [112] vertically upwards in the plane of the paper. 


If these stacking faults do have a low energy, perfect dislocations having 
dislocation lines in a {112} plane may dissociate to form extended dislo- 
cations which have a lower total energy than the original perfect disloca- 
tion. Ifthe dislocation line lies in (112), those dislocations of (4, 4, 4) type 
with Burgers vectors in (112) may extend in such a way as Ata 


(4, 4, 4) > (4, 4, 3)L-4 (4 $ 5)R, 
and the energy of the component dislocations is then less than that of the 


original perfect one. However, when the Burgers vector does not lie in 
(112) the dissociation would have to be of such a form as 


(3, 3, $) a (3, - s)L+(%, 3 —4)R, 
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in which the total dislocation energy is not reduced. These extensions 
could only occur under an applied stress. However, all dislocations of 
(1, 0, 0) type can dissociate with a reduction in energy, e.g. 


(1, 0, 0) SA (2, aaa, aLb+(s, 4; —$)h. 
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ABSTRACT 


The surface energy of a metal in the normal state is calculated using 
the model of a free electron gas in a finite potential well. A careful 
distinction is drawn between the ‘ lattice surface ’ (i.e. the surface which 
just contains all the cells of the lattice) and the surface of the well. A pre- 
vious calculation (Huntington 1951) for the terms due to electronic 
kinetic energy is repeated in a modified form ; it is shown how sums over 
wave number space can be replaced by integrations if a surface dependent 
density is used. Further a calculation by Huang and Wyllie (1949) is 
criticized in relation to the present results. 

An estimate of the electrostatic energy of the double layer at the surface 
is given. 


§ 1. InTRODUCTION 


Tue total free energy of a piece of metal is the sum of two terms, one 
proportional to the volume, the other to the surface area. The former, or 
bulk energy term, will be given by considering a large specimen in which 
the distribution of electrons is uniform throughout. The latter, or surface 
energy term, requires a detailed picture of the non-uniform distribution 
of electrons near the surface of the metal. This in turn leads to a dist- 
ribution of energy whose sum is in excess of that which would be expected 
if the density of electrons were uniform up to the surface and zero outside 
it. This excess is the surface energy of the electron gas which is the main 
contribution at the absolute zero of temperature. At a non-zero temper- 
ature there will be another contribution due to the surface term in the 
free energy of the lattice vibrations. This arises because the presence of a 
free surface changes the nature of the modes of vibration. In particular, 
modes confined to the surface (Rayleigh Waves) are introduced. 

Our calculation will be based on the Sommerfeld model of a metal. 
Here it must be observed that when surface effects are calculated a clear 
distinction must be drawn between the barrier surface (i.e. position of the 
square cut potential barrier) and the lattice surface (i.e. the surface which 
just contains all the cells of the lattice). In fig. 1, the actual electronic 


* This paper is part of a Manchester University Ph.D. thesis. 
+ Communicated by Prof. L. Rosenfeld. 
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potential, for the case of a flat plate, is plotted as a function of distance 
across the thickness of the plate. The free electron model implies some 
sort of averaging process applied to the periodic potential which determines 
the height of the potential step near the lattice surface. However, the 
actual position of the step, which must be outside the lattice surface, will 
be determined by internal properties. These require that the electronic 
charge density far from the surface must be equal to the positive ionic 
charge density and thus independent of the thickness of the plate. 

The method used to determine the kinetic energy contribution consists 
in summing the kinetic energies of individual electrons over the occupied 
states in wave number space. It will, however, be shown that summations 
over wave number space can be replaced by integrations if the density of 
states is properly chosen. This density includes an additive surface term 
which is a function of position in wave number space. The derivation of 
this density is of importance since it can be used when the volume 
containing the occupied states is more complicated than the sphere which 


Fig. 1 


Polentral Barrier 


Lattice Surface 


The electrostatic periodic potential (full line) which is replaced by the square 
cut potential well (broken line). The positions of the positive ions are 
indicated by crosses. 


is used in Sommerfeld’s model. In particular it will be used in a following 
paper for the calculation of the surface energy for a superconducting state. 

Due to the fact that the electron distribution extends beyond the lattice 
surface, a double layer will be formed at the surface. An approximate 
value of the electrostatic energy of this double layer has also been found 
and gives a small contribution to the surface energy. : 

The first attempt to evaluate the surface energy due to an electron gas 
confined in a finite potential well is made in a paper by Huang and Wyllie 
(1949), where previous work is reviewed. These authors however, do not 
distinguish between the lattice and bounded surfaces and effectively take 
them as coincident. The result differs from the one we obtain. 

Huntington (1951) has calculated the gain in energy when a large block 
is split along a surface. The model of the surface used is the same as ours 
and the same kinetic surface energy contributions are obtained. We have 
nevertheless briefly repeated his calculations for two reasons. Firstly, 
we want to derive the above mentioned relation between integrals and sums 
which will be used in our following paper. Secondly, we wish to comment 
on the discrepancy between Huang and Wyllie’s calculation and the 


present one (cf. Appendix). 
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§ 2. MopEL 


We consider the motion of free electrons in a cubic lattice (spacing a). 
The lattice is bounded by two parallel planes, a distance 2L, apart, and 
perpendicular to one of the lattice axes (z axis). We shall calculate the 
total kinetic energy of the electrons in a rectangular block having dimen- 
sions 20, 2L,,, 2L, in the x, y, z directions respectively, where x and y 
are in the directions of the remaining lattice axes. The origin of coordin- 
ates is taken at the centre of the plate. 

If L,, L, >L,, the conditions on the surfaces z=+L,, y=+L, will not 
affect the results and we shall take the usual periodic boundary conditions. 
For the surfaces parallel to the (x, y) plane the boundary condition will 
be determined by potential barriers, at z=+ #,(£,>L,), where 


¢ —T(\+dL,). ¢. 2) ne 


and d is a small positive distance of order a, while L, >a. 
The number (.4”) of free electrons in the plate is proportional to the 
number of atoms in the lattice, i.e. to the lattice volume (V=8L,,L,,L,) 


N=NV=NU1— dL, -OGL2)) 


where U=8L,L,L, is the volume of the potential well. 
The wave functions for electrons fall into two classes, symmetric or 
antisymmetric with respect to the (x, y) plane namely : 


(1) Symmetric waves 


{ 24, cos k,z exp i(k,z+k,y), |z|<f%, 


b= 4, k 22L (3) 
2A ,(—1)"-D*2 — exp y(L, Fz) expt (k,7+k,y), ; 
| cial) id Pp Y¥(Y,+F2) expt (k,a+kyy) 2<— ZL, 
(2) Antisymmetric waves 
2A, sin k,zexpi(k,«+k,y), |z|<¥, 
p,.= k 27 (4) 
F2A,(—1)"/2 exp y(¥, +2) expt (k,a+h,y), nes 
k z P y( )exp t (k,x+kyy) eels 
The wave number k is given by 
lr mr wn ) 
k., k, k)=(—, => SS 
( a? ky k,) (= > iB, > 2¢, z) (5) 


where sin 6=k,/u ; (0<8<7/2) 1 and m are any integers, and n is a positive 
integer, odd for symmetric and even for antisymmetric waves. Finally, 
the wave number y is defined in terms of the potential barrier height 
(m* is the effective mass of the electron) 


p= Tit) lm* F6 fak a eed oe 
and Yi ks ee Nei, 
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§ 3. SUMMATION OVER THE WAVE NuMBER SPACE 


The possible states can be represented by points in the ky, ky, k, space 
according to (5). Since L,, L, >L,, these points can be poceerd to lie 
on the planes k,=(n7z/2,) a5) distributed with the uniform density 
(L,L,/7*). At the absolute zero of temperature, the occupied states will 
lie in a hemisphere of radius k,,, i.e. on circular sheets of radius (k,,2—h,2)¥/2. 
Now k,,, will be determined by the condition that the total number of 
electrons (4) must equal twice the number of occupied states (2 electrons 
per state due to spin). 


Mn VEatb 
Vee aE wa Kg? — keg™ | et eke ee cae (8) 
where 1,, is the largest integer ean satisfies 
eee ee a bole (9) 
so that in general we must put 
kim =N_tt/2D,+-v7/2L,, eer ees, (LO) 


where v is a small number of 0 (1). 
Using the Euler—Maclaurin formula we can change the sum in (8) to an 
integral 
2L,b ( 


“Nm 
N = ——! on 


(ig? het) dn —Mhy?—h2)ya9+0(F2)b (a1) 


YO 


OES, 
where dn= | ‘ g 


Dh Ce ee WOW) Se (12) 
7 


yf, 


from eqn. (5). Carrying out the integration and substituting for .” from 
(2) and n,, from (10) we find 


ft aw[1—h@-atp+o(c)| -- - a9 
where An =371/8k RM AE Pees ste spa (14) 

= (3/4k,,)[(A—1)1/2+(2—A) sin-2A-¥?]_—w (15) 
and Nee eal ah eee sy (16) 


Now a direct summation shows that near the centre of the plate the density 
of electrons has no term in 1/k,,, (surface term) but is given by 


n(z=0)=(h,,3/372)[1+O(1 bth). . - . « (17) 
Since there cannot exist a free charge in the interior of the plate 
Cie) VCMT Tun ae) re (170) 
so that 
Eta Np Odie, gL {yy = we. (18) 
showing that £,, is surface independent. Thus from (13) 
te ee ed ee Bet ses soar (19) 


— dz is the value of d when p> ©. 
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We can interpret the two main terms of the expression (11) for. 4” as 
follows. The first is an integral over the hemisphere (radius k,,) with the 
wave number density 
U 
sayfa’ —z | the dk, dk, . . (20) 
i.e. p=(V/2n8)[1+(d+1/y/LZ,] . . - . + « (20a) 
so that V/27° is the density for a large specimen. The second term repre- 
sents a correction due to states having k,=0 which are included in the 
integral with weight 4 

Similarly we can salodlate the kinetic energy either by summing and 
using the Euler—-Maclaurin formula or by using the above wave number 
density (eqn. 20) and correcting for zero wave number. This gives 

pene 5a 5 kim k 1 
Fp] ee pig, RT pease Rios > 3 gin-1 2 21 
F=VS|1 [alec ita, SPAS eae “dk, +0 (rar) | (21) 


and substituting for U from (2) and d fae (19) 


p(k, ky, k,) dk, dk dk,= 


Ley? 


FT=VT-+-S (o;+09) <2, 2 2 8 ee 
where 
T=hik,,®/10n?m* 6 ee ee es 
is the bulk kinetic energy density, 
cat.) Ope eS ee eh ke 
the surface area, and 
of", * 16077771." eee ee ee othe sa eked 
=—[hk,,4/1607m*][(A— 1)1/2(14— 15A) + (8— are 15A?) sin-*A3*4)- 
(26) 


The surface energy is made up of two terms, the former (c,) being 
its value for an infinite barrier (A 00) while the latter (c,) is the correction 
for a finite barrier (A finite). 


§4. Tur ELectRostatic ENERGY OF THE DOUBLE LAYER 


Due to the penetration of electronic charge beyond the lattice surface an 
electrostatic double layer is formed at the surface of the metal. Thus there 
exists a non-uniform resultant charge distribution p(Z) near the surface of 
the metal which is negative outside the lattice surface and predominantly 
positive inside. (We take Z=z—,, so that Z=0 gives the position of the 
potential barrier.) The energy of the double layer is given by 


oy=—n{ | p(Z,)p(Z2)|Z,—Z, |Z, dZ_ . . . (27) 


(Huang and Wyllie 1949). 
The electrostatic charge density is given by 


—en(Z) dx dy dZ=—edxdydZ =X \%y | . . . . (28) 
k 
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We will use the free electron wave functions in this sum although the 
result will be used to calculate the electrostatic interaction energy 03. 
This is of course invalid, since some self-consistent solution should be used. 
However, we may expect to get the right order of magnitude even with the 
simple model and since this gives o, only about 10%, of o,, its exact value 
will not be so important. 


Normalizing the electronic wave functions to obtain the amplitude 


(A,) we find that 
l 1 


so that the charge densities are independent of surface area as would be 


expected. Substituting for %, into (28) and changing sums to integrals, 
gives 
p-¢* 5 
n(Z)= 5)" (®—k Eke dk, ZO . . . (30) 
ahs 
kn 


: Chey Lee WE EY 2k,Z4 \dk,, Z<0 
n(Z) = — | (£,,2—k,”) | sin pcg coo 2h,Z— Tee 2k, pee 80), 


v0) 


(31) 


= 
OS 
2 


17) 
g 
(Axa) 3 
(A=3) & 
a 
O- 


= = = - -2 ! 
10 2 S 7 &=2kyZ 


The variation of electron density with distance from the surface. 


For an infinite barrier n(Z)=0 for Z>0 while for 7<0 
Ue 3 cos2k,,Z 3sn2k,,Z 
n(Z)= 53] 1+ a =F | 
The function is exhibited in fig. 2, where the positions of the potential 
barrier (Z=0) and lattice surface are indicated, 
When J is finite, we evaluate the integrals approximately by taking y 
as a constant, different for the two integrals, and determining the two 


appropriate values to give conservation of charge and the correct total 
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number of electrons outside the potential barrier () per unit area) which 
can be determined. Thus 
700 —d 
| n(z) dz= NdZ 20 eel ae ee 
J —L aL 
where L >k,,~1, i.e. a point far rernoved from the surface, and the positive 
charge density of the ions is taken to be uniform. We can rewrite (33) as 


(ie [V—n(Z yaz— | n(Z)dZ+Nd=N,+Nd . . (84) 
where a 
: vg 3N i 
Ny=| 0 n(Z) d= Fer [(3A—2)(A— 1)¥/24(4—3A) A sin-1 A-12], 
(35) 
Thus 
n(Z)=N(2/5A) exp (—a/f) ZO. OS eee 
ae u IKON 6a 24 
1 6 1 24 1 
eae Rierene a, eee BN Oe - (37) 


where «=2k,,Z, f=ky/Yo. J=Kkm/Y;, and yo, y; are the two constants to 
be determined by substituting (36) into (35) and (37) into (34) respectively. 
This gives 
f=(15/16)[(3A—2)(A—1)1/2-+-A(4—3A) sin-1 A-¥/?]_ . (38) 
1/g=(3/16)[(2--A)(A—1)1/2+-A(4—A) sin-1 A-12]. . . (39) 
The approximation is good for large A. Since n (0) which is equal to 
2N/5A does not involve y in the integrands of (30) and (31), the approxi- 
mate solutions will be continuous at Z=0. The gradient is, however, 
discontinuous at Z=0, the ratio of the two values being y;/y)._ In fig. 2 
we also show the approximate solution corresponding to A=3, where the 
approximation is still good. 
Now the charge density to be used in the integral for o, is 
p(Z)= —en(Z) —sikayA | 


=[N—n(Z)]e ie bet) 


We have evaluated o, for the particular case A=3 by a series of 
approximations and find 
og=vNe?; vX5-2. 10-3 
oe 160 fe?m\ v 
(Fe aaa Oe y/ N18, . . : . (41) 
In the case of Au (A=2-97, N=5-93 1022), o,/0;~0-13. 
The calculation showed that v was not very dependent on A so that 


(41) with v5. 10-% will be used as an approximate result for all the 
metals. 


ee oe 3a 


- 
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§ 5. Tue Toran SuRFACE FREE ENERGY 


To estimate the numerical value of o we require the value of 6. Thisisa 
bulk property of the metal and can be obtained from the binding energy 
(W). The free electron wave functions satisfy the equation 


j2 
— Fp V Put Ol) P= ee Pr ~ + 2 +, + (42) 


where we have taken the potential funtion ®(r) to be zero outside the 
metal and —¢ inside it. Thus the total electronic energy is 


h2k? h2k? 
a A ‘) am? is aei(43} 


J 
(The surface term o, is omitted, since the Coulomb interaction is not taken 
into account in eqn. (42).) 


Fig. 3 


The correction term for finite barriers. The experimental points for monovalent 
metals are inside the rectangle ABCD. 


The value of 4 is now obtained by equating the bulk energy to the 
binding energy 
—W=V(T—N¢). ee wale wf (44) 
Table 1 contains values of \ based on values of ¢ as obtained from eqn. 
(44) by Huang and Wyllie. . 
Collecting our results, the surface free energy at the absolute zero o 
temperature is 


a(O)=0,+9,+63 
= 2-22 (N . 10-2)4/[1+ €(A)]+1:1 NV 10-*terg.cm-? . . (45) 


where <(A)=o,/c,, and is plotted in fig. 3. (The normal mass of the 


electron has been used.) 
The change in o due to a finite temperature 7 will in the main be due to 
the thermal excitation of the lattice waves. This contribution has been 
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recalculated by the author (Stratton 1953) and is 


o,=—9A(K) k0)N?2Z (2x) ei le ee CG} 
where “ 
Ly ea 

eas — | a ea So a pag A eae ee 47 

Ue)= 23) say 3 *—Po/ (47) 


and 6, is the Debye temperature. The function Z(x) is tabulated in 
Huang and Wyllie’s paper while the constant A(X) has been evaluated 
for three values of the elastic constant K (K=1/2s—1 where s is Poissons’ 
ratio) by the author. 

Except for an experiment on copper (Udin et al. 1949) all the experi- 
mental values of surface energy refer to liquid metals. We thus evaluate 
a, at the melting point and compare the theoretical and experimental 
values of surface energy (table 1). Since o, is in any case negative we 


Table 1 
95 = —o1(Tm) ce o exp (Tp) 
Metal] Nx10-22| A | (erg| e(A) | (erg | (erg | (erg O1L\” m | (erg ‘ * 
em~*) cem~?) | em~?) |em~-?) (oe cm~?) (erg cma 
Cu 8:5 2-18 | 830 | 0-040] 30 90 950 260 740 1100 
(1400) 
Au 5-93 2-97 | 510 | 0-025| 10 70 | 590 140 450 | 600-1000 
Ag 5-90 2-49 | 510 | 0-033} 20 70 | 600 150 450 800 
Li 4-63 2-09 | 370 | 0-044] 20 50 | 440 40 400 — 
Na 2:56 2-58 | 170 | 0-030 5 30 | 210 20 190 290 
K 18333 3°18 | 70 | 0-022 2 10 80 10 70 200-400 
Rb 1-08 3:43 | 53 | 0-020 1 10 60 10 50 —_ 
Cs 0-87 3-70 | 40 | 0-019 1 10 50 10 40 — 
Hg 4-40 3 350 | 0-025 9 50 410 20 390 460 


(The values of o,,, are taken from Huang and Wyllie’s paper while 
the second experimental value for Cu quoted is that for solid Cu near the 
melting point.) 


cannot expect a very close agreement as o(Q) is already less than the 
experimental result in all cases. Thus to get some idea of the order of 
magnitude involved we take K=1, corresponding to metals obeying 
Couchy’s relation. Thus 


o(T)=2-22 (NW . 10-24)4/3[] 4 €(A)]4-1-1 N 10-21—0-016 6,,Z(x) N 10-21 
eg Cnt tees a ed 
The various quantities are given in table 1 for monovalent metals and 
mercury taking one free electron per atom in each case. Although the 
results are probably only applicable to monovalent metals, mercury is 


included since we will deal with its surface energy in the superconducting 
state. 
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§ 6. FuRTHER REMARKS 


We can find the kinetic energy density ¢(Z) in exactly the same way as 
n(Z). This shows that t(Z)+T as Z+— oo, i.e. far from the surface. The 
behaviour of t(Z) which is similar to that of n(Z), shows that the surface 
energy contributions are really confined to the surface region. 

The calculated values of the surface energy are about two thirds of the 
measured values. This is probably all that can be expected from the use 
of the free electron model, to calculate the predominant termo,. 


It can easily be seen from the results of § 3 that the total kinetic energy 
can be written as a series 


eV lca beat Da iin Oe oe (49) 


where the second term corresponds to the surface energy. Thus when 
L,, becomes comparable with a more terms must be evaluated to give the 
total energy. However, in this case it would be better to treat the 
problem as two-dimensional. 
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(Agree eh ND) LX 
On HvuANG AND WYLLIE’S CALCULATION 


The authors effectively assume d=O at the outset. Thus their value of 
as given by eqn. (13) is 


(eee Iee (dee Dea se . se. (80) 


k 


“Mm? 


Using this value of /,,, as the upper limit in the integrals over wave number 
space and putting d=O in the density p (eqn. (20)) gives their result. 
A surface dependent k,,, however, incompatible with charge conservation 
as may be seen from eqn. (18). 

The authors do not notice this inconsistency because in their calculation 
of the o,-term they effectively use the surface independent value 
(k,,2=372N). This is the origin of the misleading remark that (in the 
interior of the metal) ‘‘ charge is conserved by their method within the 
Fermi—Thomas approximation ” (c.f. part of a letter by Wyllie quoted by 
Huntington (1951)). In actual fact, the discrepancy cannot be removed 
whatever method is used ; the question raised by Huntington about the 
adequacy of the Fermi Thomas method, used to obtain a self-consistent 
solution, is irrelevant in this connection. 
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The consistent way of using the Fermi-Thomas method is as follows. 
If —eV is the potential energy of the electron, we have 
2 
ov Eine (N= nl Zs Zee 
dZ* k 
Now if we rewrite eqns. (36) and (37) as n(Z)=NQ(k,,Z) then, including the 
space dependence of k,,, introduced in the Fermi-Thomas approximation, 
we have 
— REnteV(Z)}2mP? 
= 37h 
where EL, =h?k,,?/2m is the Fermi energy in the interior. The additional 
factor Q makes the solution of (51) very difficult. In view of the smallness 
of co, as estimated in §4 no attempt has been made in this direction. 


n(Z) Q[(3m2N)¥8Z] «Sw. (52) 
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ABSTRACT 


The matrix elements of the electron-lattice interaction, for a plate of 
finite thickness, are calculated using Bethe’s (1933) method. These are 
then used to calculate the energy of the superconducting state as described 
by Frohlich (1950). Previous calculations are extended by including 
terms proportional to the surface area. The technique of replacing sums 
by integrals with surface dependent densities has been used throughout. 


§ 1. INTRODUCTION 


Ir has long been recognized that the critical magnetic field H, which can 
destroy superconductivity in a metal cooled below the transition temper- 
ature, depends on the dimensions of the specimen. (Longitudinal 
magnetic field, cylindrical specimen or flat plate.) This dependence is, 
however, very slight, a notable increase in H, occurring only when one of 
the dimensions is of the order of the penetration depth (AZ 10->— 10-§ cm). 
The increase is due to two surface effects, entering into the free energy 
balance at the transition, which have been discussed by Ginsburg (1945). 
London’s theory gives an electromagnetic energy which is confined to a 
layer at the surface whose thickness is of order A, while we may also expect 
a change in the surface energy of the metal. It should be pointed out that 
the origin of these energies is quite different ; the former is due to the 
penetration of the magnetic field and is a macroscopic effect, the latter is 
due to the finite size of the lattice and is a microscopic effect confined to a 
surface layer whose thickness is of the order of 10-8 cm. 

The following investigation was undertaken in order to get some idea 
of the origin and order of magnitude of the surface energy difference. It 
is based on Frohlich’s (1950, referred to as F,) theory of the superconducting 
state. Frohlich shows that when electron-lattice interactions are con- 
sidered, the normal or Fermi distribution does not lead to a stable state. 
A new distribution is derived leading to a stable state (for sufficiently strong 
electron-lattice coupling) which is identified with the superconducting 
state. We have recalculated the energy of this new state for the particular 


* This paper is part of a Manchester University Ph.D. thesis. 
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case of a flat plate (using the model of I) including terms proportional to 
the surface area. In §2 we give the electronic wave functions and 
lattice waves for a finite lattice which are then used to derive the matrix 
elements of the interaction (§3). These can finally be used to calculate 
the energy of the new state by means of a second order perturbation 
calculation (§4). The method consists essentially in the replacement of 
sums over wave number space by integrals with proper densities which 
have been previously used by the author. 


§ 2. MoDEL 


We will assume that the periodic potential energy of an electron has 
the form 


Via, ¥,2)=V 12) Vay) ee 


Here V’, will have the periodicity of the lattice in the x, y directions while 
V, has the form indicated by fig. 1 in I. (Note that V,——eV here.) 
The wave functions then have the simple form (inside the lattice surface) 


P(X; Ys 2)= (Aj Uy,(2) exp th,2+ A_;, Up, (2) exp(—tk,z))va(x, y) expvk. r. 
(2) 


The first factor, which corresponds to the solution of the separated wave 
equation containing V,, has been discussed by Shockley (1939). Here 
Ux, is a function of z having the lattice periodicity while A;,_ is a constant. 
The permitted values of k, can be determined from the boundary 
conditions. The second factor is the two-dimensional analogue of the 
usual Bloch wave corresponding to V, so that 


A 


k=(k,, k,)=(n,7/L,, n,n/L,) and r=(a,y) . .. . 27 (3) 


where n, and n,, are any integers. The constants are chosen so that the 
uv are normalized per unit cell, i.e. 


[n,* ob up, Of dro—=I A le res 


where dry is the volume element of a unit cell. 


As a rough approximation we will take the values of k, and A x, obtained 
in I for a constant electrostatic potential. Thus 


2 
Ay liz 4d tee ee 
|Ax, P= |A-x, | ae 7.(;+4)| sbitby et aren 


since (4) reduces to |u;,, wt [?>—=1/a%. 

The lattice waves will be replaced by the corresponding elastic waves 
obtained when the discrete lattice is identified with an isotropic elastic 
medium (Stratton 1953, referred to as ¢). The analysis shows that the 
normal modes fall into four classes, three of which consist of mixtures of a 
longitudinal (wave number q,) and a transverse (wave number q,) compon- 
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ent while one consists of purely transverse waves. The most general 
displacement can then be written in the form 


u(7)= *'[Zq(2) exp (iq . f) Gaia ie) exp tq ar) Ga} Be(6) 


pace AES ({a> Vy)> S=unit vector in z ee 


4g Z4=4 Piles (19,2) kexp (—iq,,2)}+ ‘Plex (19,2) exp (—7q,.2)} q./%] 


1p, {exp (191.2) -Fexp (ist) ee —4P Lexplign) exp (—iq.2) 34] 
(7) 


for the three mixed types, and the 6, are time dependent normal coordin- 
ates if the Z are normalized, i.e. 


fo rLz 
29 | Zg* Lig d7=8L gliypy | ZAG Bde Ns he XS) 


(p9 is the density of the medium). The sum is to include terms for the 
four types of vibration, both symmetrical (S) and antisymmetrical (4) and 
all the permitted wave numbers q, (or q,) discussed in S. Quantization of 
the Hamiltonian gives the matrix elements of the normal coordinates : 


AN, A(N 
(9q) 9-1, veal (=): (O*) +1, sf (SX) a gees WR) 


all other elements begin zero. (N, is the quantum number of the vibra- 
tion with wave number q and angular frequency w,,.) 

The normalization condition and the known relation between the 
amplitudes F, and FY, (S eqn. (13)) enables us to calculate these constants 
in terms of the wave number of the mode. We shall only be interested in 
the results for waves of type 3 which consist of the superposition of a 
longitudinal and a transverse bulk wave, whose amplitude is constant over 
the volume of the plate. Here the results show that if we neglect surface 


terms 
; 2K i lq dus] 
Gas | Wa ak be Oe ce a en) 
( ) Sa eee tz my) 
while even if we conclude surface terms 
49 : 4 = : Us 11) 
1) 2 (Q)= Ronee gs ok ( 


In this formulae p(q,,)dq,, is the number of permitted values of ¢,, in the 
range Giz, Viz t+ 4h, for given q. This result is analogous to the one we 
obtained for electron waves (I, eqn. (29)) in that the surface terms in the 
amplitude and the density just cancel out. (K is an elastic constant 
equal to (1/2c)—1 where o is Poisson’s ratio.) 


§ 3. Tae Marerx ELEMENTS ror THE ELzcrron Lattice INTERACTION 
With Bloch (c.f. Bethe 1933) we assume a perturbing potential 
U(r)=—u(r) cradle eae ceer ome eas (12) 
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acting on an electron at r due to the lattice vibrations. The wave 
function for the system of one electron, in the vibrational field of the 
lattice may be expressed as 
iN a k)=¢,(r) II ve (64) . . . . . . (13) 
q 


the product extending over all the wave functions of the lattice waves. 
We require the matrix elements of U(r) with respect to these functions 


J={ P*(k’, N’,) U(r) Yk, Ng) dudydzd@ . . . (14) 
EU fH). 1. oe rs 
where 
J t= b*e (Z. grad V) exp (iq . #) d, dr 
x J by 9q tv, 4, IT f by, ty, ,, es ey 
p+q 
J =J be (2* . grad V) exp (—iq. f) dy dr 
Xf Phe, Oa" tn, 49, I Spi, pup Wy... . (17) 
p=q 


A somewhat lengthy calculation along the lines of Bethe’s article shows 
that the perturbing potential leads to two kinds of changes in the system 
which involve only the longitudinal component of type 3. 
3.1. Absorption 

An electron with wave number k (energy «,) absorbs a quantum of energy 
from the elastic wave with wave number q and changes to a state with 
wave number k’ (energy «,).. The probability for the process is given by 

PA(k’, Ny k, N,)= [Jqt |?2(X) 


16 “AN, V2 s : ag 
=e aa a Ax, An, “A Pa 2C(a)S(a) (?.2(X) é eo epapelat 
where 
{2(X)=2(1—cos Xt/h)/(X/h)?; XSey—e,—hwg . . . . . . (19) 
O(a)= dig et : | Uy Ve* grad Vu; ve EXP (10%) At, scene geen ee 
“14) 2Z 2 
. a sinab, é 
OS 2L, sin aa/2~ mt 
The sum in (18) extends over the four arguments 
ah bdip Pk ss Ge ee 
where the initial and the final state are connected by the relations 
MeN ogols k'=k+qie a eee 


and the condition that symmetrical lattice waves lead to transitions 
between electron states of like symmetry while antisymmetrical lattice 
waves involves states of opposite symmetry. 

C(«) has been chosen so that C(O) reduces to the constant C introduced 


by Bethe to give the strength of the interaction. The function 8 (x) can be 
written 


sin al, 
~~ 


Bie ol 


h a) 
when 5 <1 yA ee eae 


2 
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which has a large maximum at «=0. (We will neglect the maxima 
occurring when « is equal to an integral number of 7. They correspond 
to Peierls’ “‘ Umklappprozesse”’.) Now the range of 40/2 IS Gn,%/2, which 
is about 7/2. Thus near the extremes of the range the approximation 
sing «a/2~aa/2 is bad, but then S(«) is in any case small. We will thus use 
(24) for all «. 

The transition probability is small unless one of the « vanishes, when 

for example (eqns. (5) and (10)) 

~2hNg soe g aS Via ‘ 
where surface ae an been ee pPo=N_M(N,=number 
of atoms per unit volume, M=isotopic mass). 

This agrees with Bethe’s result (cf. eqn. (34) of Bethe’s article) for 
large specimens except for the factor {2[1+(2K+1)q,,/Kq,,]}-1 which 
occurs due to the different amplitudes of the electron and lattice waves 
used in our problem. 

The finite thickness (2) of the lattice thus results in an ‘ uncertainty ’ 
in the condition for the conservation of wave numbers. Thus for q,, 
such that one of the equations «=0 is nearly satisfied 


Proc C28*7(¢) OX). 
Now the ‘ effective width’ of the function S?(«) is of the order 
4a~n/L, (cf. eqn. (24)) 
Sey ANCA ira AW ed he 
In a similar manner the function 2(X) is small unless X 0 leading to 
an ‘ uncertainty ’ in the distribution of frequency 4(X/h) . tl. 


2. Emission 
An electron in the state k emits a quantum of energy which is 
absorbed by the elastic wave with wave number q and changes to a 
state k’. The probability for the wee is given by 


P,(k’, Ny’, k, Na) alesis 


a ae A, 7, (37 | B5(a)C(@) |22(). (26) 
q 
All the quantities have the meaning assigned to them for the previous 


process. 
The calculations further show that all the transverse components give 


zero contributions while the longitudinal surface waves give probabilities 
which are smaller by a factor of O(1/q?L?,) than the ones discussed above. 


§4, THE SURFACE ENERGY oF A METAL IN THE SUPERCONDUCTING 
STATE OF THE ABSOLUTE ZERO 

With Frohlich (F,) we now treat the electron lattice interaction as a 

small quantity so that perturbation theory leads to a second order change 
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in energy (Z). From a formal point of view this energy can be attri- 
buted to the virtual emission and re-absorption of vibrational quanta by 
the electrons. Thus an electron in state k emits a quantum fq,c, and is 
transferred to an intermediate state k’, where 


ke = k— q l . . é e e e e e e (27 ) 
and k,’ can have any value. The probability for the transition is however 
small unless /,’ takes one of the two positive values. 


th’ Shi toi. bolle Cane mn ee 


In the intermediate state the energy of the system consisting of the 
electron in the vibrational field of the lattice is 


Ee thnger 


The electron can then re-absorb the quanta and return to its original 
state. Thus 
Iml A RCT Tf a 
H=—822 2» [Yq | fel he) iv qi gs ee ee 


k k’ qz=0 EK TENGE, 


(Qmi=Value of g, corresponding to the Debye frequency) where the 
intermediate state is defined by k’ and q,, and f, is the probability that the 
state k is occupied. The sums go over all the possible wave numbers 
corresponding to one symmetry with the proviso that symmetrical 
lattice waves connect states of like symmetry and vice versa. This, 
together with the two spin systems which may be treated independently, 
gives the factor 8. Substituting for |J,~ |? (eqn. (26) with N,=0) and 
changing sums to integrals leads to the result 


202K VAN 2 Ta ips Se(1—fe)d 
| eee (ees PEL NEE 
(ay ae) ead i oo dk| — an eee bheyg;, 


toy | ak ak( Ae) | Pe) 2 6%) 


arr} Ee —€n TNC G, 


(where q,=k—k’) correct down to terms of O(1/k,,L,). It should be 


observed that in this result it is convenient to integrate over negative 
values of k,. This is purely formal (only positive values of k, give real 
states) and depends on the fact that we may define all the functions of 
k as even functions in k,. The integration over q,, can be carried out but 
great care must be taken to include all the surface terms that arise. 
A eres for zero wave number states (cf. I, § 3, eqn. (11)) is included 
n (30). 

Neglecting the surface term, the result agrees with the value of # given 
n F, (eqns. (2.11) and (3.1) in F,) except that the value of the matrix 
elements used in F, should be multiplied by a factor 2 (cf. Bethe 1933). 
Following Fréhlich, we can split # into the sum of two terms one of 
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which (H,) leads to a simple small shift in the energy levels while the 
other (#,) can be written as 


AFe,o,7 ( V \2C¢ fitede 
jj Ss m~ 0 , Judd 
: aa a [ae Jak Gap ag 
nr 0.2 
+P Jaw’ jak eee } See AE 
2h) \ hip og?) , 


F=C?/(3e,,Mc,?) ; Ch LINC ee ee es es ek (82) 
and we have used the free electron value of «, equal to h?k?/2m 
(<m=7k,,7/2m, the Fermi energy). The dimensionless constant F, 
introduced by Frohlich, is used in the criterion for superconductivity 
(cf. eqn. (59) below), and (c,/2) is the wave number of an electron moving 
with velocity c,. 

We now consider the change in the energy when the normal Fermi 
distribution is altered by moving the states in a shell of thickness a near 
the surface of the distribution a distance a outwards. Thus if 


where 


Higw! 
a) 
O-3 
Or2 
O-| 0-5 Re 
g 
-O-l 
OO. 2 
-0-3 
Taken from F,. 
Chk, 5 OES 1 Gey tg Ae ee ee AED) 
Geos en oe ene 
where 
4(o)— 1, —khy<x<0, et) — 0, eA, ae (30) 


f(2j=—1, —axr=<0, f(r)=—1, 0<¢<a.f(xz)=0, |2|>e, (36) 
and we take 2/k,, as a small quantity of order o)/k,, i.e. 10-3. Thus 
(k’2—k?)=2k,,(y—2), G2? wAk,,2 sin? 6/2, (ot) 
where @ is the angle between k and k’. Substituting into the first term 
of (31), the energy change will be given by 
f= — ee() r[Jaklaw tore v), G8) 
where 
b(x, y)= 20,2 sin?S/[(y—a)*—o,? sin? 6/2]. Mew et (30) 
The integration has been carried out by Frohlich in terms of a function 
P(é), (cf. fig. 1) where é=a?/o? and c is the largest value of op sin (6/2). 
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Thus ieee 
4fe 2k oO 
ee a eee ae S) ee erg ee 
P= —V ene TS PO) (40) 
where 
Oe Ent ae Oh. : sly. if By. oh nie 
Ty 2h thy Omics m? aan q L 


For the Born modification of the Debye theory, 1.€. ¢:=mi=Um 

say, and neglecting surface effects, these two values are 
o%/o,2=1/4v, 4v>1. o/o°=1, 4v<l, Ja pe oe 

where v=N/N, and N is the number of electrons per unit volume. 

However, the theory in S requires a Debye maximum frequency 
(w,,/27) to cut off the infinite spectrum (i.e. %.~>Qmi). Now we can 
rewrite the formulae for M(w,,), the total number of excited modes 
(S, eqn. (68)) as 


¢ Gee 3b | dnt 3b 9 
M(og)=P| 2 6x? (1+ nae Bt gaat) |e 


b(K)=A(K)(1872)19[(2+ (K/[1+2K ])9?]/[2+ (K/[1+2K])], (43) 
and A(K) is a constant given in 8 for K=0, 1, o. 
In analogy to the Born modification we now assume formally that 


where 


d.=Fa=In od 
in eqn. (42). Thus 
Im° 3b = 
M (Qm)=3N = Qn2 1. gee > So Ae (45) 
Im=Im(1—5/dml,) where 9,°=67?Na,. . . . (46) 


giving the dependence of the maximum wave number on the surface area. 
Using (46) we can express the surface dependence of o for the case 
4y>1 (to which we confine ourselves in what follows) as 


a od l 1 35 > a 1 3b 
cia 4y wigela oc : a pe a 


.o==o(1—b/g0,), © af 
where o >o as L,— o<. Similarly 


é=a%/o%= (1+ 2b/g,.D,) where eaatlo® Ames oe tot 


It should be observed that & itself may still be surface dependent if the 
equilibrium value of a is. However we will show that this latter depen- 
dence does not affect the surface energies. Substituting into (40) for é 
from (48) gives 


Y,=—V 


. 


3mc,2N | F = 


21/32/83 


F oe = 
pM are GEG bow cs at) 
say, %s=V(S,+8,/L,). 
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We now consider the change in the second term of eqn. (31). Let 

bene =k,2—k,8 3 t=k—Km ; Y=h' —hus C0 (50) 
and § be the angle between k and k’. 
The contribution to the change from the second term is thus given by 
V 4Fe,,f 1\2 wf ; 

—s VIL TZ sy" (ga) <3) | Moe +VUe #) ak’ dle (61) 
The multiple integral can only be partly evaluated in terms of known 
functions. The remaining numerical integration becomes rather tedious 


because the limits of integration become functions of the next variable 
over which integration is to take place. Writing the result as 


4F en. ( Wy2\ of Riv, &) ; 
we have computed the one value 
BLO ee ee Gs) 


The value of € chosen corresponds to the minimum of —P(é). 
Finally we must consider the change in the zero order (kinetic energy). 
This is 
Oa ia eee 
aD) ee a Ein 26 . 54. 
Fo= Be) Sapo | [+ | folk bere de, . (54) 


where p(k) is the eles of states in the range k, k+dk and can be 
deduced from the formula giving the total number of states having & 
less than a given value (cf. I, eqn. i 


d a 1 7 i tes i, ee 7 
p(k) dk= dks |v 5a 5 kL, += iL, | k,sin wdk) at =) 6 ee ° 
The second term is the zero wave number aie 


oy) _ k 5 
-p(t)= | 1 an +i. sin =| eee 0D) 


Substituting into (54), taking p(k,,+2)=p(k,,) and the value of U from I, 
eqns. (2), vee (15) and (19), we have 


Fiat, OAR otontt he ree a; mh | 
Shes ye (= =) 2 gti at z— sin-t “42 5 (A= 1) 
ep eae ee ae ee oe ©. (88) 


Collecting our results the change in the total energy is given by (eqns. 
(49), (52) and (56)) 
VU(S\+83 ees at ger ead 
ae 7 | — 5 (4) 9P (8) +28 


2 are {—F(4»)!9(R(v, 2 {2P(E)—EP (8) 


=V 


+Flr-+(3—2) sin“? 12—3—1)"4) } ara t(57) 
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The minimum of c%(a) with respect to variation of a is given by 
of'(a)=V(S'(a)+8'(a)/L,)=0. . . . « ~ (58) 
Thus for large L,, S’=0, which corresponds to the minimum found by 
Frohlich occurring for a value of é = a?/c? 0-3. The condition for the 
minimum to occur is 
(F/2)(4v)18>2 ; (4ys>1) >. 0 ae eves 


and should be satisfied by the superconducting metals (cf. table of 
values in F;). 
For finite L, we write a in the form 


a=a(1+<a/Z,). 
Substituting into (58) gives 
a==—s’(a)/S"(a): : 3a) oe 


However, we can use the asymptotic value a for a even in the bulk 
term since 


F(a) = V| s(@)+ 2 s'(ai)+ | 
— V[S(a)4:s(a)/L 1 


If we identify the new state with the superconducting state, then 


—S mip VLAF + AfiL1=v| =O (14 fy], (62) 


where 4/=difference in free energy per unit volume between the normal 
and superconducting phases ; 4f—difference in surface free energy per 
unit area ; H critical magnetic field for a large specimen (L, > «) 
and 


R= Af/AF=2Af/Hoy?. . . . . (63) 
For the purposes of calculating numerical values we take P(£)~ é 


when £<£20-3 and v=1 in (57). Further we use the value of b(K) 
corresponding to K=1(b(1)=0-56) and neglect the small terms in A. 


aoe 3mc?2N 
= | 8 5-2 |os, « ea Pent tae ae Oe 


> 
3mc,7N a 
12-8 [ = 11-19 |, . a aoe GD) 


w= 1 (4)1/8(F/2)1-1—1-9 
i, (4)2°0F]2)—2]0-8 ce NO) 


For mercury the values found are 


A yy?/2=6-95 X 10° erg em-%, Aoy=1600 gauss. 
Af=0-03 erg cm-?, #=5 xX 10-8 cm. 
. . (67) 
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§ 5. EXPERIMENTAL RESULTS AND DISCUSSION 


We have neglected the change in the double layer energy discussed 
in I. There it was shown that in the normal state this energy 


ers Cea oak es Se ie ae eee 52 (OS) 
Thus for the superconducting state the energy will be 


‘ a ee 
Eee a+d3 oo |+[ Cs— O,:0— om (ate as |r. 
where dashes denote differentiation with respect to k (writing o, oc 1/k8 
in general) and the derivatives are evaluated at k=k,. Thus the 
difference (4c) between the two energies is roughly 
£3 P38 (=) =3é(0?/k,2)10-3Z. =. | |. (69) 
3 k 
This is negligible as may be seen by examining the tabulated values of o3. 
The integrations over the shell will only be valid if the k spacing 
between the wave number sheets in k space is much smaller than the 
thickness of the shell, i.e. 
aera Wop lO emy eles. © 2 C70) 
[This a which is the thickness of the shell in wave number space must 
not be confused with the lattice spacing which has a similar symbol 
but reciprocal dimension. | 
We may thus expect changed values of the surface energy when 
L,510- em, i.e. when it becomes of the order of the penetration depth 
(A). It seems probable that if we included non-symmetrical distributions 
in our variational problem a lower minimum energy could be found for 
these small thicknesses. Thus our calculation shows that as regards 
the surface term electrons on the surface of the normal distribution with 
large k, values tend to be expelled from the surface while those with 
small k,do not. The bulk term however shows that all electrons on the 
surface tend to be expelled (cf. the discussion in F, with regard to the 
function L(y)). Thus a shell, thicker for large k, values and thinner 
for small k, values, may lead to a better approximation for very thin films. 
Unfortunately, the increase in critical magnetic field (/7,) only becomes 
measurable when the thickness is of the order of the penetration depth. 
In fact Ginsburg (1945) shows that 
Ae a aes +p. L, (71) 
(ae) = 1—(A/L,) tanh (L,/A) ‘ ; ; 
where H x is the value of H, as L,—> oo ie. for a massive specimen. 
Further, when L,>A, uw 
This analysis has been applied to the extensive experiments of Appleyard 
et al. (1939), on thin films of mercury. Extrapolating to the zero of 
temperature and using (72) gives \’ 1-2. 10% cm for the thickest 
film (L,=5-44.10-5cm). Using the value 4 4:3.10-° cm given by 


m 
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Lauerman and Shoenberg (1949) leaves » Y7:9.10~®. (Using the exact 
expression for H, and the above value of A gives a 8B. LOeS Ging 
We calculated p at various temperatures and thicknesses, finding similar 
curves of . against L, at constant 7 and p against 7’ at constant Lie 
The range of p is from 5.10~® to 3.10~> em, but the variation with L, is not 
monotone. Unsuccessful attempts were made to reproduce the experi- 
mental H,/L, curves with constant » and a different (constant) A. 

These experimental values of » are larger than the theoretical values 
by a factor 100. However it must be remembered that the theoretical 
values refer to thicker plates only. 

Further the theoretical value of the bulk term H,,?/2 is too large by 
a factor of about 16, so that Af is ony about } of the experimental value. 
In any case it is unlikely that A itself will remain constant as L, becomes 
of the order of 10~° cm. 

While we have not been able to give a very definite answer about the 
importance of the surface free energy in thin film experiments, we have 
been able to show on what factors it depends. This should help to give 
a more definite quantitative value when the theory reaches a greater 
degree of refinement in regard to both bulk and thin specimen. 
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SUMMARY 


The angular distribution of long-range alpha particles from the 
reaction N(p, «) has been investigated over the range 230 kev-960 kev 
proton energy. It is found that the levels in 1*O formed by resonant 
capture of protons of energy 340 kev and 1050 kev have spin and parity 
(0, +) and (1, —) respectively. Interference effects in the angular 
distribution caused by these two levels have been observed, but the 
phase differences necessary to account for these effects are not in full 
agreement with those predicted from modern dispersion theory. 


§1. LyTRopUCTION 


THE reaction *N(p, «) was first observed by Burcham and Smith (1939). 
Freeman (1950) measured the Q of the reaction as 4-96-L0-05 Mev. 
Cochrane and Hester (1949) studied the excitation curve by a photo- 
graphic plate technique and reported a resonance at about 400 kev. 
A study by these workers of angular distributions at 445 and 380 kev 
did not reveal any anisotropy greater than about 10%. Mileikowsky 
and Pauli (1952) studied the excitation curve at 90° with protons ranging 
in energy from 200 kev to 940 kev and observed a resonance about 
140 kev wide at 340 kev. 

Schardt et al. (1952) studied the excitation curve of the reaction 
15N(p, «) for proton energies up to 1600 kev. In two overlapping sets 
’ of measurements the alpha particles were observed (a) at 90° for proton 
energies between 273 and 830 kev, (b) at 138° for proton energies from 
650 kev upwards. Resonances for long-range alpha particle production 
were found to occur at 360 kev, 1020 kev and 1210 kev, of width 94, 
150 and 22-5 kev respectively. (Because of penetrability effects the 
resonance at 360 kev corresponds to a level at 340 kev.) In the work of 
Schardt et al. the N(p, cy) and “N(p, y) processes were also studied 
and appeared to have certain resonances in common with 'N(p, «) ; 
these observations are relevant to the interpretation of the present 
experiments and will be discussed in detail later. Absolute cross sections 
for all three reactions were determined by Schardt e¢ al. and an accuracy 
of +15°% was claimed ; their figures have been used in our work. 


* Communicated by the Authors. 
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The present paper describes an investigation of the angular 
distribution of long-range alpha particles for incident proton energies 
between 230 kev and 960 kev with particular attention to interference 
effects taking place between the levels formed at 340 and 1050 kev. 


$2. EXPERIMENTAL TECHNIQUE AND PROCEDURE 


Nitrogen-15 targets were deposited on a tantalum sheet in the 
electromagnetic separator of the Atomic Energy Research Establishment, 
Harwell. Proportional counters were used to detect the long-range 
alpha particles produced when the targets were bombarded with protons. 
The experimental technique is fully described in an earlier paper (Cohen 
1953). 

vite range of angles permitted by the geometry of the apparatus was 
from 30° to 160° and this complete range could be explored up to 930 kev. 
Above this energy the range of angles was progressively restricted by 
scattered protons and useful observations became impossible above 
960 kev. 

Angular distributions were taken at a series of proton energies 
(mostly at 50 kev intervals) and a study was also made of the energy- 
dependence of [(150°)/Z(90°) (the ratio of the intensity at 150° to that 
at 90°). The variation with energy of the angular distribution was 
thus related to the observed excitation curve. 


§3. RESULTS 


In fig. 1 we show typical angular distributions (reduced to the centre 
of mass system) at 380 kev and 930 kev. The former, taken about 
20 kev above the lower resonance, should correspond roughly to the 
angular distribution obtained by Cochrane and Hester at 412 kev, which 
does not depart very significantly from isotropy. The 930 kev curve 
implies a very strong interference effect in the angular distribution at 
this energy. The results at all intermediate energies exhibit this 
interference to a lesser extent. 

Each angular distribution was fitted twice assuming it to be 
alternatively linear or parabolic in cos 6. The variance ratio test was 
then applied to find if the improvement in fit by taking the parabolic 
term was at all significant. The variance of the experimental points from 
the curves fitted to them was compared, using the y? test, with the 
variance expected from the finite numbers of counts, and in no case was 
there evidence that other random errors were significant. 

Below 380 kev, angular distributions were isotropic or contained only 
a few per cent of cos 6 term. Above this energy strong negative cos 0 
terms were found, up to about 0-70 of the isotropic component. Above 
730 kev, strongly significant negative cos? @ terms were observed, rising 
to about —0-30 at the highest energies, for which the level of significance 
of the cos *@ terms is better than 0:1°%. The y? test indicated in each case 
that there was no reason to invoke higher powers of cos 0. 
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tn figs. 2 (a), (b) and (c), we show respectively J(150°)/I(90°), the 
yield at 90°, and the observed coefficient of cos? 6, as functions of proton 
energy. The beam energy definition was +5 kev at 400 kev and 
+11 kev at 900 kev. The excitation curve does not differ significantly 
from that of Schardt e¢ al., though unfortunately our technique did not 
permit us to go as high as the resonant energy near 1050 kev. The 


Fig. 1 
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Some observed angular distributions. 


(a) 380 kev: no anisotropy above the 5% level of significance. 
(b) 930 kev : (140-01) —(0-64 +0-025) eos @—(0-26 +0-04) cos? 6. 
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Fig. 2 
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(a) Variation of 7(150°)/7(90°) with proton energy. Triangles represent points 
obtained from complete angular distributions. The curves, fitted as 
explained in § 5.2, enclose a band showing the standard deviation in the 
fit. That branch marked (i) is obtained with 6=15° ; that marked (ii) 
with ¢=121°. 

(>) Excitation curve at 90°. 
(c) The coefficient of cos? @ in the angular distributions plotted as a 


function of proton energy. The curves are fitted as explained in § 5.2 
and show the standard deviation in the fit. 
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growth of the cos? @ term with increasing energy is clearly shown, and 
its association in some way with the 1050 kev resonance is strongly 
indicated. The I(150°)/I(90°) curve shows the change in the character 
of the angular distribution with energy. The full curves in figs. 2 (a) 
and 2 (c) are theoretical fits which will be discussed later. 

It is evident without any formal analysis that in this reaction there 
is strong interference between two states of opposite parity. In the 
following sections we discuss the situation in detail, with the principal 
aim of establishing the spins and parities of these states. 


$4. THEORETICAL ANGULAR DISTRIBUTIONS 
4.1. Single Levels 
The reaction we are considering is : 
DPN —_> 160 —_» BO. 

The ground state of °N has spin 4 (Kriiger 1938, and Wood and Dieche 
1938) and the shell-model predicts odd parity (Hornyak et al. 1950). 
The ground state of #?C has spin O and any reasonable model predicts even 
parity; the proton and the alpha particle are (4, +) and (0, +) respec- 
tively. Thus the spin state of the initial system can be (0, —) or (1, —), 
and that of the final system is (0, +). 

The selection rules for angular momentum and parity demand that 
those states in the compound nucleus '*O which emit long-range alpha 
particles (leaving !2C in its ground state) must be (0, +), (1, —), (2, +), 
etc. Moreover, none of these states can be formed from an initial spin 
configuration (0, —); the 1°O states that concern us must therefore be 
formed from the initial configuration (1, —). The state (0, +) is formed 
by p-protons, the state (1, —) by s- and d-protons (2, +) by p- and 
f-protons, and so on. 

The emission of alpha particles from a compound state (0, +) is of 
course isotropic ; the angular distribution expected from a (1, —) state 
in 160 is 

I(0)~1+ 2c(po/p9)/2(1—3 cos? 6) cos y--c?(po/Mo)(14-8 cos? 6), . (1) 
where the isotropic term arises from formation of the compound nucleus 
by s-protons, the term in (1+3 cos? @) from d-protons, and the term in 
(1—3 cos? @) from an interference between s- and d-protons. c? represents 
the probability of formation of 1*O through d-protons relative to that 
through s-protons, y is the phase difference between the s- and d-proton 
wave functions at the nuclear surface, and p, the barrier penetrability 
for protons of orbital momentum /. As p, is much less than py at the 
energies considered, we might expect the first term to be predominant, 
and the third term probably negligible. 

Alpha particles from states (2, +) and (3, —) would be expected to 
have angular distributions predominantly of the forms (1-3 cos? @) and 
(1—2 cos? 0-+5 cos! 6) respectively, if one makes the plausible assumption 
that their formation is mainly due to protons of the lowest possible orbital 


momentum. 
4N2 
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The angular distribution near the lower resonance is isotropic ; this 
level is therefore either (0, +) or (1, —). One may argue that the 
1050 kev level is also either (0, +-) or (1, —), since even at energies where 
the yield might be expected to be predominantly due to it, no positive 
terms in cos? 6 or cos! 6 were found in the angular distributions, although 
there are negative cos?@ terms of about —0-30. This makes it very 
unlikely that the level at 1050 kev is (2, +) or (3, —), and considerations 
of yield and width exclude spins higher than 3. Since the appearance of 
a strong cos@ term at intermediate energies makes it quite certain that 
the levels are of opposite parity, it is very probable that one of the levels 
is (0, +) and the other (1, —), and we shall next consider the angular 
distributions arising from the interaction of two states with these values 
of spin and parity. 

4.2. Interference Effects 

The angular distributions expected from overlapping levels have been 
discussed more fully elsewhere (Devons and Hine 1949, Thomson e¢ al. 
1952, Cohen 1953) and we shall merely quote the relevant result. Let 
r,, 0,, refer to the amplitude and phase of the complex Breit-Wigner 
denominator of the level (0,-+) and let Ae’* represent the complex 
amplitude of its inherent probability of formation. Let the quantities 
19, 59, Be'® be the corresponding quantities for the level (1, —) as formed 
by s-protons. Then if we put (B/A)=z, (r,/r.)=r, (8—a)=¢, and 
(5,—5,)=6, the expected (energy-dependent) angular distribution is 

1(8)~1 + (polpi)@*r*{ 1+ 2(po/po)*/?e(1—8 cos? 8) cos y 

+(Po/po)c?(1+-3 cos? 8)} 
—2(polpy)*?xr{ cos ($—8)—2(po/po)/7¢ cos (S—3+-y)} cosd. . (2) 

Here p, c, and y have the same meaning as in eqn. (1). c, y, # and ¢ 
are expected to be nearly independent of energy. 

In the next section we shall make use of eqn. (2) to interpret the 
experimental results. 


§5. COMPARISON OF THEORY WITH EXPERIMENT 
5.1. Identification of the Levels 


The appearance of negative cos? @ terms at the higher energies makes 
it very likely that the (1, —) spin assignment should be made to the 
1050 kev level. The only other plausible explanation would be that the 
cos’? term arises from interference between the 1050 kev level (assumed 
to be (0, +)) and the 1210 kev level (assumed to be of even parity and 
non-zero spin). The relative intensities of the levels as found by 
Schardt et al. (loc. cit.) would give rise to a cos?@ term of about 0-20 
if the 1210 key level had spin 2. However, the work of Kraus et al. (1953) 
on the reaction *N(p, «y), which also is resonant at this energy, makes 
it almost certain that the spin of the 1210 kev level is (4, +). Calculation 
shows that if this is so, interference terms from this level cannot be 


greater than 10% in the energy range considered; thus one cannot 
account for our results in this way. 
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5.2. The Fitting of the Angular Distributions 


If we accept the assignment (1, —) to the upper level, the level at 
340 kev must be (0, +) formed by p-protons. With these assumptions, 
the variation with energy of the coefficient of cos? @ was fitted by least 
squares and gave 


2¢ cos y=(0-80-0-06), 22=(0-06,+0-03). 


The fit is shown in fig. 2 (c). The negative cos? 6 term is at some energies 
certainly more than 0-25, and because of this it follows from eqn. (2) 
that cos y must be more than 0-34. 

To explain the observed cos@ term, and the observed values of 
1(150°)/Z(90°), a least-squares fit was made to the latter quantity for 
energies between 500 kev and 820 kev, in which range 5 (the phase 
difference between the two resonance denominators) does not change 
by more than 5°. The above quoted value of 2c cos y was taken. The 
quantities determined from the least-squares fit are 2? and F, where 


F={cos (6—8)—2c(p,/po)? cos (6B—S+y)}. . . . (3) 

The fit is shown in fig. 2 (a). Extrapolation to higher and lower 
energies involves a strong variation of 6 and in consequence a variation 
of F. A knowledge of y itself then becomes necessary ; because of our 
condition 0:34<cos y<1 the extrapolation can be made without great 
uncertainty. 

Some disagreement between theory and experiment, of the order 
of 10°% at the higher energies and 1°% at the lower energies, might be 
expected because of the effects of the level at 1210 kev. If we use only 
our knowledge that cos y>0-34, then 6=(20°+ 20°) or (116°+20°). Ifone 
makes the extreme assumption that cos y=1, one finds 6=(15°+-5°) or 
(121°-L5°) ; because of the possible effects caused by the 1210 kev level, 
one cannot decide between these alternative extrapolations, which are 
both shown in fig. 2 (a). 

5.3. Hxcitation Curves 


The extrapolated interference patterns are such that a level at 1050 kev 
of 150 kev width would appear, if established by observation at 138° 
to the incident proton beam, to be at about 1025 kev if 6=20°, and at. 
1030 kev if 6=116°. Although the magnitude of this shift is uncertain, 
its existence would suggest that the 1020 kev resonance observed at 138° 
in 15N(p, «) and the 1050 kev resonance in *N(p, y) (see Schardt et al., 
loc. cit.) represent the same level. 

The value of x? (the relative inherent excitation of the two levels) 
obtained from the J(150°)/Z(90°) curve is (0-033--0-006); this is in 
reasonable agreement with the figure (0-06,-+0-03) required to account 
for the observed cos?@ term. It is worth noting that the value of 2 
determines the excitation curve of the reaction at 90°, provided the 
penetrabilities are known. (At 90° interference between levels, of the 
type discussed by us, vanishes completely as may be seen by putting 
cos 0=0 in eqn. (2).) It is thus possible to derive yet another value of x? 
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from the observed excitation function ; using the published results of 
Schardt et al. (loc. cit.) we find x2=0-12, which seems significantly higher 
than our values. With this value the relative heights of the resonant 
peaks are well fitted, but the ‘ saddle’ in the 90° excitation curve in the 
neighbourhood of 450 kev is higher than one would expect. This was 
ascribed by Schardt et al. (loc. cit.) and by Fowler and Thomas (1953) 
to constructive interference of the type treated by Wigner (1946). 
It seems quite possible, however, that this discrepancy, as well as the 
anomalous value of x?, simply represents a shortcoming of conventional 
penetrability calculations when the conditions at widely different energies 
are being compared. 


§6. FurTHER Discussion OF RESULTS 
6.1. Application of Modern Dispersion Theory 


The phases ¢ and y may be predicted from the dispersion theory of 
Wigner and Eisenbud (see Thomson et al. (loc. cit.) and Seed and French 
1952). The predicted value of y is 84° with little dependence on energy 
and nuclear radius; this is in disagreement with the observed condition 
cos y>0-34, y<70°. The predicted value of ¢ is 25° for an interaction 
radius a of 5-:05x10-% em, and 16° for a=4:22x10-% cm, with little 
dependence on energy. One of our values for ¢ does not disagree with 
this, though it is evident that too much reliance should not be placed on 
the theoretical predictions. 

Fowler and Thomas (1953) have suggested that the level at 340 kev 
is (1, —) formed by s-protons, because of the large width for proton 
emission. From the values of Schardt et al. for the width and the cross 
section at resonance we infer J’,=(91--0-5) kev or (3:3-L0-5) kev. We 
may reject the first alternative as being rather high even for s-protons. 
If we accept the value 3-3 kev and correct it for the effects of level shift 
(Lane 1951, and Thomas 1953) we find 


I’, (true)=12 kev for p-protons, and 3-3 kev for s-protons. 
We now use the criterion given by Teichmann and Wigner (1952) : 


SY 


ra 


3h 
I’, (true)< ( =a) . 2kp, 


where .=reduced mass of system, 


a=interaction radius, 

k=wave-number of incoming proton, 

p,=penetrability coefficient. 
This requires I’, (true) <50 kev for s-protons, and <9-5 kev for p-protons, 
assuming in our calculations that a=5-05x10- em. Thus the Wigner 
criterion is just violated if we assume p-protons are responsible for 
forming the level at 340 kev. 

The calculation, in particular the correction of I" > for level-shift, is 

peculiarly sensitive to the assumed penetrability coefficient. In our 
calculations, we have assumed that the Coulomb potential is valid right 
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down to the nuclear boundary, a condition that is unlikely to hold. 
It may well be that I’, (true) is in fact less than Wigner’s upper limit, 
but that our knowledge of penetrabilities is not exact enough for us to 
use this criterion in marginal cases. On this view it would seem that there 
1s NO strong argument against our supposition that the 340 kev resonance 
is (0, +-) formed by s-protons. 


6.2. Isotopic Spin 

In the above discussion we have tentatively identified the alpha- 
emitting level formed at 1020 kev with the gamma-emitting level formed 
at about 1050 kev. From the yield of the radiative capture process 
?N(p, y) (Schardt et al., loc. cit.) we may infer that I, is (150+22) ev 
at 1050 key. Now the maximum expected width for gamma-radiation 
of this energy is 2000 ev for electric dipole, and 45 ev for magnetic dipole, 
using the formulae of Blatt and Weisskopf (1952). For any higher order 
of radiation it is of course much less. It is thus almost certain that the 
gamma-emitting level is (1, —), and that therefore it is indeed identical 
with the alpha-emitting state. One cannot, however, quite exclude the 
possibility that it is (1, +) and proceeds to the ground state by a magnetic 
dipole transition. 

If isotopic spin 7’ is a good quantum number, the alpha-emitting state 
must have 7'=0 (Adair 1952). Wilkinson (1953) has pointed out that if 
this is so, and if one identifies the alpha- and the gamma-emitting levels, 
as we have done, there is a clear violation of the isotopic spin selection 
rule AT 0-0 forbidden for electric dipole transitions (Radicati 1952 and 
1953). If the selection rule is to hold, then we must postulate two distinct 
levels, separated by less than 30 kev, and of the same width and spin, 
though possibly of opposite parity. This seems unlikely, and it is 
probable that isotopic spin is not a good quantum number for this level, 
which perhaps is not surprising at such a high excitation of the compound 
nucleus. 

§ 7, CONCLUSION 


The results provide yet another example of ‘interference’ between 
nuclear levels. It seems reasonably well established that the level 
formed at 340kev is (0,-+) and that the resonance at 1020 kev 
corresponds to a (1, —) level near 1050 kev. We seem to have a clear 
exception to the isotopic spin selection rules. Of the two observed 
phase-differences between coherent proton-waves in this reaction, one 
possibly agrees and the other probably disagrees with the predictions 
of the theory of Wigner and Hisenbud. 
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ABSTRACT 


It is pointed out that the 8-06 mev (1—) state of 4N is of 7—1, 
corresponding with the 6-10 Mev state of 4C. The electric dipole transi- 
tion to the 2:31 (0-++) 7=1 state of 14N is therefore forbidden by the 
isotopic spin selection rules. It is shown that the strength of this for- 
bidden electric dipole transition is less than 0-7°/ of that of the allowed 
electric dipole transition to the (1+) 7’=0 ground state of 4N; this 
implies that the contamination of the 8-06 Mev state by a state of 7’=0 is 
less than 2°% in intensity. 


§ 1. INTRODUCTION 


THIS paper continues the series (Wilkinson and Jones 1953—I ; Wilkinson 
1953 a—II) in which the stringency and general applicability of the 
isotopic spin selection rules is being investigated. The present work 
—as I—concerns the rule (Radicati 1952, Gell-Mann and Telegdi 1953) 
that electric dipole transitions in self-conjugate nuclei must be accompanied 
by a change of unity in the isotopic spin, i.e., transitions with 47’'=0 are 
there forbidden. This rule is seen to operate in several instances (in 1°O 
(one example)—see I; in 1°B (four examples)—Jones and Wilkinson 
1953 and forthcoming paper; in “N (one example)—present paper ; 
in 2°Ne (one example)—forthcoming paper) in suppressing such transitions. 
without change of isotopic spin. It is unlikely that this suppression is due 
to chance fluctuations in the electric dipole matrix elements since such 
fluctuations appear to be rather small as a rule in light nuclei (Wilkinson 
1953 b). That the inhibition may in fact be ascribed to the isotopic spin 
rule is made probable by the observation (in IT) that electric dipole transi- 
tions without change of isotopic spin are not discouraged in nuclei of 


qi | 


§ 2. StarEs oF 14N 


When protons of 550 kev bombard 1%C a state of width 32-5 kev is 
formed in 14N at 8-06 Mev (see Ajzenberg and Lauritsen 1952 for most 
references to level characteristics). This width is about 12% of the Wigner 
limit (Teichmann and Wigner 1952) if we assume s-wave proton capture 
and approaches or exceeds the single-particle reduced width if we assume 


p-wave proton capture. Since s-wave interactions in light nuclei show 


oe BE Ol a Se ee ee ee eee 
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large reduced widths while p-wave interactions show relatively small 
widths (A. M. Lane—private communication) this observation makes it 
very probable that we are dealing with s-wave protons and that the 
8-06 mev state is therefore (0—) or (I—). This conclusion is strengthened 
by the observation that the emission of gamma radiation by this state is 
isotropic with respect to the proton beam (Devons and Hine 1949). This 
level de-excites by radiation to the (1+) 7’=0 ground state of “N with a 
value of (2J+1)I,=3lev. This corresponds to a value of 
(2J-+-1) | M |2=0-09 where M is the matrix element relative to that of the 
single particle model (Weisskopf 1951); this value is appropriate to an 
allowed electric dipole transition (Wilkinson 1953 b) and so demands that, 
for the 8-06 Mev state, 77=1. Now the first 7=1 state of 14N is (0+) and 
is at 2-31 Mev, so we might expect to find the companion to the 8-06 Mev 
state of 14N in 4C at about 5-7 Mev. There is a state well established in 
40 at 6-10 Mev: Benenson (1953) has shown that it is formed by s-wave 
neutron transfer in the reaction °C (dp) “C and so must be (0—) or (1—) ; 
but Thomas and Lauritsen (1952) have shown that this state emits 
gamma-rays to the (0+) ground state of 4C and so it must be (l1—)—an 
observation supported by the fact that these workers also find that the 
internal pair-formation suggests an electric dipole transition. Thomas 
and Lauritsen also remark that the difference between 6-10 Mev, the 
excitation of the state in 14C, and 8-06—2-31 Mev, the separation of the 
7'=1 states in 4N, may be quantitatively correlated with the large reduced 
width of the 8-06 mev state of *N—the Wigner—Thomas shift (see Thomas 
1952). It is therefore highly probable that the 8-06 Mev state of 14N is 
(l—) 7=1. This is made certain by the work of E. Milne (private 
communications from R. G. Thomas and H. H. Woodbury) on °C (pp’) #C 
which shows the state to be of J=1. 

We must now note that the transition from the 8-06 Mey (l1—) 7=1 
state to the 2-31 Mev (0+) 7’'=1 state is an electric dipole transition that, 
on the basis of energetics alone, we should expect to be about 0-4 times as 
strong as the electric dipole transition to the ground state. But it is 
forbidden by the isotopic spin selection rule since 47'=0.* The history of 
this possible 5-75 Mev transition is a curious one, during which its reported 
intensity has steadily dwindled. Lauritsen and Fowler (1940) and Fowler 
and Lauritsen (1949) reported the intensity to equal that of the ground 
state transition ; Carver and Wilkinson (1951) (see Barnes, Carver, Stafford 
and Wilkinson (1952) for details) showed that this was not so and reported 
radiation of 5 81+0-25 Mev in about 7°% of the intensity of the ground 
state transition but entered the caveat that this might be contamination 
from the 6-14 Mev gamma-rays arising in the reaction 1°F (pa) 160 ; 
Woodbury and Fowler (1952) and Woodbury, Day and Tollestrup (1953) 
showed that 7% was indeed too generous an allowance and placed an 


upper limit of 2% on the relative abundance ; Hicks, Husain, Sanders and 
anne has also been noticed by Gell-Mann and Telegdi (1953) and Benenson 
Jodo). 
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Beghian (1953) also failed to find the transition (corresponding to a limit 
of about 4%). i 

We may now understand the absence of this transition as an isotopic 
spin effect as remarked above, and it becomes of considerable interest to 
try and find it or place a closer limit on its occurrence in order to be able to 
make as strong a statement as possible about the contamination of the 
T'=1 state at 8-06 mev with a 7=0 component such as would enable the 
transition to take place. (We may suppose the 7’=1 state at 2-31 Mev 
to be relatively pure since it is certainly removed by several Mev from any 
suitable (0-+-) 7’=0 contaminant.) 

It is appropriate to remark here that the discovery of such transitions 
of determinable radiative width that violate the isotopic spin selection 
rule is the best way of actually determining isotopic spin impurity as 
opposed to setting a limit on this quantity. This is because the unin- 
hibited electric dipole transitions seem as a rule to show a constancy of 
(2J-+1) | M |? to within a factor of about 2-5 (Wilkinson 1953 b) and so the 
measurement of an inhibited radiative width would determine the isotopic 
spin impurity with comparable precision. The other form of violation of 
the isotopic spin rules, the formation of, say, a T7=1 state from 7'=0 
particles such as is observed in ®Li («y) 1°B (Jones and Wilkinson 1953 and 
forthcoming paper) is less helpful since uninhibited reduced widths are 
subject to much greater fluctuations (of the order of a factor of 100) 
and the correction for the penetration of the Coulomb barrier is also often 
rather uncertain. This second form of violation can only properly be 
interpreted as placing a lower limit on the isotopic spin impurity (through 
the use of the Wigner limit for the reduced width). 

It was for this reason that we undertook a more extensive search for the 
forbidden transition than those referred to above. 


§ 3. THE Reaction C (py) 4N 


We bombarded a thin target of electromagnetically-separated '°C with 
protons of 560 kev and observed the gamma-rays at 90° to the proton beam* 
with the aid of a 1” cube of Nal(TI) and an EMI multiplier type 6262. 
Eight runs were made, separated by standardizing irradiations using 
sources of 137Cs (662 kev), ®°Co (1-172 and 1-332 Mev), ?C (4:43 Mev) and 
16Q (6-14 Mev). 

The combined irradiations are shown in the figure. 

The standard deviations may be estimated from the number of counts per 
channel and are everywhere smaller than the radius of the circles. Peak 
A is due to the 8-06 Mev gamma-rays produced in the direct transition to 
the ground state of “N; it is due to pair creation with subsequent 
escape from the crystal of both annihilation quanta (the two-quantum 
escape peak). 
ee te ee ee eee ee en ee ee 

* The angle of observation is of no consequence since, as we have seen, the 


reaction is induced by s-wave protons. 
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It had been reported by Woodbury and Fowler (1952) and Woodbury 
et al. (1953) that, in addition to this principal ground state transition, there 
is a weak transition to the state at 3-948-L0-015 Mev* giving a gamma-ray 
of energy about 4-1 Mev. This transition is well seen in our spectrum as 
peak B, due, as peak A, to pair creation with escape of both annihilation 
quanta ; at C we have the step due to capture of one annihilation quantum 
(the one-quantum escape peak), and at D the final step due to capture of 
both quanta (the no-quantum escape peak). At E we have the two- 
quantum escape peak due to the unhappily-ubiquitous 6-14 Mev gamma-rays. 
from fluorine contamination (see the calibration distribution Y below it); 
at F we see signs of the one-quantum escape peak of the same radiation. 


o 


re) 


@ 


BR 


TENS OF THOUSANDS OF COUNTS PER CHANNEL 


(e) 


oO 8 16 


24 32 40 
CHANNEL 


Pulse distribution due to the gamma-rays from ™N detected in a Nal(TI) 
crystal. For the meaning of the arrows see the text. Distribution X 
is part of that due to the 4:43 Mev gamma-rays from !C following 
N(pa)?C; distribution Y is part of that due to the 6:14 Mev gamma. 
rays from 1%O following 1°F\pa)!®O. (Distributions X and Y do not 
have the same ordinate scale or zero as the main distribution.) 


At G we have placed an arrow in the position expected of the two-quantum 
escape peak from the forbidden 5-75 Mey transition that we seek. There is. 
certainly a significant rise at about this position and we may be tempted 
to identify it with this forbidden transition. If we do this then the 
intensity of the supposed transition is about 0-7% of that to the ground 
state. But there is an alternative explanation of the peak G. ™N 
possesses a state at 5-08+-0-02 mMevt which combines directly with the 
* Thi: 4 . Car * =e 5 i! = 
Bone saat as ae 7 eat ements by Thomas and Lauritsen (1952) 


+ This energy derives from measur i 
; surements of Thomas and Laurit 
*C(dn)'*N and Woodbury et al. on !C\py)!N. eS 
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(1+) ground state and also shows a competing branch of intensity 26° to 
the (0-+-) 7 =1 state at 2-31 Mev (Woodbury et al. 1953) ; no corresponding 
state is known in MC so it is probably 7=0. This makes it probable that 
it is (1-+-) since J=0 is excluded by the transition to the 2-31 Mev (O-+-) 
T’=1 state which also makes J>2 improbable, while (1—) would give 
two electric dipole transitions of which that to the ground state would be 
isotopic spin forbidden. This parity assignment agrees with Benengon’s 
(1953) tentative conclusion from °C (dn) !4N stripping while Woodbury 
et al. suggest J=1 or 2. If our suggestion of (1+) is correct we must 
expect a number of transitions to take place to this level from our (1—) 
level at 8-06 Mev—by the energy alone about 5°% of the number to the 
ground state. Indeed Barnes et al. (1952) reported a possible gamma-ray 
of energy about 5-2 Mev in about 1-5°/ abundance while Hicks et al. (1953) 
reported a 5-1+-0-2 Mev line in 4% abundance. These figures, corrected 
for the branching ratio observed by Woodbury et al. are well consistent 
with the expected abundance of about 5°, for the transition to the 5-08 Mev 
state. Now if this 5-08 mev line is really there its one-quantum escape 
peak will fall at 4:57 Mev, rather close to the principal two-quantum 
escape peak at 4-73 Mev from our forbidden transition. The position of 
this one-quantum escape peak is indicated by the arrow at G’ which fits 
the observed peak about as well as that at G corresponding to the forbidden 
transition. The principal two-quantum escape peak of the possible 
5-08 Mev line would fall at 4:06 Mev at the position indicated by the 
arrow D’ and so be largely hidden by the no-quantum escape peak from 
the 4:1 Mev line. If we interpret the peak GG’ as a cne-quantum escape 
peak it implies a relative intensity for the 5-08 Mev line of 2°% which agrees 
well with the earlier estimates. A slight argument against accepting 
the peak GG’ as due wholly to a 5-08 mev line is that the form of the 
successive peaks B, C and D is in excellent agreement with what is 
expected of a 4:1 Mev line (the nearby 4-43 mev line of C provides a good 
measure of this—see distribution X of the figure); if we subtract from 
the BCD group the two quantum escape peak at D’ appropriate to the 
one-quantum escape peak G’ the remaining BCD group becomes a little 
mis-shapen. A 2% abundance of a 5-08 Mev line implies, using the 
branching ratio of Woodbury et al. (1953), a 2-7% abundance of the 
preceding 3-0 Mev line. We have sought this line but not found it. The 
search is made difficult by the fact that the photo-peak of this line is 
hidden by the two-quantum escape peak of the 4-11 Mev line, while the 
principal two-quantum escape peak almost coincides with the Compton 
edge due to the 2-31 mev line. A 3% abundance would have escaped 
detection. Hicks et al. have reported a 3-05-+-0-1 Mev line in 15% 
abundance; this appears to be excluded by our results and by those of 
Woodbury et al. 

We must also bear in mind the possibility that this bump GG’ is indeed 
the two-quantum escape peak of a gamma-ray of about 5-7 Mev, but one 
that follows a transition from the initial 8-06 Mev level to one of the 
known states at about this excitation. 
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§ 4, THe Isoropic Sprn Impurity 

The present results, therefore, are not conclusive. We have seen that it 
is very reasonable that the 5-08 Mev line should exist and that there is some 
evidence for it so we must leave our 0:7% as a limit rather than as a deter- 
mination of an abundance. This corresponds to a contamination of the 
8-06 Mev 7’=1 state with less than 2° of 7’=0 in intensity ; in the nota- 
tion of Radicati (1953) «,2(0)<0-02. 

If the evidence in favour of the 5-75 Mev line were acceptable we would 
say «,2(0)~0-02. This would be a rather large contamination but perhaps 
not unreasonable in view of the fairly high density of states at this 
excitation in 4#N. A higher contamination has been observed at 13-09 Mev 
in 160 where «,2(0)>0-03 (Wilkinson 1953.c). (As always, when we are 
concerned with a self-conjugate nucleus, we must bear in mind that these 
results in themselves are explicable in terms of charge-symmetry and the 
charge-parity rules alone (Kroll and Foldy 1952) though the identification 
that we have made in § 2 between the 8-06 Mev state of ##N and the 
6-10 mev state of !4C tells of course, for charge-independence.) 


§ 5. OTHER GAMMA-RAyYS 

Woodbury et al. (1953) report the energy of the cascade gamma-ray as 
4-0+0-1 Mev ; Hicks et al. (1953) give 4:-45-+-0-1 Mev and hypothesize a 
new level in 14N in explanation. This is a serious discrepancy. We have 
measured this energy as 4:11-+-0-02 Mev from the data presented in the 
figure and also from a further series of 14 runs interleaved with calibrations 
from the nearby 4:43 Mev gamma-ray from !#C made by 1°N (px) #?C; 
these runs were made with high dispersion in order to examine the 4-1 Mev 
peak alone in more detail. This energy agrees perfectly with what is 
expected of a transition between the 8-06 and 3-95 Mev states. Woodbury 
et al. (1953) have reported the branching from the 3-85 Mev state to that at 
2:31 Mev to be twice as probable as the direct transition to the ground 
state. This would imply that the 3-95 Mev gamma-ray should be present 
in one third the intensity of the 4-11 Mev line and it would therefore be 
rather difficult to separate from that line. We have seen no sign of the 
3-95 Mev line and believe that its intensity may be slightly less than 
Woodbury et al. claim, though we do not feel strongly on this point. 

We have, incidentally to our accurate measurement of the 4:11 Mev 
line, observed the last element of the cascades and determined its energy to 
be 2-307-40:015 Mev. Thomas and Lauritsen (1952) from a study of 
°C (dn) *N find 2-310-+-0-012 Mev as the energy of this transition ; 
Woodbury et al. find 2-32-++0-02, 2-35-+.0-04, 2-32-++0-04 and 2:32-L0-02 mev 
yielding an average of 2-323-+-0-015 Mev. These estimates are consistent 
with our own and may be averaged to give an energy for the state of 
2-313-+0-007 Mev. 
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SUMMARY 


This paper is concerned with a rigorous treatment of the dielectric 
polarization of a model consisting of unpolarizable point dipoles fixed 
on the sites of a lattice but free to rotate. The dielectric constant of 
this system is expanded in inverse powers of the temperature as far 
as 7-3, It is shown that for small dipole moments the Debye and 
Onsager approximate formulae both overestimate the polarization, the 
Onsager result being the more accurate of the two. Previous treatments 
of this problem did not make adequate allowance for the long-range 
effects of boundary field polarization. 


$1. INTRODUCTION 


One of the simplest statistical models that can be used as a basis for a 
theory of dipolar interaction in condensed polar media is the dipolar 
lattice. In this model each molecule is represented by a point dipole 
fixed at a lattice site but free to rotate in the electric field of all the others. 
Although the partition function cannot be evaluated in closed form, 
considerable progress can be made with a high temperature expansion 
(Van Vleck 1937, Pople 1952) and the model has been used for a discussion 
of the cohesion of polar liquids. 

The present paper is devoted to the calculation of the polarization of 
the dipolar lattice model in a uniform electric field, the aim being to 
expand the dielectric constant in inverse powers of the temperature by 
a method analogous to that used in the free energy treatment. We 
shall only consider the problem of a lattice of permanent dipoles and 
shall not be concerned with modifications (such as inclusion of polariza- 
bility or quadrupole interactions) which may make the model more 
realistic. The results for the simple model are of considerable interest, 
however, for, being rigorously derived, they can be used as a significant 
test of the approximate formulae of Debye (1929) and Onsager (1936) 


This model for dielectric polarization has been investigated previously | 


by Van Vleck (1937) and Rosenberg and Lax (1953). These authors, 
however, do not make adequate allowance for the effect of boundary 
ee 
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field polarization. Kirkwood (1939), in an important paper on this 
topic, showed that the dielectric constant « of a system consisting of N 
permanent dipoles in a volume V is given by 


4nN e+2y.. M 

er ae 
where M is the total moment induced in a spherical specimen by any one 
of its dipoles maintained in a fixed position and orientation p. Because 
of the reaction field due to the external boundary, M is not concentrated 
near uw. Kirkwood therefore introduces a local moment p which may 
be defined as the moment of a sphere surrounding p which is large com- 
pared with molecular dimensions but small compared with the whole 
specimen. In terms of this moment, « is given by 


e—l= 


(1.1) 


47N 3e we. we 


Now the procedure of Van Vleck and Rosenberg—Lax is equivalent to 
a direct attempt to obtain a high temperature expansion for y.M. 
As the coefficients are obtained as rapidly convergent lattice sums, it 
seems possible that they are really obtaining the local moment pw rather 
than the total M. The details of the analysis later in this paper indicate 
that this is so and the results have to be modified accordingly. 


§2. High TEMPERATURE EXPANSIONS FOR THE DIELECTRIC CONSTANT 


Before obtaining the rigorous expansion for the dipolar lattice model, 
we shall give the corresponding expressions for the approximate formulae 
of Debye and Onsager. The Debye formula, corresponding to putting 
M=p. in (1.1) can be written 

Bese re accuse ee ot (21) 
where x is the non-dimensional quantity 47Nyu?/3Vk7'. We shall find it 
convenient to expand the quantity 3(e—1)/(e+2) rather than «—1 in all 
formulae. 3(«—1)/(e+2) is 47 times the resultant moment of a spherical 
specimen per unit volume per unit field before insertion of the dielectric 


material. K 
The Onsager formula, obtained by putting up in (1.2) can be written 


3(e—1)/(e+2)= A+ ge) (1+ ge+o%)®@—1 440}, 2. (2.2) 
or, expanding for small x 
B(e lie 2a ser Oa). a. » (2:3) 


The two formulae are seen to agree in giving no term in 2? but differ in 


the coefficient of x. 
The treatment of the dipolar lattice model will be based on the exact 


Kirkwood equation (1.2). If we write x’ for 47Ny.p/3VkT we have 
B(e— 1) /(e- 2) =a — 2024 Oe), 1... (2.4) 
so that it only remains to find x’ in terms of z. 
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Suppose that we have N dipoles, fixed in position at points R,...R, 
of a lattice, but free to rotate. The dipolar interaction energy is then 


given by 


W=2 Wj; ye (2p 

i<j 

where ‘ 
10 =e pgeltgp: ah ek ol oe rr 


Here Greek suffixes are used as tensor suffixes and ¢{} is the dipole 
interaction tensor 


D)— Ro{R,70,,—3R,,Rija}, 
p j { ij “ap y] ues (2.7) 
R,;=R;—R,. 

Now suppose that v; is a region surrounding the site R; which is of the 
required intermediate dimensions. Then we may write 
u. p—pe= a* U;- U5, ° . e . ° e 5 (2.8) 
HAD 


where >* indicates summation over lattice sites in v; and the bar represents 
an average over the canonical ensemble for the whole specimen (in the 
absence of an external field). Now 


N 
Be. w= {us uw; exp \- E w/k} dw / Jexp {= E wylkt} dw, (2.9) 
k<l k<l 


Jdw having been written for integration over the angular co-ordinates 
of all dipoles. The numerator and denominator of (2.9) can now be 
expanded in powers of x. The denominator is (47)*{1+O(2?)} and the 
leading terms occurring in the numerator are (for 1+7) 


B+ pj dw=0, 
ic + By Wis deo= | within mighty do= (40) =0, - + (2.10) 


fe - PjWp, dw=0. 


It follows that (uw. p/u?)—1 is at least as small as O(a2). The only non- 
vanishing integrals in the next term of the expansion have integrands 
Pej» Wj Wn (kKAt, 7) and we find that 
(we. w/e?) —1=(4ar)-* (p?k2T2)-1 De® D ju + PejW jz W jy, dw. . (2.11) 
(Ai) k 
It should be noted that the sum over k in (2.11) is over all the lattice 


sites in V (except 7 and j) whereas the sum over j is only over those in 
v,; (except 7). Integration over angle gives 


(y . pu?)—l= dp (pt/PT2) S* Te, | |... (2,12) 
At) k 
en €op = 8Ry Rg °{3(Ryz « Ry,)?—By2Rj2}- . . . (2.13) 


The evaluation of the double summation in (2.12) requires some care. 
For any given value of j in v, the k-summation is over a large region 
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round both i and j. Thus k-summation (again for fixed j) converges if 
extended to infinity (by inspection of (2.13)) and the sum will be a 
function of R,,. Since the region », itself is macroscopic, we can then 
carry out a second infinite summation over j7. As a result we may write 


(2.12) as a doubly infinite summation, the k-summation being carried 
out first. 


(ue. p/?)—l=35 (4/k2T) 2 ares? € JF), pee, (2.84) 
HAD E 
(Infinite summation in this context means summation over sites extending 
infinitely from the neighbourhood of i.) 

It should be emphazised that this clear separation as a consecutive 
doubly infinite summation is only possible because of the introduction of 
the two regions v;and V. If both summations were over V this argument 
would break down, for there would no longer always be a large region 
of k-summation outside both i andj. This is the point that has not been 
sufficiently considered by Van Vleck (1937) and Rosenberg and Lax (1953). 
They assume, in effect, that the double series on the right-hand side of 
(2.14) gives (uw .M/2)—1 (eqn. (1.1)) and includes the contribution of the 
boundary field polarization. 

It is convenient to add the term j=i to (2.14) and subtract again 


(ye. /e2)—1 =F (uke) ) {Beeps Gk) — Sele me ee (2515) 


We next show that the double summation in (nes is zero. This is on 
account of the following two results : 
le Jexp ep. is convergent for all j(€{}<% being asymptotically propor- 
tional to R;,° for k far from both i and j). 
2 Berg ab ap Aap — 0. (oll We ee geek eee (2a) O) 
Jj 


‘This follows directly from (2.7) and the symmetry of ae lattice. 
From the second of these results it is clear that EEE @) converges 


uniformly for all p (that is for all finite k- Rteaconen We therefore 

have sufficient conditions to interchange the order of summation (Jeffreys 

and Jeffreys 1950, § 1.125) and the double summation vanishes. 
Evaluating the remaining term in (2.15) using (2.7) we now have 


Bet =1— FMT) ERG +. 


—1—b22+0(x'), ee 2 (2:17) 
where b=2(V/4rN)? J Ree, ee woe es (2 8) 
kK #4) 


is a number depending on the type of lattice. For the simple, body- 
centred and face-centred cubic lattices 6 has the values 0-10642, 0-09196 
and 0.09154. 
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§3. Discussion 
The high temperature expansions for dielectric polarization obtained 


in this paper are 
Debye formula 3(e—1)/(e+-2) =x 


Onsager formula =x—0-22222¢3+.... 
Simple cubic lattice =2—0:3286423+.... 
Body-centred cubic lattice =a«—0-314182°-+.... 
Face-centred cubic lattice =—7—0-31376z'-+. . 


For reasons discussed in the previous section, these aifice from the 
figures given by Van Vleck (1937) and by Rosenberg and Lax (1953). 

It should be noted that the expansions for the three lattices differ 
only slightly. This suggests that the lattice formulae should, in fact, 
give a good approximation to the polarization for any isotropic arrange- 
ment of dipoles as in a liquid. When the accurate results are compared 
with the corresponding expansions of the approximate formulae, it is 
clear that, for small dipole moments at least, both the Debye and Onsager 
equations overestimate the polarization, the latter being the more 
accurate of the two. The dielectric constants of actual polar liquids 
will also be influenced by other factors such as polarizability and other 
short range intermolecular forces, but the results of this paper suggest 
that the Onsager equation may be substantially in error even for molecules 
interacting as permanent point dipoles. 


The author is indebted to Dr. 8. F. Boys for discussion on this subject. 
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ABSTRACT 


A method is described for the realization of very high pressures at 
liquid helium temperatures. First results are given of its application to 
studies of superconductivity under pressures up to 44.000 atm. Bismuth 
is shown to become superconducting under pressure. 7", for tin is shown 
to vary approximately linearly with volume up to 17 500 atm. The value 
of dT./dp for thallium is found to be negative, in disagreement with 
previous work. 


§ 1. INTRODUCTION 


ReEcENT theories of the origin of superconductivity make it desirable to 
obtain more experimental information on the effect of varying the pressure 
on the properties of superconductors, to supplement the growing mass of 
data on the effect of varying the isotope in a given superconducting metal. 
Indeed, the pressure (or volume) and the mass number would seem to be 
the only variables whose effect might be capable of immediate theoretical 
interpretation. In particular, it is of interest to determine the influence 
of pressure on the transition temperature, 7',, and on the value of the 
quantity dH,/dT at T, (where H, is the threshold field at temperature 7’), 
and also to find whether the application of high pressures can destroy 
superconductivity, or cause it to appear in metals not normally super- 
conducting. 

The present paper describes some experiments of this kind which we 
have carried out at pressures up to about 40 000 atm.—over twenty times 
higher than the maximum pressures previously employed in similar 
investigations. The increase in the working range of pressure has been 
made possible by the use of a technique due to Bridgman adapted in a 
special way for use in low-temperature enclosures. 

In a further paper the implications of some of the results obtained will 
be examined from the general standpoint adopted in the theories of 
Frohlich (1950 a, b, 1951) and Bardeen (1950 a, b, 1951 a, b, c), in which 
superconductivity is attributed to interactions between electrons va their 
interactions with the lattice. 
a UEEEEE EEE EEEEEEEEE REESE 


* Communicated by the Authors. 
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§2. Summary OF Previous WoRK 


The earliest systematic study of the effect of compression on the 
properties of superconductors was carried out by Sizoo and Onnes (1925), 
who applied pressure to wires of tin and indium by compressing helium 
into the vessel in which they were mounted. For both metals the 
transition temperature (7',), at which the resistance fell to half its residual 
value, was found to be depressed. However, the maximum pressures used 
were no greater than 300 atm., and the maximum observed depression of 
T, was about 0-005°. The later discovery by Keesom (1926) of the 
solidification of helium under pressure threw doubt on the results obtained 
by Sizoo and Onnes for all but the lowest pressures used, and removed the 
last hope of using a transmitting fluid for the application of conventional 
high-pressure techniques to experiments at liquid helium temperatures. 

An advance in technique which enabled appreciably higher pressures 
to be reached at low temperatures was made by Lazarew and Kan (1944), 
who generated pressures of about 1750 atm. by allowing water to freeze 
in a thick-walled ‘bomb’. Electrical leads leading out of the ‘ bomb ’ 
made it possible to measure the resistance of wires mounted in it as they 
were cooled, under pressure, to liquid helium temperatures. Of the 
metals (thallium, indium, tin, mercury, tantalum and lead) subsequently 
investigated by Kan, Sudovstov and Lazarew (1948, 1949), all but thallium 
showed a depression of 7’, under pressure. Using the same technique, 
Alekseyevski and Brandt (1949, 1952) investigated the non-superconduct- 
ing metal bismuth and its superconducting compounds RhBi,, NiBig, 
KBi,, LiBiand Au,Bi. No signs of superconductivity could be observed in 
metallic bismuth at pressures up to 1050 atm. Compression lowered the 
transition temperatures of LiBi and Au,Bi but raised those of the other 
compounds named. The absolute magnitude of the quantity dT'./dp 
in these experiments ranged from 5-8 x 10~1!! deg. dyne~1! cm? for tin to — 
almost zero for RhBi,. The maximum observed change in 7’, was there- 
fore about 0-1°; in most cases the total change was much smaller. The 
fact that the transitions observed were fairly sharp was taken to imply that 
the pressure was uniform throughout the volume of the superconductor. 

A number of other experiments, on specimens under non-uniform 
distributions of stress (such as linear extension), have also been carried 
out by the authors named above, but we do not discuss the results here. 


§ 3. EXPERIMENTAL METHOD 

In order to reach much higher pressures it was decided to adapt for 
our purpose the technique due to Bridgman (1935 a, 1949, 1950, 1952) 
by which quasi-hydrostatic pressures up to about 400 000 atm. may be 
generated in thin solid specimens. This technique depends on the fact 
that a small area of a massive block of metal is able to withstand stresses 
far greater than the normal breaking stress of the bulk material, because 
of the support afforded by the surrounding metal. In applying this 
princip'e to the study of the shear strength or electrical resistance of 
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metals under very high pressures, Bridgman interposes a disc-shaped. 
specimen between the plane faces of two truncated conical bosses, each 
machined on massive blocks of steel or Carboloy. As the blocks are 
forced together, very intense stresses are developed over the small area 
of contact, and the specimen and compressing surfaces deform in such a 
way that the pressure in the specimen approaches a hydrostatic distri- 
bution locally, with a value which decreases radially from the centre to 
the perimeter. The actual form of the variation is difficult to estimate, 
but it was believed by Bridgman that a large central part of such a 
specimen would be under fairly uniform pressure. This question is 
discussed elsewhere (Chester 1953), and we shall only mention here certain 
indications arising from the present results which suggest that this is 
true for the metallic specimens so far investigated by us. (The values of 
mean pressure quoted in the present paper are equal to the total force 
applied in a given experiment divided by the area of the specimen under 
stress.) 

In order to avoid the necessity for cooling to the working temperature 
the rather bulky apparatus which would be required to transmit a 
sufficient force to the specimen, we have mounted our specimens in a small 
self-contained device, which we refer to here as a ‘clamp’. Pressure 
may be generated in the spec’men at room temperature by applying a 
force to the clamp by means of a hydraulic press, and after a suitable 
adjustment the clamp may be removed from the press and transferred 
to a small cryostat, with the state of stress in the specimen still maintained. 
Apart from the simplification in technique wich results, there is an impor- 
tant advantage in the fact that the force is applied to the specimen at 
room temperature, and not at the working temperature in the liquid 
helium region. With this procedure the specimen, particularly if of a 
soft metal, remains in a fairly well annealed condition throughout. 
If the force were applied at liquid helium temperatures, annealing would 
not occur and the specimens stressed in the way described would be, if 
not actually disrupted, severely ‘ cold-worked ’—a condition which is 
known (Hilsch 1951) to affect the superconducting properties profoundly. 

The design of the clamp is illustrated in fig. 1. The specimen, in the 
form of a thin disc, is set between the truncated compressing cones C. 
For the application of pressure, the clamp is supported by its rim R in the 
cylindrical jig J, a plunger G is inserted in the axial hole H in the screw W, 
and the whole unit placed between the rams of the press. When force is 
applied between the plunger and the base of the jig, the pistons P of the 
clamp—and therefore the specimen—are put under compression, whilst 
the outer wall of the clamp is put into tension. The various members 
can now be ‘ locked ’ in their state of stress by tightening the screw. 

In practice itis not possible to preserve the full stress because of friction in 
the screw-thread, and the following procedure is adopted to ensure that the 
fraction preserved is as large as possible and is accurately known : When 
a suitable value of the force is reached, the screw is tightened progressively 
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and the hydraulic press relaxed at such a rate as to keep constant 
the strain in the clamp—as indicated by a resistance strain-gauge cemented 
to its outer wall. As the press is relaxed further, the screw becomes more 
difficult to turn and after a certain stage no further tightening is attempted. 
At this point the press is relaxed completely and the fraction of the original 
strain which has been retained (in practice never less than 90%) is noted 
from the indication of the strain-gauge. The clamp is now transferred to 
the cryostat. No change in the pressure applied to the specimen is to be 
expected on subsequently cooling the clamp to the working temperature 
because its load-bearing members are made wholly of one material, and 
therefore contract equally. 


Schematic diagram of clamp (for explanation of symbols see text). 


The superconducting transition of the compressed specimen is observed 
by a magnetic method. This is to be preferred to a method depending 
on the measurement of resistance, because it avoids the difficulties of 
electrical insulation at high pressure, and is less likely to lead to the 
spurious effects associated with the existence of filaments of abnormally 
high critical field—as reported by Lazarew and Galkin (1944) for inhomo- 
geneously strained superconductors and by de Haas and Voogd (1930) for 
inhomogeneous superconducting alloys. The method best suited to the 
geometry of the clamp is to wind secondary coils 8, in series, round each 
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cone and to apply an alternating field along the axis of the clamp by means 
of a solenoid coaxial with it. The e.m.f. induced in the secondary coils 

_is amplified, rectified and displayed on a micro-ammeter. Because of the 
large demagnetizing factor of a transverse disc, reasonably sharp transitions 
are to be expected only if the alternating field is of small amplitude. We 
have found also that it is not possible to make use of the increased sensiti- 
vity to be expected of high-frequency measurements because of skin 
effects in the material of the clamp. In consequence, most of the work 
has been done with alternating fields of amplitude about 0-1 gauss and 
frequency 30 c/s. 

Preliminary investigations, made with clamps of various steels, indicated 
that the proximity of ferromagnetic materials complicated the transition 
of the superconductor and made accurate determination of 7’, difficult. 
All the quantitative work to be described has therefore been carried out 
with a clamp of beryllium-copper—a material which remains weakly 
paramagnetic down to liquid helium temperatures. For experiments 
up to mean pressures of 16 000 atm., the circular area in contact with the 
specimen was about 5 mm in diameter, and the thickness of the specimen 
after compression was usually less than 0:01 mm. Some rough work up 
to higher pressures (about 40000 atm.), has been carried out, using 
conical tips of smaller area, of ‘ Vibrac’ steel. 

Because of the existence of a radial pressure gradient in the specimen it 
was expected that for most superconductors under pressure the outer part 
of the disc would become superconducting at a higher temperature than 
would the centre part, and it was thought that the superconducting rings 
thus formed might tend to mask the transition of the centre. To avoid 
this complication a small wedge of a metal, known to be not superconduc- 
ting under the conditions of the experiments, is inserted into a suitable 
cut in the circumference of the disc before compression, to prevent the 
formation of complete superconducting rings in the periphery. Aluminium, 
normally a superconductor below 1-2°K, has been found suitable for use 
with the metals so far studied. (An independent test was carried out to 
ensure that aluminium did not become superconducting under pressure 
above 2°x—-the lower limit of temperature of the present experiments. ) 
This precaution was not, of course, necessary in the exploratory work on 
non-superconductors. 

A Helmholtz pair located outside the cryostat, with the axis of its coils 
in the plane of the specimen, permitted transitions to be observed in known 
magnetic fields. The demagnetizing factor of the disc in this direction 
appears to be sufficiently small to be ignored, because the form of transi- 
tion curve obtained—though not of course its location on the temperature 
axis—is not sensibly altered by varying the field (see fig. 4). 

Temperature control in the region 2° to 4-2°K was effected by varying 
the pressure over a bath of liquid helium in which the clamp was immersed, 
the temperature being deduced from the value of vapour pressure. Above 
4-2°K a charcoal desorption method was used and temperatures measured 
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by means of a carbon resistance thermometer of the type described by 
Clement and Quinnell (1952). 


§ 4, EXPERIMENTAL RESULTS AND DIscUSSION 

The results on the metals studied are now presented and discussed. 

separately. 
Tin 

Tin, according to the Russian workers, shows the largest change in f bee 
under pressure. It is the most extensively studied superconductor,. 
having its transition temperature in a range where very accurate control 
of temperature is possible, and it was therefore the obvious choice for a 


first detailed investigation. 

Superconducting transitions were observed by cooling through the 
transition temperature under mean pressures of 1, 10 500 and 16 000 atm. 
The results are shown in fig. 2, where the alternating flux in the secondary 
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Transitions of tin in zero field at differing pressures : (a) 1 atm., (b) 10 500 atm., 
and (c) 16 000 atm. mean pressure. Flux in arbitrary units. 


coils is plotted, in arbitrary units, against temperature. The width of the 
curve at 1 atm. provides a reasonable measure of the limit of resolution of 
the method for a given alternating field, and we then interpret the 
spreading of the curves at high pressure as due to the pressure gradient in 
the specimen. Inspection of the curves suggests that over about three- 
quarters of the area the pressure is uniform within our limits of resolution 
(the small drop in flux at the high-temperature end of the curves is pro- 
bably due to metal extruded beyond the area of the compressing cones). 
By assuming a linear gradient over the remaining area we can arrive at 
estimates of the maximum pressure in the specimen. For a mean pressure 
of 16 000 atm., for instance (curve c), we estimate that the pressure 
over the central three-quarters of the specimen lies between 16000 and 
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19 000 atm., which we write as 17 500-1500 atm. We believe the results 
may fairly be stated as follows : 

At latm., the value of 7’, determined by us agrees with the normal 
value (3-73°K), within experimental error. 

At 11 500+1000 atm., 7',.=3-21-+40-07°xK. 

At 17 500+1500 atm., 7.,=2-9+0-1°K. 

It is important to emphasize that the curves represent reversible 
behaviour and the effects shown are therefore not due to work-hardening 
(or ? cold-working ’). After relaxing the pressure the curve obtained with 
a given specimen is indistinguishable from that obtained before the 
application of pressure. 

Combining these with the results of Kan, Sudovstov and Lazarew (1948), 
and making use of Bridgman’s values (Bridgman 1949) for the compressi- 
bility of tin, it is possible to plot 7, against 4V/V, where AV is the change 
in volume due to the increase in pressure. This is done in fig. 3, where it 
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Relation between transition temperature and change in volume for tin under 
pressure. The areas enclosed by the dotted lines indicate the estimated 
experimental uncertainty in the present work. 


may be seen that the variation is linear, within experimental error, and that 
the earlier and the present results lie on the same line. 
The relation derived from the graph is 


J’, (°K)=3-73-+29 (AV/V) 
where the coefficient of 4V/V may be in error by +10% at the lower 


pressures and by +20% at the higher pressures. 
We have also observed transitions under the highest of these pressures 
for a number of differing values of magnetic field (in the plane of the disc). 
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The curves obtained are shown in fig. 4. A field of 50 gauss depresses the 
foot of the curve by about 0-37°. The corresponding depression at 
latm. (see for instance Shoenberg 1952) would be 0:35". Thus the 
quantity dH,/dT at the temperature 7’, remains substantially constant 
under pressure. At a given temperature H, is of course lowered by 
pressure. 
The magnitude of the quantity Ho, equal to (H,)7—o, is of theoretical 
significance since it is related to the difference in energy between the 
normal and superconducting states at absolute zero. Our experiments 
have not yet been carried to a sufficiently low temperature to determine 
the form of the curve of H, against 7’, but we may point out that if it has 
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Transitions of tin under 16 000 atm. mean pressure in differing magnetic fields : 
(a) zero field, (b) 23-6 gauss, (c) 47-2 gauss. Flux in arbitrary units. 


the approximately parabolic form characteristic of tin at 1 atm., then 
the present results imply that H,j—and therefore the energy difference 
between the normal and superconducting states—is lowered by 
compression. 
Lead 
We have confirmed that 7’, for lead is lowered by compression, but have 
not yet obtained detailed results on this metal. 


Thallium 

Because of the anomalous behaviour reported for thallium by the 
previous workers (Kan, Sudovstov and Lazarew 1949), using resistance 
measurements, we have investigated it by our method up to higher 
pressures. 

In two runs under mean pressures of 11700 atm. and 13 400 atm. 
respectively, no signs of superconductivity could be observed above 
2-35°K although an uncompressed specimen began to show superconduc- 
tivity at 2:39°k. The whole transition curve was displaced by pressure 
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to a lower temperature, 7’, being lowered by about 0-06° under 13 400 atm. 
The change is much smaller than that observed for tin, but is in the same 
direction. This result disagrees with the observation of Kan, Sudovstov 
and Lazarew ; we may mention that the behaviour of thallium in the 
experiments of these workers was stated to be very sensitive to the 
method of mounting the specimen. 


Bismuth 


Bismuth, though not normally a superconductor, is obviously very close 
to being one. Hilsch (1951) has shown that very thin films of bismuth 
deposited from the vapour on to a surface at liquid helium temperatures 
exhibit superconductivity, with 7’, at about 5°K (although annealing at 
room temperature destroys this property of the film) and a large number of 
alloys containing bismuth have been found to be superconductors. 
The positive value of d7’,/dp found by Alexeyevski for the bismuth-rich 
compounds RhBi,, NiBi; and KBi,, together with the vanishingly 
small electronic specific heat of the pure metal at low temperatures 
led him to suggest that bismuth might be a ‘ virtual’ superconductor 
and that sufficient pressure might bring its transition temperature into the 
liquid helium region. Another interesting pointer arises from the curve of 
atomic volume of the elements plotted against atomic number; the 
superconducting metals are found to occupy a fairly well defined range of 
atomic volumes intermediate between the peaks and troughs of the curve 
(see for instance Mendelssohn 1952). Bismuth at ordinary pressures 
would appear to have just too great an atomic volume to take its place in 
this group. 

A number of runs were made under pressure with samples of bismuth of 
purity stated to exceed 99-:999°%. No signs of superconductivity could be 
observed (down to 2°K) at mean pressures less than 20 000 atm., but at all 
pressures tried between 20 000 atm. and 41 000 atm., a superconducting 
transition was observed at about 7°K, the value of 7’, not varying by 
more than 0-1°K throughout the range of pressures covered. This 
behaviour was reversible ; that is, after relaxation of the pressure at 
room temperature the specimen was again non-superconducting. 

It is of interest to consider this result in conjunction with the phase- 
equilibrium diagram for bismuth at high pressures determined by Bridgman 
(1935 b), in which polymorphic transitions are shown at room temperature 
at about 25 500 atm. and 27000 atm. In view of the uncertainty about 
the distribution of pressure in our experiments, it seems reasonable to 
associate the onset of superconductivity with the changes to more close- 
packed crystalline forms of bismuth discovered by Bridgman. 


Calcium and Strontium 


Strontium satisfies the criteria for superconductivity arising in the 
theories of Frohlich and Bardeen, and although calcium does not satisfy 
these criteria, both metals are exceptional in that their electrical resistance 
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at room temperature increases with pressure (Bridgman 1949, 1952)— 
a property which would be expected to favour the appearance of super- 
conductivity at high pressure. 

We have subjected samples of calcium of purity 97% and strontium 
of purity 99-85% to mean pressures of 44 000 and 42 000 atm. respectively, 
and examined them down to 2:1°x. No signs of superconductivity 
could be observed in either metal. 

The work is to be extended to other superconductors and non-super- 
conductors, and to higher pressures. 
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CXXXV. A Note on the Range-Scattering Method of Mass Estimation 
of Ionizing Particles 
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[Received August 4, 1953] 


SUMMARY 


The applicability of the method of Menon and Rochat to the estimation 
of the mass of singly-ionizing particles arrested in photographic emulsion 
is examined and modifications necessary to take account of (a) dip and 
(6) the variation of scattering constant with cell-size and velocity, are 
discussed. The ‘sliding cut-off’ method of eliminating large unrepre- 
sentative signals is described in detail. 


§1. InrRopUCTION 


THE range-scattering technique permits of very easy and rapid identi- 
fication of particles which happen to come to rest in the emulsion. The 
scattering method has the additional advantage of yielding estimates 
of mass which are independent of variations in the degree of development 
of the plates. 

In the method as described by Menon and Rochat (1951), the mass 
estimate depends on the experimental value of a statistical parameter P, 
which is a logarithmic function of the mass. /P is defined as 

P=1:37 log R+ 2:37 log «, 5 ages ae cies 
when RF is the residual range, in microns of the particle, and «. is the 
experimental value of the mean absolute angle of scattering per 100 
microns, measured along the faster half of the trajectory. For Ilford G5 
emulsions, the calibration measurements of Menon and Rochat using 
identified particles, yield the following expression. 

eras lOO (i eta. ae eo. (2) 
where yz is the mass of the particle expressed in Mev. 

This method has been used to estimate the masses of particles K1 
and K2 (O’Ceallaigh 1951) and has been employed in subsequent work on 
charged heavy mesons (Crussard et al. 1952, Menon and O’Ceallaigh 
1953, Friedlander, Harris and Menon 1953). 

Certain simplifying assumptions are made in using the parameter P 
to estimate mass. For instance, it can be shown that P depends on 
K*37, where K is the scattering constant, but in the method as it exists, 
the variation of K with'cell size and velocity has been neglected. These are 
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without significant influence in the application for which the method 
was originally conceived, but must be avoided in critical work such as 
mass estimates of particles of unknown mass, particularly when such 
measurements are carried out on long lengths of track. The method of 
taking into account such variations of the scattering constant, and a 
few modifications in the experimental procedure are considered in the 
following sections. 


§2. EFFECT OF THE VARIATION OF THE SCATTERING 
ConsTANT ON P 


Following Gottstein ef al. (1951), we may write for singly-charged 

particles (Oe Kp, 981? oH Kp, »Bs'? (3) 
$s pv p(B2—1) > . . . . . 

where (®), is mean absolute angle of scattering measured in cell-size s, 
s being expressed in units of 100u, p, v=momentum and velocity, 
B=(1—f?)-1?, »=mass of the particle in Mev and Kz ., expressed in 
(degrees x Mev)/(100 .)1/?, is the value of the scattering constant appro- 
priate to Bands. We have, by definition 


1 /& 
“es = —1/2 
% 100 eas ie (D).s dB, - + e e . (4) 


where B, and B, are the values of B corresponding to residual ranges 
Rand f/2 respectively. Hence 
P ce 1 ic Ks, Aa; dB 5 
100 u(B,—Bp) 5 Ll) (5) 
Since K;z ,) is not strongly dependent on B, it will be sufficient to regard 
it as constant and equal to K,z ,, the value corresponding to 3R/4. 
Integrating and neglecting powers of (B,—B,) greater than the first, 
we have with sufficient approximation 


a - 8 K s) 
%100u— ST ane 5 In| 3 =| ———— ae . ce . . (6) 
We assume as valid re a vee charged ne the range-velocity 
relation (B—1)=1-3975 x 10-*(Bju)°o, 
Substituting (7) in (6) and rearranging we have 
578 V8 8 
iy 22 43. 4K op oy oP oS * | (40 re ) . ° . . (8) 


If now we set K,p y=K, regarding it as a constant and independent of 
B and s, and also neglect the small corrective term Ky, »/4u which is 


- equivalent to setting B,+1 a 
Bile 
we have pore? — 43-4 K R~O578/e 5, . (9) 


which reduces to 
log p= 2°37 log (43-4 K)—1-37 log R—2:37log&. =. ~. (10) 
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Introducing the definition of P given in eqn. (1) we obtain 

P=2-37 log (48-4K)—logp. . . 2. . . . (11) 
Substituting K—=26-2 (degrees x Mev)/(100 pu) in (11) we obtain (2) 
which is Menon and Rochat’s expression. 

Recently, extensive calculations of stopping power for Ilford C2 nuclear 
research emulsions have been presented by Vigneron and Bogaardt (1951) 
and Vigneron (1953). We find that the results of Vigneron (1953) may 
be very well fitted by the range-velocity relationship 


B—1=1-460 x 10-%(Rju)j™58 SS (1) 
for 1:01<B<1-16. This would correspond to a modified value of (11) 
P=2:-32 log (41-4K)—logu. . . . . . (113) 


These results are in good agreement with those of experiment over the 
limited region hitherto examined (B<1-043), but there is a need for further 
experiments with faster particles. 

To deal with the case of G5 emulsions, and to allow for variation in 
moisture content of the emulsion, it would seem sufficient to modify the 
value of the constant in (12) choosing always the value 0-568 for the 
exponent. 

It is thus clear that the value of P depends on K287, and accurate 
estimates of mass require exact knowledge of the value of K. Further- 
more, its variation with s and B cannot be neglected. 

A theory of the variation of K with s and B has been given by Moliére 
(1947, 1948, 1951), and is discussed at length by Gottstein et al. (1951), to 
which reference may be made. 

It is shown there that the value of K,; ,) is determined by a parameter 
Q,, where Q,/s is dependent only on the velocity of the particle. Figure 2 
of Gottstein et al. gives 2,/s as a function of f?, and fig. 1 of the same paper 
gives K(, , as a function of 2,. The value K=26-2 is based on Menon 
and Rochat’s (1951) calibration point marked 5 in this figure, and corres- 
ponds to a value 2,~500, and thus to protons of residual range about 
1mm. Various calibration measurements that have been carried out in 
different laboratories suggest that the theoretical slope (e.g. shown in 
fig. 1 of Gottstein et al. 1951) may be taken with reasonable confidence as 
representing the variation of K.z ,, with 2,,. 

The method of calculating the value of « from (7) is therefore as follows. 
An approximate value of » is obtained from (1). This enables us to deter- 
mine (B—1) and hence, f? corresponding to 3R/4, by use of (7). The value 
of Q,/s is then found from fig. 2 of Gottstein ef ul. From fig. 1, we may 
obtain the appropriate K(y . for the value of 2,/s used in making the 
experimental estimate of &199,- The value of the small corrective term 
Kg, »/4 can be found from the approximate mass choosing K,, .=26. 
A corrected value of 1 can then be obtained by substitution of these values 
in (7). Usually it will not be necessary to repeat this procedure. 

The effect of neglecting the variation of Ki; ,) may be quite serious, 
and for particles of mass ~1000 m,, on long lengths of track, may possibly 
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be greater than the statistical uncertainties associated with the experi- 
mental determination of &199,,. 


§ 3. StanaL/Notse Ratio AnD Norse ELIMINATION 


In the application of the method by Menon and Rochat, the cell-size 
s was chosen so that the ratio, signal/noise=p, was at least 4. This 
requirement is probably over stringent, since with good instruments and 
experienced observers, the noise level tends to remain very constant. 
The optimum value of p can be shown to be ~2-6. Thus a value of p~2 
may well be used and the noise level determined for the conditions of the 
experiment. This may readily be done using the method of cross products 
suggested by Moliére (1951) the experimental application of which is 
described by Menon et al. (1951) Appendix I. Small accidental variations 
in noise level will still be without significant effect on the accuracy of the 
results. We may set up an equation for the minimum cell size s such that 
pe, € being any value fixed in advance. The solution is 
log s=3[0-578 log R+log {9./(p2—1)}—P/2:37] . . (13) 
where g=constant and ~0:15 micron for good instruments, skilled 
observers, and well lined up tracks. 
A plot of s on log-log paper will be linear, and may be used to find the 
desired value of s for selected p. The effect of changes in P on the choice 
of s is not great, and we may well use the same value of s for all particles. 


§ 4. CORRECTIONS FOR ‘ Dip’ 


It is necessary to allow for the effect of ‘ dip ’ on the measured values of 
Rand «. We require the value of these parameters in the unprocessed 
emulsion. Distinguishing these by the subscript w we have 

sure os | (14) 

“=a sec—/2 f.. 
R and « are the experimental values and are referred to the plane of the 
glass (x, y), the optic axis coinciding with the z axis. R,, and @,, are the 
values referred to a plane containing the mean trajectory of the particle, 
the normal to which makes an angle 0,, with the z axis. We assume that 
the effect of processing is a uniform linear translation in the direction of z 
and hence that tan @,—/S tan @ where S is the shrinkage factor. 


§ 5. ELIMINATION OF LARGE ANGLES 


It is found occasionally in multiple scattering measurements, that large 
unrepresentation values of second difference occur which would have little 
effect on the value of the mean if the sample were of sufficient size. They 
can, however, lead to unreliable estimates of the scattering parameter, 
when they occur in a small sample of second differences, obtained from a 
limited length of track. It is customary therefore, to employ a ‘ cut-off’ 
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procedure, to eliminate these large values of second difference. In the case 
of high energy particles, where the loss of energy along the track is 
negligible, it has been found convenient to employ what is known as a 
‘4D cut-off’ i.e. an angle is ‘ cut-off’ if it is greater than four times the 
mean of the second differences obtained without including it. 

In the Menon—Rochat method, this procedure is not permissible, since 
measurements are made on a section of track along which the energy of 
the particle, and hence the multiple scattering, is changing rapidly. In 
these circumstances, it is necessary to apply a ‘ sliding cut-off procedure ’ 
i.e., an angle «; is ‘ cut-off’ if it is greater than ¢,o,9 where &» is the mean 
angle obtained after excluding «,;.c; takes a sliding value, increasing 
continuously as one proceeds from FR to R/2 and has a value c,, at the 
range R.~—0-72 Rk, to which corresponds the measured mean angle of 
multiple scattering & (Menon and Rochat 1951). R is the true total 
range. Figure 1 of the present paper shows c, plotted as a function of range 


R, expressed as a fraction of the total range R. 


for a ‘4D sliding cut-off’. In this case c,,—4 and c; has a value 3:3 at 
R and 4:95 at R/2 and varies continuously in between. The calculation 
involved is based only on the energy loss by the particle along its trajectory. 

The application of the cut-off procedure is comparatively simple. The 
range R; at which the suspect value «; occurs is determined; R; may be 
expressed as a fraction of the total range R. The value of c; corresponding 
to R, is read off from the graph in fig. 1. «; can then be excluded if it 
exceeds 6,89 where & has been obtained without the inclusion of «;. 

It is clear, however, that the exclusion of an angle at a range near K 
would affect the mean to a lesser extent than the exclusion of an angle near 


4P2 


as é 
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R/2. This difficulty can be partly overcome by a substitution procedure. 
After an angle «; has been excluded, it is replaced by a value «;,x»), the angle 
to be expected at a range R; when the determined mean angle along the 
track is &. A new mean value & may then be obtained after inclusion of 
this value ex). In marginal cases, it is advisable to ‘verify that the 
excluded «, is greater than c,x, where « is the mean value finally obtained. 

The suitability of 3-5D, 4D, 4-5D and 5D sliding cut-offs, corresponding 
tO C.~=3-5, 4, 4:5, 5 respectively, are being tested in the calibration 
measurements of Friedlander et al. (1953). 


§6. CONCLUSIONS 


The variation of K,., is probably well described by curves of the type 
given in fig. 1 of Gottstein et al. The method described in this note is 
thus probably adequate to allow for the variations of K:, , with Q,. 
Some doubt still exists, however, as to the precise value of the scattering 
constant for a given value of 2,. This uncertainty may lead to serious 
errors in absolute mass determinations. 

In the calibration measurements of Menon and Rochat (1951) (point 
marked (5) in fig. 1 of Gottstein ef al.) only large angles which definitely 
appeared to be individual ‘ single scatters ’ were eliminated ; the procedure 
may thus be described as a ~8D cut-off. The position of point (5) is in 
better agreement with curve (a) which gives the theoretical values of the 
scattering constant for a ‘ no cut-off procedure ’. 

In the determinations of the mass of heavy charged mesons by Menonand 
O’Ceallaigh (1953), values of the scattering constant given by curve (a) 
of fig. 1 have been employed. We may remark that the estimates quoted 
may require revision on completion of the calibration measurements of 
Friedlander et al. (1953), which will give absolute values of the scattering 
constant for appropriate cut-off procedures. 


REFERENCES 


CrussaRD, J., MORELLET, D., TREMBLEY, J., and Orxkry-Lecourtors, A 
1952, C. R. Acad. Sci., Paris, 234, 84. 

FRIEDLANDER, M. W., Harris, G. G., and Menon, M. G. K., 1953, Proc. Roy. 
Soc. (in course of publication). 

FRIEDLANDER, M. W., Harris, G. G., Menon, M. G. K., and O’CEeatuaian, C., 
1953 (work in progress). 

GortsTEIN, K., Menon, M. G. K., Mutvny, J. H., O’CEatuaicn, C., and 
Rocuat, O., 1951, Phil. Mag., 42, 708. 

Menon, M. G. K., O’CEALLatGH, C., and Rocwart, O., 1951, Phil. Mag., 42, 932. 

Menon, M. G. K., and O’Cratuaian, C., 1958, Proc. Roy. Soc. (in course of 
publication). 

Menon, M. G. K., and Rocuart, O., 1951, Phil. Mag., 42, 1232. 

MoutmeE, G., 1947, Zeit. f. Naturforschung, 2A, 133; 1948, Ibid., 3A, 78; 1951, 
unpublished. 

O’CEALLAIGH, C., 1951, Phil. Mag., 42, 1032. 

VicNERON, L., 1953, J. de Phys. et Rad., 14, 145. 

VienrERON, L., and Bocaarpt, M., 1951, C.R. Acad. Sci., Paris, 233, 941. 


a! 


ar 


D. B. GAYTHER Phil. Mag, Ser. 7, Vol. 44, PI. 51, 


[ 1297 ] 


CXXXVI. The Decay-in-Flight of a Heavily-Ionizing r-Meson 


By D. B. GayTHER 
The Physical Laboratories, The University, Manchester* 


[Received August 29, 1953] 


ABSTRACT 


A r-meson has been observed to decay in the gas of a cloud chamber 
operated without a magnetic field. The t-meson was apparently created 
directly in an energetic nuclear interaction. Its energy at creation was 
about 100 Mev, and it lived for a proper time of (7-8-41-0) x 107° see in 
the illuminated volume of the chamber. 


SINCE the discovery of the 7-meson by Brown et. al. in 1949, several groups 
using photographic emulsions have established that this unstable particle 
decays into three charged particles which are probably z-mesons. The 
(-value, or energy release, lies in the range (70-75) Mev and the mass of 
the r-meson is about 975 m, (assuming the mass of the charged 7-meson to 
be 276 m,). These results have been obtained from measurements on a 
small number of decays at rest. our examples, however, of the decay of a 
z-meson-in-flight have been reported (Leighton and Wanlass 1952, 
Annis and Harmon 1952, Van Lint and Trilling 1953). In each case the 
7-meson decayed in the gas of a Wilson cloud chamber. The ionizations of 
these 7-mesons were indistinguishable from the minimum value. Recently, 
in an experiment on the Pic-du-Midi (2867 m), an event in the gas of a 
Wilson cloud chamber has been observed which is interpreted as the 
decay-in-flight of a heavily ionizing 7-meson. 

The cloud chamber has an illuminated volume of 40x 30x 20 cm and 
contains two lead plates ; the upper plate is 4 cm thick and the lower plate 
1-5cm thick. There is no heavy material above the chamber. An 
arrangement of enriched BF; proportional counters (see Barker ef. al. 
1952) is used to expand the chamber ; it selects nuclear interactions pro- 
duced in the upper lead plate. No magnetic field is used. The lower lead 
plate is used to examine the penetrating power of the shower particles. 

The event in the gas is shown in the plate (Plate 51) at level L. The 
heavily ionizing track A produces tracks 1, 2 and 3. Track A is copunctal, 
within errors of measurement, with the origin of the large interaction in 
the upper plate. Therefore, particle A must have been produced in this 
interaction and to have subsequently interacted or decayed in the gas at 
level L. No visible recoil track is produced and it therefore seems rather 
unlikely that particle A interacted with a gas nucleus. It is interesting to 
consider, therefore, whether the event could have been a decay process. 


* Communicated by Professor P. M. 8. Blackett, F.R.S. 
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The only known unstable particle which decays into three charged 
secondary particles is the 7-meson. 

The secondary tracks, 1, 2 and 3, could all have been produced by 
m-mesons. For example, track 1 has an estimated ionization of less than 
twice the minimum value but the particle failed to penetrate 1-6 cm. of 
lead, thus its mass was less than 400 m,. The measured angles between 
the various tracks and the estimated ionizations are given in rows (2), (3) 
and (4) of the table. ¢ is the angle between track A and a given track ; 
the error in each ¢ is about 1-5°. The angles between the projections of 
the secondaries on to the plane perpendicular to the line of flight of A 
are denoted by @;, etc. 


Summary of Data on the 7-Meson Decay 


(1) Track number A 1 2 3 
(2) ¢ (deg.) — 127-5 137 91-5 
(3) 6 (deg.) — 6,,=32 6,,=158 (057 =-110 
(4) Estimated ionization 

(< minimum) >8 <2 3-5 <2 
(5) Calculated ionization 

(< minimum) 12:5 2-1 4-6 1-8 
(6) Estimated momentum range 100-131 98-127 53-68 110-142 

(Mev/c) 
(7) Calculated momentum 

(Mev/c) 108413 105415 57+10 118+18 


Using the data in the table, the mass of particle A can be estimated in 
the following manner, assuming it to be a t-meson and the three secondary 
particles to be z-mesons. Particle 1 stopped in 1:6 cm of lead and there- 
fore had a momentum of less than 127 Mev/c. Assuming conservation of 
momentum at the decay and using the upper limit for the momentum of 
particle 1, upper limits for the momenta of particles A, 2 and 3, can be — 
obtained and are given in row (6) of the table. Track 3 is long and well 
illuminated and has an estimated ionization of less than twice the mini- 
mum value; its momentum is therefore greater than 110 Mev/e. Thus 
lower limits for the momenta of the other three particles can be obtained ; 
these values are also given in row (6) of the table. By using the momentum 
ranges for all the particles in the relativistic form of the conserva- 
tion of energy equation, it is found that the mass of the unstable particle 
lies within the range (960-1020) m,. It is thus consistent with the 
known mass of the 7-meson. The errors in each of the mass limits due to 
angle measurements are unlikely to be greater than about 15 m,. 

If it is now assumed that particle A was a 7-meson of mass 975 m,, then 
the momenta of the secondary particles and of particle A can be calculated. 
The values are given in row (7) of the table and the corresponding calcu- 
lated ionizations are in row (5). The agreement between the visual 
ionization estimates and the calculated values is very good, and it may be 
concluded that the event in the gas was produced by a t-meson. 
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The origin of the interaction with which the 7-meson is associated has 
been determined to within (1-2) mm by tracing back the tracks of several 
of the penetrating shower particles. The line-of-flight of the t-meson is 
copunctal with this origin. The path length of the 7-meson in the lead 
plate was (2-5+-0-2) cm and its energy at creation was (107-7) ev. It 
appears, therefore, that the 7-meson was produced directly in the high 
energy interaction and is probably not the secondary of a short-lived 
particle. This confirms the conclusion of Ceccarelli et. al. (1952) who 
observed a t-meson in a photographic emulsion which emerged directly 
from a star. The interaction in which the slow 7-meson was produced 
contained 10-15 penetrating particles most of which were probably mesons, 
thus the total energy of the interaction was probably greater than 20 Bev. 

The 7-meson lived for (7-8-E1-0) x 10-"° sec in the illuminated volume of 
the chamber and is the longest lived yet reported. The potential traversal 
time, or the time the 7-meson would have taken to traverse the space 
between the lead plates if it had not decayed, is (1-6-L0-2) x 10-® sec. 
These times have been measured from fiducial surfaces drawn just below 
the upper plate and just above the lower plate in the manner suggested 
by Wilson and Butler (1952). No direct estimate of the lifetime of the 
7-meson has yet been made but Herz (1953) has found a lower limit of 
10-° sec, based on observations of three decays at rest in photographic 
emulsions. 
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IN a paper on the “ Relativistic Hamiltonian System of a Particle and 
Relativistic Heisenberg’s Equation”, Tzou (1948) introduces the more 
symmetrical relativistic form of the Hamiltonian equations 
dxt# 0x dr, OM 
dr Cn, > dr sOxt 
and also the relativistic modification of Heisenberg’s equation} 


(u==1, 2; 3,4), «i atueney 


© 10, #1 SIA 00h eee 
where Q=Q (x", 7,): 


In §4 of Tzou’s paper the motion of a free electron is discussed (along 
similar lines to those of Dirac (1947 a)) and it is stated that the results have 
some deviation from those obtained by Dirac (1947 a) using the non- 
relativistic Heisenberg equation. The results, however, are equivalent as 
will be shown in this note. 

For a free electron ¥ has the operational form 


Hi —04f a1, 7) 7 
and the constant value He = = Ng ne oa lk 
This gives us the wave equation 

KH p=—ov( —m, ty) b= mgt. 5 ke oy A 
If we write H =—te Brrr; so late, ope an 
eqn. (4) becomes (1B"rr,,—myc)b=0 + ay ee cr 
where the f" are operators satisfying 

B* B’4-2" Bt 250"... a ee 
Then a =[x", #]=—ic PB". Voge SBS) een 


Therefore we see that 6" is the 4-vector (or world) velocity operator 
corresponding to the u-component of world velocity. Tzou, however, does 
not treat the time component of this vector, 64, as an operator, but as a 
scalar; and it is on account of this that the apparent deviation of his 
results from Dirac’s results appears. 


* Communicated by R. A. Newing. 


} Tzou (1948) writes this dQ/dr=[#, Q], and uses p, in place of z, for the 
4-vector momentum, 
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Equation (6) is the same as Tzou’s eqn. (18) and therefore leads to 

Tzou’s eqn. (31) for the world-velocity of the free electron 

dat he (dp 1) 

sas 1 se op -1 

aie 1, 14 — ae), exp ( ptr) H ee ed) 
Tzou now compares this result with Dirac’s result [which is for dx'/dt 
(j=1, 2, 3)], by assuming the classical relativistic relation 
dt/d7=[1—(v*/c?)]-¥/2 which enables one to write eqn. (8) as 


p2\ —1/2 dat ee 
(1- =) Teens 
p2\ —1/2 
and with p»=4, pt=—(1— =) : 


This makes f* a scalar, instead of an operator. It is clearly wrong, for [%, 
then, does not satisfy the necessary condition (§4)?—=1. 

The relation between dt and dr, appropriate to this problem, is obtained 
from eqn. (9) with w=4 (x4=ict) 


dt he (dp* A ) 
6 1 ¥f—1 
wa 14H 14 — (a =), exp ( 21 F7 1) c CRD) 


By treating §* as an operator we can reproduce Dirac’s result starting 
from eqn. (6). 


Multiplying eqn. (6) by —f* gives 
(—7P4B"2r,,+ BAmoc)p=0 


or (18) Btor,+-p,+-P4m,chp=0 (7=1,2,3) . . . - (11) 
since* T=(P1; Pa Ps, Up 4). 
If we write at—=19) 4, p=p* 


eqn. (11) is the same as Dirac’s wave-equation for the free electron 
(Dirac 1947 b) 
(a’p;+ pat p?m c)p=0 (j=1, 2, 3) 
(clearly the «’s and p? satisfy the necessary commutation relations). 
Put w~=4 in eqn. (8) 


_ at a4 
ic icB 
p4=—dt[dr 
and, remembering that (f*)?=1, we get 
d di\-ld 4,4 iO 
Wi == (=) A; =——8f a 5 e « * . ° ( 2) 
j 
and Bee pt Date (ose 
CY eee Mee a (1S) 


which agrees with Dirac’s result and shows that a=(«!, «?, «%) is the 
3-vector (or space) velocity operator. 


* In Tzou’s notation this would be 7,=(p1, Pa, Ps> Pa): 
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By eqn. (12) er 
j J 
ne ——inpst ees (4 eB 7): ren C)) 
d, 

Since in = pre _ A pl 
and PY H+ H pt=—2e rn, 
thdat _ 2 
then oy Tt at oe UCT 4 


cen, . . . . . . . e ° (15) 
for 74=ip,—tH/c, where H is Dirac’s relativistic Hamiltonian. 
Now z,, commutes with # so dz,/dr=0 and we get 


dn, ray EE =—=()). 

dt dr (16) 

AERC he opt ars —o 

dias tie aoe 

so that by differentiating eqn. (15) with respect to t 

ihdral da! 
s+ = > ery ein es Aa 
2 di* ~~ dt sea Sk Lu 


which is Dirac’s differential eqn. for «’. 

Equations (13) and (17) together reproduce Dirac’s results, thus showing 
that these results are not changed by using the relativistic Heisenberg 
equations. 
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Angular Distributions of Gamma Rays from 1°F (p, xy) 10 


By J. E. SanpErs* 
Cavendish Laboratory, Cambridge+ 


[Received August 31, 1953] 


MEASUREMENTS of the angular distributions of the 6-14 Mev (y,) and the 
combined 6-91 and 7-11 Mev (y.+-y3) radiations at the 874 kev resonance 
of the reaction 1°F (p, «y) *O have already been reported (Sanders 1952). 
Separation of the y, from the (y+) radiation was achieved by using a 
deuterium-filled gridded ion chamber as the gamma ray detector and 
analysing the photoproton pulses on a multi-channel kicksorter. 


* Communicated by the Author. 
+ Now at A.E.R.E., Harwell. 


Correspondence 1303 


Observations of the intensities of the y, and (y2+ys3) radiations were made 
at various angles 6 between the incident proton beam and the direction of 
emission of the radiation. 

Similar observations have now been made at the resonances at 935 and 
1381 kev. CaF, targets of between 10 and 20kev thickness were 
bombarded by the proton beam of the Cavendish electrostatic generator. 
If the results are fitted to expressions of the form I(9)=1--a cos? 6, the 


values found for the coefficient a are those shown in the accompanying 
table. 


Resonance a(y,) a(yo+ys) 
935 kev 0-00 +0-05 0-0 +0-1 
1381 kev —0-14+0-03 0-58 +0-13 


The results at 935 kev are in agreement with the isotropic angular distri- 
bution of the total gamma radiation at this resonance (Day et al. 1950) and 
with the assignment 1+ to the compound state of 2“Ne involved (Seed and 
French 1952). 

Combining the distributions at the 1381 kev resonance using the intensity 
ratio y4/(y2t+ys3) at 2=90° of 8-7 found in the present work, the angular 
distribution of the total radiation will be J(¢)=1—0-06 cos? 6, in good 
agreement with the measurements of Day et al. On the basis of the near- 
isotropy of the total radiation, it has been suggested (Chao 1950) that the 
compound state of 2°Ne involved is a 1+ state. The anisotropy of the 
component radiations revealed by the measurements described here makes 
this assignment unlikely, while a 2~ state would be consistent with the 
observed distributions. 
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The Velocities of Fundamental Particles 
By L. L. WuytTeE 


DirricuttiEs in the theory of particles may be connected with the role 
of time (Heisenberg 1951, Whyte 1952). The application to the various 
particles of the concept of velocity should therefore be re-examined. 

The presence of a linear velocity is a necessary inference only where two 
events have been observed which are localized and distinct both in space 
and in time: the departure and arrival of an entity which may be one 
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particle or the head of a pulse. This condition is not met in any of the 
experiments which purport to determine the velocities of particles lighter 
than nucleons. 

For example, all the ‘ determinations of electron velocity ’ (Thomson, 
Wiechert, Kirchner, Perry and Chaffee, Dunnington, etc.) used stationary 
conditions, i.e. continuous electron beams and static or periodic deflecting 
fields, and did not isolate two distinct events. Such experiments cannot 
discriminate between oscillation or vibration periods of an extended system 
and transit times of localized entities. The space and time coordinates of 
the departure and arrival of a particle and the presence of a transit time 
cannot be inferred from observations using static or periodic fields for they 
provide too little information, a fact inadequately recognized in the 
quantum-mechanical discussions. Moreover no electronic apparatus now 
in use provides compelling evidence for transit times varying with the 
voltage, for none meets the above condition. 

The only type of experiment which can justify the concept of electron 
velocity is one replacing static or periodic fields by separate signals, 
localized both in space and in time, of the departure and the arrival of a 
particle or of the head of a pulse. Such signals can be provided by the 
use of mechanical rotors, deflectors charged at successive moments by a 
single pulse, pairs of counters recording single particle lags, or similar 
arrangements. It may now be possible to apply such methods to electrons, 
positrons, and mesons, and it is of importance that this should be done. 

Experiments meeting the above condition have been carried out with 
neutrons and atomic or molecular beams, for example by using inertial 
motion in a rotating system, and the velocity concept confirmed. But this 
proves nothing regarding lighter particles, whose properties are in several 
respects different. 

Recent publications on Geiger counters have employed the phrases: 
“transit time of an electron”? (Stevenson 1952), ‘‘ over most of its path 
an electron is moving at nearly constant velocity”? (Ramsey 1952), 
though referring to overall time lags (relaxation times) of processes 
involving many electron paths and collisions. This looseness of expression 
obscures the fact that no measurements have yet been made which must 
be interpreted as transit times of electrons. 


I wish to thank Dr. D. Gabor for helpful discussions. 
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The Hxcitation Function for the Reaction *Be(p, p2n)’Be 


By D. C. Satrer and L. Brrp 
Atomic Energy Research Establishment, Harwell, Berks. 


[Received September 11, 1953] 


Tue excitation function for the reaction ®Be(p, p2n)7Be has been investi- 
gated for proton energies up to 155 Mev. The internal proton beam of 
the Harwell cyclotron was scattered by a thin tungsten target and 
refocused on to a ‘stacked foil’ target of beryllium and carbon plates 
(2 in. X 3 in. X ¥ in.) and brass absorbers, which was placed in a position 
diametrically opposite the tungsten target and below the median plane of 
the cyclotron magnet. This method of irradiation, first used by Hintz 
and Ramsey (1952), selects a narrow band of proton energies from the 
normal wide spectrum of the internal beam. The carbon plates were used 
to monitor the proton flux through the stack by means of the 
2C(p, 3p3n)7Be reaction, for which the excitation function is known 
(Dickson and Randle 1951). 

The energy spectrum from the 0-005 in. thick tungsten scatterer at a 
radius of 47% in. had been investigated previously (Dickson and Salter 
1953), and so the radial position of the ‘ stacked foil ’ target for maximum 
intensity of the scattered beam was already known. The upper and lower 
limits of energy of the protons entering the stack were 156-8 and 155-3 Mev 
respectively ; they were calculated from the radial positions of the scatterer 
and the targets (47, in. to 47+ in.), and the mean magnetic field along 
the scattering path. The range-energy curves of Aron, Hoffman and 
Williams (1949) were used to calculate the effective energies, and the 
spread of energies, of the protons traversing each plate. The range 
straggling of the protons was not an important source of energy spread. 

For the irradiation, the height of the main beam was adjusted to 4 in. 
using a ‘ chopper target ’, and the stack was set 2 in. below the median 
plane of the beam, so that no part of the beam could strike either the stack 
or its attachments. The stack was irradiated for 2+ hours and the active 
targets were loaded into identical sample holders, which completely 
enclosed them. 

The decay of each target was followed for fifty days using a scintillation 
counter with a sodium iodide crystal; at the end of this period the 
background count was less than 65% of the lowest target count observed. 
After a few days the logarithmic decay plots became straight lines, showing 
a half-life of approximately 53 days, characteristic of 7Be decay (Segre and 
Wiegand 1949). A source of 7Be, chemically extracted from carbon after 
cyclotron irradiation, was counted each day with the active targets to 
check the performance of the counter. Since the same active isotope, 
7Be, was being counted in the beryllium and carbon targets in identical 
geometry, the absolute efficiency of the counter was not needed. 
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The errors shown in the figure are relative, and arise from errors in the 
120,(p, 3p3n)7Be excitation function and statistical counting errors. In 
addition to these, there is an overall uncertainty of 8%, due to uncer- 
tainties in the !1C half-life, the 7Be branching ratio and the y-ray efficiency 
ratio of the Geiger counter used in the C(p, 3p3n)*Be work. 


oX 1027 (cm?) 
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0 40 80 120 160 
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Excitation function for the *Be(p, p2n)7Be reaction. The energy spectra for 
three beryllium targets are shown at the bottom of the graph. 


The curve is in excellent agreement with the one point previously 
reported (Randle, Dickson and Cassels 1951) and should be of use for the 
monitoring of high energy proton bombardments of beryllium targets. 


We should like to thank Mr. J. M. Dickson for his advice and interest 
throughout the experiment, and Dr. A. Blainey for providing the beryllium 
targets. 
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The Angular Distribution of the «-Particles from the “Li(p, «) Reaction 


By D. K. Cartwricut, L. L. Green and J. C. Witnmorr 
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[Received September 24, 1953] 


Tue angular distribution of the 7Li(p, «) reaction has been the subject of 
several studies in recent years (Martin et al. 1949, Talbott ef al. 1951, 
Hirst and Uebergang 1951). 

The angular distribution of the «-particles is of the form 

Y(0)= Y90(1-+ cos? 6+B cos? @) 
and the value of B at various energies is still rather uncertain. In the 
course of the setting up and checking of a piece of apparatus to study 
angular distributions we had occasion to investigate this reaction. 

The «-particles were detected in a proportional counter which subtended 
a solid angle of 0-15°% at the target, and the targets were less than 20 kev 
thick. The bombarding energy was electronically stabilized to 1 kv and 
the stabilizer was calibrated using the 1°F(p—y) resonances (Lauritsen 
et al. 1950). The statistical errors on each count were of the order of 
1% (about 20,000 counts in both the proportional counter, and a monitor 
counter at 90° to the beam). Counts were taken on both sides of the 
forward direction, and forward and backward at 90°. 

The results of A and B obtained by means of a‘ least squares ’ analysis, 
together with their probable errors, P, and Pp, are shown in the table 
for three different bombarding energies ; also shown are the results of 
Talbott et al. and Hirst and Uebergang in this energy region. 


a—Talbott, Busala and Weiffenbach (1951). 
b—Hirst and Uebergang (1951). 
c—Present authors. 


* Communicated by Professor H. W. B. Skinner. 
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Hirst and Uebergang do not give probable errors for B, but they should 
be of the same order as for A, and certainly not smaller. It will be seen 
that the probable errors of the results presented here are smaller by a 
factor of two to five than those of Hirst and Uebergang, who also 
used an electrical counting method. On the other hand, although the 
probable errors of Talbott et al. are only slightly larger than ours, these 
authors state that their errors are probably rather optimistic, as their 
results are compounded from several runs, the results of different individual 
runs showing wide variations. Hence the results presented here define 
the values of B, and of A, more closely, in this energy region, than 
heretofore. 


We are grateful to Professor H. W. B. Skinner for helpful discussions. 
One of us (D.K.C.) thanks the Department of Scientific and Industrial 
Research for a maintenance grant. 
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